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ABSTRACT 

The earth's magnetosphere is coupled to the lonosphere 
by iarge scale electric current systems. Knowledge of the 
geometry and characteristics ort these current systems can 
lead to an understanding of the solar-terrestrial 
interaction through which energy from tne interplanetary 
medium is transmitted into the magnetosphere and eventually 
dissipated. In this thesis is developed a comprehensive 
model of the magnetosphere-ionosphere current systems whose 
Character is consistent with present knowledge of 
ionospheric electric field and conductivity distributions. 
The model is designed to reproduce the average magnetic 
perturbation pattern at high latitudes as derived from 4 
study of 42 days of magnetometer array data recorded over 
three months in late 1971 and early 1972. The model is 
representative of winter conditions in the northern 
hemisphere. However, it is shown how the model current 
system must be modified to allow it to be used to describe 
magnetosphere-ionosphere coupling in the summer montis. The 
model studies show the importance of considering the 
Magnetic effects at the earth's surface of the Birkeland 
current sheets and the north-south “eS aNCAN Currents 
which link them. The major current systems in the model are 
the eastward and westward electrojets; both of these systems 
are fed by downward field-aligned current in the noon 


sector, wita nost OL Withe current ‘Jreturcning: “to ‘the 


iv 


| wth sottied 2unimee oes rok arinwoh stile nok wea 


Vv ( pe Ree Lai Dine _ a: f 
i, Beth : % 1 ; yn Wy ; 
ee eee Be!) rie aes Leh aaa 


a 


a f 
Fe 
9 ae ae 


Che ay hgh of Rie & eae o% ae re ie ae 
eer ee i % ls aa Pia y ef a Marana Tabi ; 


peaboaee aa iad "30>, tie oooh wile Apt ‘i shee’ + re act 
Lae stise tani <P ioe BS a & 


Kaen ene idamoed a Sad): weet ae 8. een : abi ei ; * 


sit aobliatapos. ie iano ‘atthe: al ahs 
pace aber « dup aN d sougavaatye oe ge avainan balk 
on uibietabana 2 Peg Abe 7 blogatiab: wae ay a | 
sendaga ae Ses aes pF oubmoR’ has ipitiat ds, gugthete sam 
stgenpel ws gysniss; aml ifs jaxyis je. ‘aneenenl et Ri win’ 


ono chia ‘ae. 2 Di wanted seit | ems alii 


ak Loko sit viewer gliee. Dns eke eight = 
arettace eas Gs 2 OL tak Debt dle 10 evs 
tnezwed Lele ae tou (Gms, Bot we een 
edyztees (Set en aah ibe oe) 33 Wotin ov Deteehes ou aoe 
. ; J : tah ic 
alg pare hha ta SOUR sIOGeL anf woe | se baie 


See 


Diels ae 


fs 
ore f § cay! a 
(Sehesti: vlgedgeodos Aynca-diaos oxy hire amet, th ea308 


ae fbn wg. as Sao LE Siwartue 1¢tee Sal othe anne - totes 
- eient ena’ weit ‘te ete bRobL as? yeis Laeissou tar troncasle ag 
aot. mee xt “Fass Birnipiis- btenh Bravavot 1k . 


j 


humisdas® ony to ssbpaye eds 7s bit 4a dteeRhd 


ait ee seincutes ereere? odo at 7 is dgitw 


ie | 1 t J 7 oe A - ie 


Magnetosphere along field lines in the pre-midnight 
quadrant. A small portion of the current in the pre-noon 
quadrant flows to low latitudes, nee the noon meridian, 
and returns to high latitudes in the dusk sector. A small 
portion of the electrojet current diverges up magnetic field 
lines at the dawn and dusk terminators which separate the 
sunlit and dark: ionospheres. In-this context, it is shown 
that, at quiet times, the electrojet resulting from the 
ionospheric electric field being superposed on the sunlit 
ionosphere (vize, the UV electrojet) produces significant 
Magnetic perturbations at the earth's surface. Finally, it 
is shown that the model can be perturbed in the evening 
sector, in a Manner consistent with known substor, 
phenomenology, to produce magnetic perturbation patterns in 


excellent agreement with observations of polar itagqnetic 


Ssubstorms. 
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CHAPTER 1 INTRODUCTION 


1.1 The Magnetosphere 

Early in this century it was believed that the earth 
WaS surrounded by the vacuum of interplanetary space, and 
that the magnetic field was normally that of a pure dipole. 
However, the work of Birkeland (1908,1913) led to the 
suggestion that solar flares caused bursts of solar gas _ to 
be emitted by the sun and that, at these times, the earth's 
piagetic field would carve out a cavity in the streaming 
solar gas. Chapman and Ferraro (1931) considered the 
interaction of the solar gas with the earth's magnetic field 
configuration under the circumstances associated with solar 
flares, and they were able to estimate the position of the 
boundary between the earth's magnetic field cavity and the 
interplanetary medium. Biermann (1951), noting that comet 
tails were always directed away from the sun, independent of 
the direction of motion of the comet, inferred that there 
waS a continual streaming of matter away from the sun. 
Parker (1958,1963) developed extensive theoretical arguments 
in which it was assumed that a continuous solar wind existed 
and, with the advent of space probes, it has been proven 
that interplanetary space is filled with a plasma streaming 
away from the sun. Due to this continuous presence of the 
solar wind, it became apparent that the magnetic cavity, 
termed the magnetosphere, is a permanent feature of the 


€arth's environment. 
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A schematic drawing illustrating most of the features 
of the magnetosphere is shown in Figure 1.1. It is apparent 
that present knowledge has lead to the realization that the 
magnetosphere iS a complicated configuration of particles 
and fields. It is not the purpose of eee thesis to provide 
a detailed understanding of the entire magnetosphere, but a 
brief description of it is in order. 

The earth with its magnetic field acts as an obstacle 
to the flow of the solar wind. Average solar wind bulk 
velocities are of the order of 300 to 500 km s-t, with 
number densities of 5 to 10 cm-1!, and temperatures in the 
range Of 710% “OK to 10% °K. These vaiues give’ an Aliven 
velocity of 50 to 100 km s~1!, and a sound velocity of 100 to 
200. km s~t so that the solar wind is thus both supersonic 
and super-Alfvenic. Consequently, a standing 
mgnetohydrodynamic (MHD) shock wave is set up, typically 
about 14 earth radii (Re) upstream from the earth, as shown 
in Figure 1.1. Immediately interior to this is the region 
cailed the magnetosheath, which is a region of thermalized 
plasma. The surface aiong which the solar wind plasma 
pressure is balanced by the earth's magnetic field pressure 
is called the magnetopause. 

The solar wind biowing past this cavity distorts it 
such that it is blunt on the Sunward side of the earth, and 
extended in the anti-solar direction to form the so-calied 
magnetotail. It is now known that the magnetotail extends 


beyond the orbit of the moon {about 60 Re) (Ness et al, 
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Figure 1.1 A schematic drawing of the magnetosphere 
(after Heikkila, 1972) modified by Rostoker 
to include the plasma mantle. 
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1967), although its)-actual length: is not known at’ the 
present time. Axford and Hines (1961) proposed the idea of 
plasma convection within the magnetosphere, but their 
magnetosphere was a closed, teair-drop shaped cavity. They 
Suggested a "viscous-like" interaction between the solar 
wind and the magnetospheric plasma in the neighbourhood of 
the magnetopause. This interaction would result in plasma 
flowing down-tail, and eventually returning along the sun- 
earth line to give a closed circulation pattern. This motion 
is shown in Figure 1.2. Close to the earth, there exists a 
region of stably trapped particles, the plasmasphere, which 
is bounded by the plasmapause. Beyond the plasmapause, on 
the nightside of the earth, there is a region of plasma 
symmetrically distributed on either side of a mid-plane, and 
extending along the length of the tail (see Figure 1.1). On 
average, this so-called plasma sheet has a thickness of 
about 5 Re at the center, and widens to about double this at 
the flanks (near the magnetopause). The magnetic field along 
the central plane of the plasma sheet is weak and normal to 
this plane. This region of weak field is called the neutral 
sheet at higher latitudes. Outside the plasma sheet are the 
tail lobes, in which the plasma density is very low (<C.01 
cm-3) compared to densities of about 0.1 to 10 cm~3 in the 
plasma sheet itself. 

Figure Veg depicts convection mapped onto the 
equatorial plane. However, if this motion is mapped to 


ionospheric altitudes, it corresponds to anti-sunward 
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Figure 1.2 Streamlines of plasma convection flow in the 
equatorial plane of the magnetosphere (After 
Axford and Hines, 1961). 
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Gonvection. across . the polar. cap. regions, and sunward 
convection at lower latitudes. In a collisionless plasma, 
both ions and electrons would follow the convection pattern. 
However, in the lower ionosphere, ions undergo many 
GOLltisions “With neutral atoms, so. that the ions .remain 
almost stationary in comparison with the electrons. This 
Pazeterecntial motion. of, electrons..and ions, produces Hall 
current in the direction opposite to that of the convecting 
plasma. The convection pattern of Axford and Hines’ then 
corresponds approximately to westward flowing Hall current 
in the pre-noon sector, at sub polar cap latitudes, . and 
eastward flowing current in the post-noon sector. This basic 
convection pattern is an oft-repeated observation (see the 
review by Cauffman and Gurnett (1972), and the references 
therein) and the currents associated with it correspond in 
Hauge *pagt.to those to be discussed later in;this thesis, 
Although the basic convection pattern derived by Axford 
aide Hines. (j1961)-shas withstood the test of extensive 
observations, the concept of the closed magnetosphere has 
not. Observations of solar electrons and low energy solar 
protons (identified by their energy spectra) at high 
latitudes, indicate an essentially uniform flux of these 
particles (Stone, 1964; Fennel, 19/33 Vampola, wasi> 
PeD2AruAdnand BUCEOWS, wl 9/0. yeanhongs others) .. it. hasysbeen 
concluded that, since solar particles have free access to 
Ehevhagh latitude regions ..of. the, .earth,. that .the high 


latitude field lines arepiairectlys..conneeted. to «the 
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Pocerplanetary | Magnetic “field (IMF). I+ is therefore 
believed that the magnetosphere is "open". Dungey (1961) 
proposed magnetic field line merging as a means by which the 
magnetosphere could become open. Figure 1.3 is a schematic 
depicting this process. The IMF merges with the earth's 
dipolar field at the magnetopause, and the merged field 
dines then are-convected in the anti-sunward direction by 
mieomsOLdi Wilt. [hiss COnvection: /of field... lines occurs 
because, in the solar wind, field lines are "frozen-in" to 
the plasma. This convective motion results in field lines 
moving across the polar cap pen the dayside to the 
nightside,and returning at sub-polar cap latitudes. 
Recently, Rosenbauer et al, (1975) have suggested that solar 
wind plasma may enter the magnetosphere directly in the 
polar cleft (Figure 1.1) to form the plasma mantle streaming 
away from the earth. The existen e of the plasma also 
indicates a direct coupling of the interplanetary mediun 


with the magnetosphere. 


1.2 Magnetosphere-Ionosphere Interactions 

Particles in the plasma are subject to drift motions 
due to electric fields, pressure gradients (Siscoe, 1966), 
magnetic field gradients and magnetic field curvature (see, 
for example, Roederer, 1970). Some of these drifting 
particles precipitate into the upper atmosphere, where due 


to ionization and excitation of neutral atoms and molecules, 


a region of enhanced conductivity is created. It is in this 
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Noon-midnight meridian section of the open 
magnetosphere showing field line merging at 
the nose of the magnetopause. The numbers 
represent the motion of field lines, with 
motion proceeding to the higher numbers. 
Note that in this figure, the north pole has 
been placed at the bottom of the figure. 
(After Hess, 1968), 


Ah 


UOI}OOIID 
PUuIM IelOS 


~~ F 
ad! 
rye 
ne 
of 
ada % 
$22 wher 
: 
whit Pot 
hee . 
he a 


aor 


= 4 
EMnge a 
ge 
+ ’ 
A Ee 


ey 


twalter gent a a 


ky A 


any 


S 


ee 


_ 
* 


12 


region that the major horizontal ionospheric currents flow. 
If it is assumed that the conductivity along field lines is 
sufficiently high to permit neglect of any potential drops 
along field lines, and further, that the field lines are 


vertical, then the equation 
Bre 
Ee DO) eee 1a, 


states a generalized Ohm's Law for the ionosphere. Here, 
P 


conductivities («see below), Ens is the electric field 


and 3 are height-integrated Pedersen and Hall 


perpendicular to the magnetic ‘field, 8 ; and is the 
height-integrated current density. The height-integrated 
conductivities are derived from the height dependent 
conductivities, which, for the ionosphere, and assuming a 


Single ion type of mass M, are given by 


nett /me net; /r9 


. se . a 
On oat a ae 7 PS gis (ste 
and 
Gn ean net a a: 
~~ = pe 


1 ta) t,7 Reger: 


where (4 and @; are the electron and ion collision periods 
with neutral atoms, Sl is the ion gyrofrequency, and @ the 
electron gyrofrequency. Figure 1.4 shows the variation of 
On and OF with altitude, It is clear that the Pedersen 


conductivity predominates at high altitudes, and this fact 
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Figure 1.4(a) 


Figure 1.4 (b) 
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will be employed in Chapter 4. In the ionospheric E region 
(altitude range 85 to 140 km), electrons may gyrate freely, 
(437, > 1). but ions are impeded due to collisions, so that 
2TS 1.6 From equations 1.2 and 1.3, the Pedersen current is 
Carried mainly by ions, the Hall aun by -electrons. “At 
higher altitudes, around 160 km (the F region), both ions 
and electrons gyrate and drift essentially freely, so that 
the two terms in Or, (equation 1.2) almost cancel. However, 
there remains a significant Pedersen sapere ee 

As well as the horizontal currents which flow by virtue 
Sethe horizontal ~electric fields and Hall and Pedersen 
conductivities, field-aligned currents also flow. Birkeland 
(1908,1913) first introduced the concept of field-aligned 
currents, but this work did not gain immediate acceptance. 
Fejer (1961) and Kern (1962) reintroduced the concept, and 
Bostrom (1964) made extensive use of field-aligned currents 
in his magnetospheric substorm models. 

That, in reality, field-aligned currents must flow, has 
been shown by Vasyliunas (1968). Consider the non-height- 
integrated form of equation 1.1, including field-aligned 
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electrostatic, E=-Vo ,where ri) is the electrostatic 


potential. Therefore, 
E-J = -V-(¢J) 1.6 


Integration of this over the volume of the ionosphere under 


consideration, and noting from equation 1.5 that E.J > 0, 


fd 


Jerd = -[o-@d)dv =-[bJ-da > 0 


The final integral is over the surface of the system, and 


1.7 


Since it is non-zero, indicates that currents are flowing 
into and out of the systen. 

Field-aligned currents were accepted before direct 
experimental evidence through in situ measurements was 
forthcoming. zZmuda et al (1967) and Armstrong and 
Zmuda (1970) presented evidence from polar Orbiting 
satellites which indicated that field-aligned current sheets 
flowed into and out of the auroral regions. These results 
were based on the existence of level shifts in the east-west 
component of the measured magnetic field. Many similar 
studies of data from polar-orbiting satellites now exist 
(Zmuda and Armstrong, 1974; Sugiura and Potemra, 1976; 
Iijima and Potemra, 1976), and have confirmed that indeed, 
large scale field-aligned currents connect the outer 
Magnetosphere to the ionosphere. Anderson and Vondrak (1975) 


have reviewed the data concerning field-aligned currents. 
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1.3 Large-scale Current Systems 
1.3.1 The Sq and L current Systems 

Observations of geomagnetic disturbances have led to 
the description of a number of large-scale ionospheric 
current systems. AS global magnetic variations became 
availabie, it became possible to map the quiet-day 
variations in decai, and a solar Variation (Sq); was 
isolated. Basically, this variation is caused by ionospheric 
currents which are driven by an atmospheric dynamo. E-region 


plasma is wind-driven with a horizontal velocity V across 


ar 
frecmemagnettC- ileld, B,) iresuiting in the production of an 
peer etc LCL d WIV xe Neon eld draves a height=integrated 
current, I. As well \av-uspace. charge, electric field, Ez, 


develops to give a total electric field, — + ¥V x B. The form 


of the Sq current system can be found from the solution of 


7s = 10 


where S- is the height-integrated conductivity tensor. (See 


aa 
Aas 


for example, Akasofu-and Chapman, 1972). Figure 1.5 is a 
plot of the Sq current systems for December 
solstitial,equinoctial, and June solstitial seasons, as well 
as the yearly average during IGY (1958) (Matsushita, 1969). 
The current direction is counter-clockwise in the northern 


hemisphere, in these figures. 
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Another geomagnetic fluctuation present at all times is 
that found to be associated with atmospheric tidal notions, 
driven by the moon's gravity field. The perturbations due to 
this so-called L current system are approximately one order 


of magnitude less than the Sq variations. 


1.3.2 The Dst and _DP1 Current Systems 

When the Sq and L magnetic variations are removed from 
ene data, /the so-called.disturbance field (D)\remains. This 
field can be represented as the sum of two components, Dst 
and DS. The Dst component is independent of longitude, and 
DS is the difference between D and Dst. The Dst component is 
believed to be due to a westward flowing ring current in the 
equatorial plane and several earth radii (Re) from the 
earth. Many equivalent current systems consistent with the 
DS component have been put forth in the literature. The DPi 
system (Obayashi, 1967) is the one associated with magnetic 
SuustOmis. Wigure | 120 | 1S seamwpOLlar DLOt of jtnisvclassical 
system. Note that relatively intense current flows both 
eastward and westward near 70°N latitude, corresponding to 


the location of the auroral oval. 


1.3.3 The a and DP=2 Current. Systems 
DULANGEti Mes OL very Low magnetic activity, there is 4 
polar cap magnetic variation which has been called Sf” by 


Nagata and Kokubun (1962). Subsequent authors (Kawasaki and 


Akasofu, 1967; Feldstein and Zaitsev, 1967) revised this 
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Figure 1.6 
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An equivalent ionospheric current system. If 
the current between contours is 105 A, this 
System would reproduce the geomagnetic 
disturbance (ane the DP1 type. (After 
Piddington, 1969). 
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System tO exclude any possibility of magnetic disturbance. 
Figure 1.7, from Kawasaki and Akasofu (1967), shows the 
horizontal magnetic disturbance vectors for the Ss system 
for May 8,1964. In the same paper, these authors verified 
Nagata and Kokubun's observation that x was much stronger 
in the summer months, ee concluded that the sf magnetic 
Variation was due to an ionospheric current systen. 
Subsequent work (Kawasaki and Akasofu, 1973) has shown that 
the data of Figure 1.7 may be reasonably modelled by a 
combination of field-aligned and ionospheric currents. 
Nishida (1968) proposed an additional current systen, 
DP-2 (Figure 1.8), but Akasofu et al (1973) have expressed 
the belief that the DP=-2 current system is in fact due to a 
combined effect of expansion of the auroral oval and an 


enhancement the Sy system. 


1.3.4 The Electrojets and the Interplanetary Magnetic Field 
Harang (1946) published observations of the diurnal 
variation of magnetic disturbance vectors at high latitudes 
(53. 9°N $to 7a. 2°No Geomagnetic. Latatude) «it 1s clear from 
his work and ali subsequent work that there are two main 
ionospheric equivalent currents, the eastward and westward 
electrojets, which can account for the majority of magnetic 
observations at hignhyelatitudes. The eastward current 
produces a positive (i.ee., northward) north-south component 
of magnetic field, the westward current a negative north- 


south component, and these signatures are apparent in 
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Figure 1./7 A polar plot of the horizontal geomagnetic 
onent at high latitudes on the extremely 
quiet day May 8, 1964. (After Kawasaki and 

Ak £ 7). 
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Harang's publishea data, witheastward current flowing fron 
near noon to approximately Rees ie and westward current 
flowing from nocn into the pre-midnight sector. The region 
of current reversal, in the pre-midnight sector, has become 
known as the Harang discontinuity. 

It has become clear in the past few years that the 
magnitude of auroral zone magnetic perturbations is governed 
py. the interplanetary magnetic field (IMF) and an 


partacular, by the sign of ther. vernticad. jue and 


z) 
azimuthal (By) components of the IMF. As well, the average 
Magnitude of the westward and eastward electrojets decreases 
from Summer to winter (Meng and Akasofu, 1968; Langel and 
Brown, 1974, Friis-Christensen and Wilhjelnm, iPS bey ke This 
seasonal effect is not the same for both electrojets when 
examined in terms of the role of the IMF, as will be seen 
below. 

iD etehNsetOtiatieldglineymerging, at vsynotssurprising 
that. the IMF vertical. component  has:»a omarked effect,;.on 
ionospheric currents, Since the merging rate is enhanced 
when B,<0. Akasofu (1977) has given an excellent review of 
the relationships between By and various ionospheric 
phenomena. He also points out eee gEeate caressa mustiibe 
exercised in interpreting scme of the relationships which 
have been claimed to exist. For example,the AE index,which 
is a measure of world-wide auroral activity, has been 


correlated with Bp. The AE index is computed from magnetic 


records from a number of auroral zone magnetic observatories 
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which are approximately uniformiy distributed in longitude 
(Baivis hhandeqsugiurae1966¢-4 fAllens+and eKsoéhl, 1975). The 
horizontal component records are superposed, and the curves 
defining the upper and lower envelopes of these superposed 
records are defined as AU and AL respectively. Then the AE 
index is simply the distance between AU and AL. i.e., 
AE = ({AU{ + JAL{ 

The AE index presupposes that at least one observatory can 
monitor the westward electrojet to give AL, and at least one 
observatory can monitor the eastward electrojet, to give AU. 
It is assumed that these electrojets are most intense in the 
auroral zone, where these observatories are located. 
However, Holzworth and Meng (1975) have shown that Bo 
affects the size of the auroral oval, so that when Bb, is 
northward, the radius of the auroral oval is reduced. Thus, 
Akasofu points out, a strong correlation between B, and AE 
may be, at least in part, only due to the correiaticn 
between B, and the size of the auroral: oval. Hirshberg and 
Colburn (1969) have done a similar analysis between Kp anda 


B Kp is a 3-hour index of the level of woridwide 


2° 
geomagnetic activity introduced by Bartels (1949) (see 
review by Rostoker, 1972). This index is computed ina 
complicated manner from data from sub~auroral zone 
observatories, and Akasofu's criticism may also apply. 
However, Hirshberg and Colburn (1969) showed that there was 


a good correlation between Kp and the magnitude of Boe when 


BL <0. In any case, negative B_ causes an enhanced merging of 
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field lines with a resulting increase in convection rate and 
thus enhanced electrojet current strengths. 

The east-west component of the IMF (By, ) aiso plays >a 
regulatory role in: ionospheric current dynamics. Solar 
magnetic field lines are carried away from the sun by the 
solar wind, and due to the rotation of the sun, take on an 
Archimedes spiral configuration. Wilcox and Ness (1965) 
detected the interplanetary sector structure in which, on 
average, four sectors of alternating magnetic field 
direction (away from the sun, and towards the sun), were 
found for every solar rotation. When the IMF is directed 
away from the sun, it corresponds to a positive By s 
svalgaard (1968,1973), Mansurov (1969), and Friis- 
Christensen et al (1972) have shown that the By~component 
affects the polar cap magnetic field. In particular, when 
By? 0, the inferred polar cap equivalent current is 
‘counterclockwise (looking down on the north pole) and the 
reverse direction when By <0. However, during winter months, 
the magnitude of the polar cap magnetic variations are 
extremely smail, approximately two orders of magnitude less 
than during summer months (Svalgaard, 1973). 

Returning to the seasonal differences in the eastward 


and westward electrojets and the role of B,, Langel and 


Brown (1974) have found that, during the summer, iAH} and 


{AZ{ in the westward jet region are larger than in winter 
when By <0. During away sectors (By 0). (AVits Bsevléess Gin 


summer than in winter, while {AZ{ is approximately the same 
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in winter and in summer. In the eastward electrojet region, 
{AVH{ and Zot are both greater in the summer. 

The observations concerning the role.of the IMF upon 
ionospheric current systems are presented here for 
completeness. In this thesis, data have not been ordered 
according to the polarity of the IMF. However, an argument 
is developed in Chapter 4 which depends, in part, on the 


above considerations. 


1.4 Ionospheric Current Modeis 

The uitimate test of any interpretation of data is the 
construction of a model consistent with that interpretation, 
and the comparison of computed results with the data. A 
large number of current models have been devised to describe 
specific current systems, although many of these have 
consisted ~of, “only ao0onospheric currents. Kawasaki and 
Akasofu (1973) had reasonable success in modelling the Sy 
current system, uSing a model consisting of fielid-aligned 
current (along (dipote?  freid™ tines) s cand “a ~flat= eGarta 
approximation for the ionospheric currents. Maeda and 
Maekawa (1973) undertook a numerical study of modelling 
Giopal polar Lonospheric “currents. They “considered that 
current would be driven by i) an atmospheric dynamo electric 
tread. i) +a. polarization. electric “fieid,” arising ~ from 
Conductivity gradients; “and 1739" an externally applied 
Srectric) f1e1d,, 100cuced by, the solar wind. They did not 


include the effect of field-aligned currents, although the 
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existence of such currents as sources and sinks was pointed 
out. Three different models of conductivity were employed, 
as described by Fejer (1953). Maeda and Maekawa found that 
the geomagnetic variations produced by one type of 
Polal@ZatjiOn 2 electrics iiveld, in | Which a. /positive,point 
charge was placed near dawn, and a negative point charge 
near dusk, are very similar to Sy. tase Charge \distreibutaon 
is Similar to that used by Kawasaki and Akasofu (1973), 
although these authors used charge sheets. However, Maeda 
and Maekawa do not give a detailed comparison between their 
calculations and real data. 

Several globai current models have been developed by 
Posibera, .et oli ati 27 5o\en Iny their investigation,,itiwas 
assumed that field-aligned current sheets flowed into and 
DprOoristhes aulopaluotvcoleateite borders usuch»thatucurnen: 
was downwards on the poleward side, from midnight to noon, 
and upwards on the poleward side of the oval from noon to 
Midnight. The field-aligned current at the equatorward 
borders was of the opposite sense, and of reduced intensity 
by a factor of approximately 2. For the model developed for 
relatively quiet periods, the auroral oval was assumed to be 
an annulus with an enhanced conductivity structure with 
respect to the region outside the annulus. The ratio of 
Pedersen to Hall height-integrated conductivities was two 
everywhere. To determine the ionospheric current flow, the 


equation 


V-L = J) sin je 
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was solved, assuming that the field-aligned current, iy 6 
was vertical, (i.e., x =90°), and where I is the height- 
integrated ionospheric current density. This model yields 
both eastward and westward electrojets confined to the 
auroral oval, as shown in Figure 1.9 . AS well, some current 
flows across the polar cap, from about midnight to noon, and 
there is weak current at sub-auroral latitudes. However, the 
choice of field-aligned current distribution was somewhat 
acuttrary 20 that thesexact distribution was | not «known at 
that time. As a result, the currents derived from this model 
do not appear too realistic, inasmuch as eastward current is 
showh just past dawn, and magnetometer evidence does not 
verify this (see Chapter 3 of this thesis for example). 
Although this model was not tested against ground-based 
magnetic data, it does demonstrate an approach that could be 
quite fruitful in approaching ionospheric current flow. As 
‘well, these authors demonstrated for the first time that 
east-west aligned field-aligned current sheet pairs, if 
unbalanced along a meridian, can produce ionospheric current 
flows that resemble those which are believed to exist. 
Nothing has been said concerning the mapping of 
these various current systems to the outer magnetosphere. AS 
this mapping is required to understand the possible driving 
mechanisms for the currents, it will be included in Chapter 
5 of this thesis, where some possible current generators are 


described, 
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FACGULne 4.9 A model ionospheric current pattern. The 
current) strength between two adjacent 
current contours LS 2h ay re (After 
Yaestinaravetiai 19756) 
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1.5 The Objectives of this Thesis 

The data analysis presented in the thesis was 
undertaken in the hope that it would lead to a model of 
ionospheric and magnetospheric currents which would be 
consistent with all known observations in the high latitude 
regions of the earth. This goal was, in retrospect, perhaps 
too ambitious. However, working from ground-based magnetic 
data, and using average characteristics of the ionospheric 
electric. field and Zonospheric, ,conductivities,..a current 
model has been developed, which, although it does not 
reproduce the observations exactly, does reproduce in a 
Statistically significant, way, (many features of the data. 
Phaiseisethesfirstatime that. a. current. model» has been..so 
rigidly tested against real data, and the prime value of 


this thesis lies in that test. 
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CHAPTER 2 DATA PROCESSING AND ANALYSIS 


2.1 Data Handling 

From late in 1969 to early in 1972, the University of 
Alberta intermittently operated a number of 3-component 
magnetic observatories. In the centered dipole systen of 
coordinates, these observatories all lie along approximately 
the 300°E meridian. Table 2.1 gives the positions of these 
observatories, as well as the observatory or station 
fMNMemonic, and Figure 2.1 shows a map of the location of 
these stations. As well as these observatories, data from 
Resolute Bay (RESO; 33..0°N, 2906 Oe) and Newport 
(NEWP355.1°N,3C00.0°F) were available to extend the coverage 
of the line of stations both northward and southward. At 
each of the sites shown in Figure 2.1, data were recorded 
digitally on 7-track magnetic tape, with a sample rate of 
‘one data point in each of 3 channels every 1.92 seconds, 
timing Of) the date ‘being ‘considered accurate to s41s 
MEdsabech, 4972). The system has a dynamic range of +1000nT 
with a +1nT sensitivity. As well as recording the magnetic 
data, a timing signal from radio station WWVB was recorded 
for 2 minutes every 7.5 hours. Kisabeth (1972) has described 
the instrumentation in some detail. 

For computer compatibility, thej| field tapes were 
converted to 9-track station master tapes, containing data 
for one station only. As this was being done, the WWVB 


timing data were interpreted by a special computer code, and 
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Figure 2.1 


A map showing the location of the University 
of Alberta magnetic observatory sites. Data 
were also available from Newport, UsS.A., 
and Resolute Bay, Canada (not shown) along 
the -Same meridian “as. the University =e. 
Alberta observatories. 
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PenenGg ein tornation, salong ewithe the ~station, mnenonic. was 
Gewehenn aSued@ labele forayeachaeblocks/ofvsadata (1 hour, 5 
Minutes, 32 seconds of data). In some instances, the WWVB 
timing information was not machine decodeable. Then, a Cal 
Comps plot was’ made ,of.,thee,WNVB. ¢+signal,/.and. it was 
interpreted by hand. 

When all data had been transferred to 9-track magnetic 
tapes, a so-called event tape was created for specific 
intervais of time. Each event tape contains a series of 
files, with each file containing alimwdacaywetoOm they tame 
interval for one station. Thus, each event tape contains the 
data for all the stations for a given interval of time. The 
availability of event tapes minimizes the number of tapes 
Which must be manipulated to study a specific period of 
time, or event. 

Atethe; beginningjot(thisrthesis study,” the) «data. were 
available on 9-track master tapes; data acquisition did not 
playeets. oni¢gicant, rolesing thisethesis. 

The magnetometer sensors at each of the observatories 
were aligned in the direction of local magnetic north, local 
Magnetic east, and the vertical (H,D, and Z respectively, 
where H is positive northward, D is positive eastward, and 
Z is. positive downward). To facilitate the interpretation 
of the data and the modelling of the magnetic field sources, 
the centered dipole coordinate system (an orthogonal system) 
was used, aaa all data in this thesis will be presented in 


this system. This required that the measured horizontal 
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components ( H, D) of the magnetic field be rotated into the 
centered dipole system through the transformation 

xt cos@ -=sin® H 

yt sin 0 cos O/7>\D 
where the primed coordinates represent the centered dipole 
system and the unprimed coordinates the local magnetic 
system. Here, G is defined as the difference between the 
hagnetic declination OL the station and the dipole 
declination. Table 2.2 lists 9 for the observatories used in 


this thesis. 
2.2 Data Processing 


2.201 Latitude Profiles 

The primary suite ofedata that wassused in this study 
was from the time period from Day 332,1971, (November 28, 
moi) to Day 24, 199 2,eWanuenye24, 1972) The data wiil be 
presented in the form of latitude profiles. A latitude 
profile ic constructed from) data’ points “acquired ‘at all 
Statcons vate oa wediven “tame, | plotted as: a -runction of the 
latitude of the station. This method of presentation has 
been chosen for several reasons. A latitude profile depicts 
aelargesanount Of, information in ai singie frame; that Js, 
the latitudinal behaviour of the three components of the 
Magnetic field at a given time is shown on one plot. Also, 
guemany cases, a latitude profile can be interpreted readily 


in terms of overhead ionospheric currents {see below). These 
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two features greatly facilitate the interpretation of large 
amounts of magnetic data. An example of a latitude profile 
PSovsnuOWl ia Frqure 222 (aj) This isa model profile’ Gi.e., it 
is the result of calculating the magnetic field due to a 
hypothetical current system), for an ionospheric current of 
S-eein latitudinal extent, andi20° in longitudinal extent, 
and with uniform current density. Figure 2.2(b) shows a 
Simiiar profile, but for a case where the current density is 
distributed in latitude according to the expression 
J(X) = exp(-3x?) where x is the latitudinal distance from 
the centre of the current. Both of these current systems are 
connected tc the magnetosphere by field-aligned currents 
with closure in the equatorial plane as shown in Figure 
2.2(C). Note that the north-south or X'-component peaks 
directly under the centre of the electrojet, and that the 
positions of the poleward and equatorward edges of the 
current system are provided by the Z-component extrema. The 
profiles shown in Figure 2.2(a) and 2.2(b) are calculated 
for ionospheric current flowing in a westward direction. For 
current flowing eastward, the signs of all the perturbation 
components are reversed. 

As well as such an east-west current system, current 
may flow in a north-south direction in the ionosphere, 
linkiag the two anti-parallel field-aligned Birkeland 
current sheets. An idealized picture of such a current 
system is shown in Figure 2.3. Kisabeth (1972) has computed 


Hatatude profiles for short ({i.e., limited in longitudinal 
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Figure 2.2 (a) 


Figure 2.2({b) 


A model latitude profile for. a currentyot 5°? 
latitudinal extent and 20° longitudinal 
extent, and of uniform density. The total 
GUurrent,1S 10¢ WA “and the) protijte: 15) Ssaken 
along a line 5° east of the central 


meridian. (After Kisabeth, 1972}. 


Situlernero figure e2z.2{aj, but for 4, ecurrent 
density distributed according to 
J(x) = exp(-3x2), where x is the distance 
from “the (Latitudinal center of thersysten. 
(After Kisabeth, 1972). 
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system 
involving an ionospheric current system 
flowing in the north-south direction. {After 
Kisaberas e157 2 ye 
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extent) north-south systems; such a profile is shown in 
Figure 2.4, along with the perturbations due to the 
contributing systems. As well, Figure 2.5 shows the magnetic 
perturbations for a very long (1809 in longitudinal extent) 
north-south current system. It will be observed that the 
general features are similar to those shown in Figure 2.4. 

limes, .AnSstrucrt i yeweat,: this point ito. consider.«<the 
Magnetic perturbations due to certain other hypothetical 
currents as well. 

For a current system which, unlike that shown in Figure 
2e2(c) does not flow along lines of equal latitude, the 
datitude profile produced along a meridian line differs from 
aeeacashOwn in Figure 2.2(a),.| For such a. “tilted" «current 
system, the perturbation in the X'-component is decreased 
while the perturbation in the Y'-component is increased in 
Such |.aaeway wethat .'2 + Ye ys ta constant.ssrrgure 2.6 shows 
the effect of constraining the ionospheric current to flow 
along the path described by @=@ rice, -9,) (1 + Cos y) 
where @ is the colatitude, Y the longitude (measured 
counterclockwise from midnight), and @O and 0, 
correspond to the colatitude athe (97 =180°, | and Ge 
respectively. For the profile shown, the poleward boundary 
has values of (o ee = re 20) while’ for ‘the 
equatorward boundary, ee eye Cyt) Ie lib t= 
apparent that AyY' follows AxX'. 

Later in this thesis, reference will be made to "net 


field-aligned currents", that is, currents flowing into the 
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three-dimensional 
as shown in Figure 2.2. 
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contributions or the component 


systems are shown as well, 


current 
This systen 


current 


, i 
+ eh ee ic Anadmtonnr — ene wha 
« ta pais “ . , at “ 


i asset | { i , 

1 ) 4 i} 
j om Lae a a Pres AGS i 

Pee ty ae S- wy) i es 


=% y 
e j 
| , 


Sous 


= 
— 


wt / / 7 i 4 
0 ii Li ait ie ; ee 
rer py ae hy ae om wet, Agia * ee 

ay J wn =i ve ua Pi) | 

1 Pe To aay we Se br = / 
ae : 

vy res 

ma ' i 


i ee Sh Ta 


He bag yl ri della la aaa ee - ee 


Be LR ROL ay: Gee Pee YS Mae 
Limite! oh kod . 


We: as 
cle 


=" t 


7 


v 
peste ainbaeigng 


4 hie rr + 

i ; ‘ m i : > ay; 4 

h | te Pig 2 heh ware er Somos aes 
, a ‘ 7) Aine | 


: , if © ee im Js. 
* a ae if : : - a va 
: A . a y seu ; ” ¥ ° 

ae yt S a ae 


48 


N-S CURRENT SYSTEM 
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ALTITUDE OF !ONOSPHERIC SHEET = 115 Km ~ 
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Frgure 2.5 Model “Latitude, “profile -of “a very long 


north-south current system. The system lies 
between latitudes 65.5°N and 69.5°N, and is 
1809 ain longitudinal extent. At latitude 
67.59, the height-integrated ionospheric 
curcvent. density, “2S 1.0. Amr *, and) the, total 
Current in theivsystem is 738 x: 10° A. The 
profiles are located at: | 

a) the central meridian (909); 

b) 45° east of the central meridian; 

c) 90° east of the central meridian. 
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Figure 2.6 


A nodel latitude profile. f£0F  CHlerer, 
constrained to flow along the path described 
by 9 = 6, +Z7@,-O)Ci-cos 7) , where 
@ is the colatitude, ¢y the longitude se aud 
@O, and @O, correspond to the colatitudes at 
P=16800* “ands 100) respective. For this 
profile, the poleward boundary has 
(6), Op) =2415%,209)) and the eguatorward 
boundary has (6,7 @,) "=" (20°, 25°) 4) A tee a ee 
10° A as flowing between Longitudes OF re 

he--4profile 1s” calculared at 
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ionosphere along a magnetic field line at one longitude 
which fiow up the field iine into the magnetosphere at some 
other longitude. The magnetic perturbations which are 
produced by such currents are characterized by a level shift 
in the east-west component of the magnetic field across the 
latitudinal extent of the net current flow. Figure 2.7 is 
the latitude profile for net downward field-aligned current, 
gietrabated, uniformly over 5° of Latitude. and"20° of 
longitude. 

When recording magnetic data, one cannot easily 
separate the field perturbations due to ionospheric currents 
from those induced in the conducting earth. It has been 
useful in the past to attempt to compensate for this effect 
by modelling the earth by a sphere which is perfectly 
conducting up to some depth below the surface of the earth, 
and is an insulator above that depth (Bostrom,1969; 
‘Kisabeth,1972). This is only an approximation to the real 
Situation, but has the advantage that the model is time 
independent, i.e., the nature of the induced field is not 
dependent upon the time variation of the external field. If 
an earth of finite conductivity is used, then one must be 
concerned with the skin depth which is frequency dependent, 
as well as the phase difference between the magnetic 
snauction  ficlidsrand lithe! dnduced electric field (which is 
also frequency dependent). For an infinitely conducting 
Material, the skin depth and phase difference are constants. 


Kisabeth (1972) has presented results of this approach for 


dibne Mawiet ce ae pert) al aim wii ne 
ivan his Page Oroa pHH rl bass wake) . ad | 


eis nbdy aheRy aca eaibg <i saan, ome - 
hee. betet, 6 gay ee uae roimgaat “ae apaaiaia 
SAY OSE? o Vint pasa ae Pe saan ; 


J 


er. Ais exmyee er x ailivacagth sai wt toe 


2 


ede an boned ye flea enue ban 10 es 


ate oh 


4 ; er Sit I . 
; F as a | ie A OHO Shan vA 7 om, Ao , Way 
. ‘ . +e om ae a/ . ; inet 
. ees j 
- P. 4 } ee 4 | 
e Lassa TOono So ey! ‘ ¥: ee i "5.8 tee, a si i ~ . : 
a aie) ; nodaae i 
e f . | 
wee te sa 
SHI Ied De OMG OMG e OT MD sieheies Sesh aka © 
s > & hee oT ah al ie % gl sas e ie ™ A " 


ph tie AE ie DR Page ay a, 
i. Bt 


“ vii pes oie NR de Pee i ie ced ‘aie ea a 4 


{ 


2 Pos tise as ine i erie! i ANAS: ie 


we 


i? ya fe pth fetid ssaeb 


“aeriypiotetty tata, "pall di me apne 


| F t ie 
4 ; f 
: i * ‘ ‘oe ae q 


4 aabas' wg anc SN ROS, ns wate oe Ey 


| ae Lift 
i we ag ’ ae Ze. es uy, 
mat eh, Leahou. Sr? ar ' ada haarp ca Bats ae! beat 
4: DN 


, ie the : oe us ist ae At Le yey "s | 
jon 2b Bist He el ad? aay ony ake dt) 
22 anda, At fits pee) Butt! ae: gettante, ink? san i 

J hed. hi a fe oe! aA ila id 


a at - Bifi. wed | Boge’ at theme ian ata 
a ae a ie: 0 { t od”, alt 
rncabnyok = ‘ake dhe de waueet ‘ikea att, 


nl Paes Ths ee ea on FF wae coe WAY Rows 


“el ee a Be. dogedan, estoy. ‘west 


fag Sore 
e 


ait “ie bie ‘aad. wg eden 
- ’ eo Pon “sil aS ae e at cues 


guys Seid: | RAS * 7 A OF act 


-y Mae 


Liye ie ay 'e oe herr es 7 
warraorticncis ane. pbaauer dia. fie pee ann 
a 


202 dowels, eke, a 


ga bhai 02 seni 
an fe hott te ao 


View, 7 - er. ih 1A j . - : * iN ae a2 
a | NY oh?) se ; zi oe 
ot ’ | + 


. wry rk eet "hy 


re 
Fa 
> 
— 
be 
wm 
Fa 
uJ 
ood 
ge 
col 
[ ae) 
aa | 
uJ 
mH 
Ww 


oO 
ie) 


DISTRIBUTED UNBALANCED FIELD ALIGNED 


CURRENT 

LONGITUDE: 0° TO 20° Zz 
LATITUDE: 70° TO 75° eS 
ALTITUDE OF LOWEST PENETRATION [15 Km eo) ; 
DEPTH OF SUPERCONDUCTOR 600 Km S 
TOTAL CURRENT 10° A e 
eee 7 
-. iL 
| Bee 


15S 60 65 70 75 80 
CENTERED DIPGLE LATITUDE 


1S. 
to 


85 


10° 30° 


jeeeog, 


(NT) 


wooo” 


Se Be DB Se DS eed 


FIELD SINTENST Tey 


55 60 6s. 70 75 S01 ees “55 600.) 65 70 75 80 gS 
CENTERED DIPGLE LATITUDE CENTEREOU DIPOLE EA rue 
Figure 2.7 Model latitude profile for a downward 


field-aligned current distributed uniformly 


over 5° of latitude and 20° of longitude. 


A 


total of 106 A is flowing, and the profiles 


are taken at: 

a) 10° east of the central meridian; 
b) 20° east of the central meridian; 
c) 30° east of the central meridian. 
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daitterent depths of an infinite conductor, and these are 
reproduced in Figure 2.8 . It will be noted that the x'- and 
Y*-components are enhanced by the induction and the Z- 
component is reduced. Although this approach involves a 
Rather risky approximation to the earth's’ conductivity 
structure, insufficient detail about the conductivity 
structure under the University of Alberta magnetometer line 
is known to warrant the use of techniques to separate the 
external and internal components of the perturbation field. 
However, it is believed that the techniques used by Kisabeth 
provide at least a first order correction to the data making 
any inferences drawn more reliable than they would have been 
were not some correction for induction taken into account. 
The latitude profiles of real data presented in this 
thesis are smooth curves drawn through the data points. 
These curves are meant to serve as an aid to the eye and do 
hot necessarily represent real data between the labelled 
points. However, the curves are drawn based on experience 
wath a very large number j/-of ‘stich profiles.) Since the 
discussion of the data will be directed more towards the 
Statistics of the samples rather than the individual events, 
it is believed that it is reasonable to discuss the profiles 


as drawn. 


Zecee POlar plots 
Although for analysis of the data it has been found 


that latitide profiles represent the preferred method of 
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FIGuUrEe 12078 


Model latitude profiles showing the effect 
on the computed components of the 
perturbation: magnetic “-fieid’ “ot placrrame 
Superconducting sphere at different depths.’ 
The panels on the left show the total field 
and taose on the right the field due to the 
induction “only. - These profiles dre tomeanu 
east-west current loop, spanning the 
latitude range’ \65°N (to (PS°NS A tora aes 
106) F is flowing over. 20° of. longitude... the 
profiles .are taken ~5° ‘east of the centras 
meridian. (After Kisabeth, 1972). 
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data presentation, frequently in auroral physics, current 
SyVeceissarecrdiseplayed Immpolarg plots. ane this#iithesi ssi) this 
method of displaying the data and model results will be used 
Oh occasion. In these plots, the horizontal magnetic field 
vector has been rotated through 90° so as to represent an 
equivalent horizontal current, and the equivalent current 
vectors have magnitudes proportional to the strength of the 
magnetic perturbation. Note that the equivalent currents do 
not necessarily represent real currents; some portion of the 
observed magnetic perturbations arises due to field-aligned 
currents which link the ionosphere to the outer 
Magnetosphere. However, this presentation permits easy 
comparison of the results of this study with some earlier 
work in which eguivalent ionospheric current systems have 


been used. 
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Appendix II describes briefly Backus-Gilbert linear 
inversion theory, and more detail may be obtained from the 
published papers onl Tthis topic! (Backustianhd® Gilbert, 
1967, 1970), and a recent review by Parker (1977). However, 
it should be emphasized that considerable care must be 
exercised when interpreting the results of an inversion and 
in particular when using the results of a "model that fits 
the data". As pointed out in the Appendix, one must be able 
to forward model the data before it can be successfully 


inverted. In the case of inverting magnetic latitude 
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profiles, one must have a current model which wili reproduce 
the observations of the magnetic field, subject oniy to 
finding the correct latitudinal current distribution. 

fo demonstate this, a latitude profile was generated 
for a hypothetical eastward flowing current with a 40 
latitudinal extent (colatitude 20° -to -24°): ‘and’ <a 209 
longitudinal extent. (Figure 2.9). The current was uniform 
in latitude and ciosed along field lines via an equatorial 
ring current. The profile thus obtained was inverted for the 
current density, subject to different modei parameters. 

In the ideal case, the geometry of the current systen 
is exactly known, as are all the currents whose fields 
eontribute “to the protive. If the data of Figure 2.9 are 
inverted, and the correct parameters (1.e€., the correct 
current boundaries) are supplied to the computer code, then 
the result is as shown in Figure 2.10. The solid line is the 
‘average <J(®8)>, and the broken line is the so-called 
flattest model (see Appendix II). In this example, the 
average agrees exactly with the known current distribution, 
as there are errors neither in the data nor in the supplied 
parameters. The flattest model, which is only one of an 
gnfanite set, ofjgmodelsye which Pnaya fit’ the data has) Heen 
constrained to have zero current intensity at the 
latitudinal extrema, and in order to fit the observations, 
this constraint introduces spatial oscillations into the 
nodel. 


The results of a second inversion of the data is shown 
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Figure 2.9 
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A model latitude profile demonstrating the 
inversion ot magnetic detdaa wo. Obtain 
height-integrated current density (current 
PRCCHSLtV me As, a. sunctton,. (OL), ,tatituac. JA 
total “cusrent of 108A owas distravuted 
uniformly over 4° (669N to 70°N) of latitude 
LOE da current Intensity /ori2.22 Am7?. 
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Figure 12810 
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Result of inverting the data shown in Figure 
2.9. The solid line is the curve of <J(6)>, 
the average current intensity. The dashed 
line 1schie WCUrrent) “intensity (for- the 
filattest model. In this example, the current 
limits were set to 66°N and 70°N. 
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Barer equine” 225114 Pesthwsy Case "the, “constraints ‘on “the 
latitudinal exten nave been relaxed, and this allows the 
Peattest model 10 £34 the data while still petmitting ‘zero 
current at the poleward and equatorward borders. In this 


case the model approximates the known current distribution. 
However, the current estimate no longer yields an exact 
current density, because by supplying the computer code with 
Seratcitudinal ‘extent that? 1S widemevthan ‘that used - to 
generate the data, the inversion puts current into a region 
where in reality there is a; current, thus. causing. an 
overall degradation cf the result. 

Finally, the data of Figure 2.9 were inverted assuming 
toat no field-aligned current or ring current contributed to 
be profrite. Although this is not physically realistic, it 
eyes to demonstrate the effect of uSing an incorrect 
Bode, 101, alternatively, ‘Of serrors ini the data.) For *these 
Pesuits (Figure 2.12), <J(8)> °is xreduced because yless 
SurrenG? .2s*4requzred =-to produce the magnetic field 
perturbations aot the field-aligned and fring current 
Contributions are ignored. Lhe flattest model in this 
example consists of large amplitude spatial oscillations, 
and bears little resemblence to the known current intensity 
distribution. 

These examples serve to demonstrate several features ODL 
the inverse PEODPemMs mtn the: first | place, each’ of the 
fhettest models Ofmcurrent edistribution, if —used@*inY the 


forward problem, will generate AX' and AZ profiles in good 
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Figure 2.11 
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Wi ore Lorri gure 2490." The ™constraante’ on 
the current limits has been relaxed to the 
latitude range 65.5°N to 70.5°N. 
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In this example, 


assumed that there was no field-aligned 
OF lind Current contributing to’ the magnetic 
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agreement with the hypothetical observations (Figure 2.13). 
However, each of these models differ in character 
demonstrating the non-uniqueness of the problem. Choice of 
an incorrect forward model can lead to grossly incorrect 
models that fit the data, although the unigue average, 
<J (9) > (1.e., unigue for a given set of parameters) in all 
cases gives a not unreasonable current distribution. Errors 
mieche “data do not. grossly affect the result of ‘the 
inversion of the latitude profiles. Errors in the data may 
arise from statistical considerations or systematic data 
collection errors. However, if an incorrect forward modei is 
chosen, or the parameters are incorrectly specified, then, 
relative to that model, the data are in error as well. If 
the data of Figure 2.9 had been real data, and had been 
inverted using incorrect model parameters, then the result 
of Figure 2.12 could have been obtained. The fit to the data 
is very good (the root mean square relative error is less 
current distribution, this model might be acceptable. During 
periods of strong magnetospheric activity, ionospheric 
conductivity and electric fields may become quite intense in 
spatially localized regions, and in such cases, one might 
hot expect a smooth current distribution. However, the data 
suite used in this thesis consists of hourly averaged data, 
and the averaging process will tend to smooth out intense 
latitudinally localized auroral variations. This assumption 


together with the knowledge of how spatially oscillatory 
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Figure 2.13 


Comparison of the model data from Figure 2.9 
with that calculated from the flattest model 
current distribution. 

a) Results from the flattest model o£ Figure 
Pack Wakes 

b) Results from the flattest model of Figure 
Zale 

c) Results from the flattest model of Figure 
PME 
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models may arise justifies the search for smooth current 
distributions, and this criterion was applied to the 
inversion results to test the correctness of the forward 


model. 


Although the latitude profiles used in this study 
consist of data points which are averages, and although 
profiles from a given hour have very similar qualitative 
Characteristics from dayweto rday,)srthere rors ,istill «san 
muavordablesvaniabiibicy whichsis; difficult tojsquantaty. pain 
order to obtain what might be called a set of generalized 
latitude profiles which would facilitate modelling of the 
source current systems, the method of superposed epoch 
analysis(SPEA) was used (Chree,1912,1913). That is, latitude 
profiles for the same hourly interval on each day were 
‘examined for a common feature, the latitudinal reference 
point. Values of each curve on each profile were then 
obtained at angular distances of +5, +10, +15 degrees from 
the reference point. In this study, the reference point was 
the peak value of the negative X'-component in the pre-noon 
sector and the peak value of the positive X'-component in 
the post-noon sector. in the noon sector, the X't-component 
was not well defined, and the reference point was chosen as 
the base of the Y'-component ramp that was always observed 
in this sector. Averages were then computed for each 


component at each relative latitude. The resulting profiles, 
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one for each hour of the day, are what will be called the 
SPEA latitude profiles. Note that the SPEA profiles involve 
a fair degree of subjectivity. That is, curves were drawn 
through the data points of the original averaged latitude 
profiles, and as pointed out earlier, these curves are drawn 
Subject to a certain degree of Subjectivity , although 
Prem ecurves, do" spass) through, (cachn, data, “point, and are 
constructed? with the Characteristic behaviour oftypically 
observed profiles studied over several years in mind. 

There is obviously a severe problem in presenting large 
data suites in which subsets of the data demonstrate 
qualitatively similar features. The actual number of data 
Poincs used in thes” thesis 2s of “the order of °5.4°x 107, and 
it would be unreasonable to present each of these. The 
. superposed epoch technigue allows one to show a large amount 
or —-data ji concisely,- although it must be borne in mind that 
the method is designed to enhance certain features of the 
fata.) wif taddi tion, the SPEA may mask other features of the 
data. For example, consider the latitude profile shown in 
Figure 2.6. When it is digitized as described above, in 5° 
steps away from the AX' peak as reference point, the AZ 
SceECtamedom Not tall = at: d7qrttazatlon points. in fact, the 
resulting profile has Z-component extrema 10° apart (Figure 
2e1k)- giving the suprescion..Or an electro jet that 1s) 10° *in 
Tatrtidina iy extent mkieun ssOther WOrdS, da * 5% "digitization 
increment is insufficient to resolve features less than or 


equal to 5° apart. However, the SPEA profiles will be used 
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SPEA PROFILE 


DEGREES 


The latitude profile that results from 
digitizing shiqures, 2.6 at S° intervals away 
from the peak in the X'-component. Note that 
the 442 extrema appear more widely separated 
than in Figure 2.6, although the character 
of AX' and AY' remain essentially the 
same. 
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primarily to provide estimates of the magnitude of the X'*- 
component peak value, and the size of the level shift in the 
¥'-component. Neither one nor the other of these are 


affected by the 5° digitizaticn step. 
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CHAPTER 3 RESULTS OF DATA ANALYSIS 


For this study, latitude profiles from hourly averaged 
values (centered on the haPirhour)tof (X's s¥t wand 27. iwere 
Uovrelated itor Hhevtperiod from Day.) (332, TORU attcoO, Daye 2a, 
iowa pounce esthemwapsolute:s fields not imcasured ;) (baseline 
values were required from which to measure the tnagnetic 
Peiaq perturbations. “To “this end, values of the three 
components at each ernneed one determined over the hourly 
interval 0500-0600 ywniversal time (UT) (i.-e., approximately 
2100~2200 magnetic local time (MELT) jon) days*<during which 
the observed magnetic field exhibited only small 
PescurbpaLLlonus Guring stherdayn(a iso-Gablied, "quiet" day). “For 
the data Suite in this study, baseline values were chosen on 
Day (tio, Sao ee Se hose) aiand Way) Joho) Ze a iKhp=7=). tana 
these baselines were used over an interval of +2 weeks, 

Single station magnetograms for the oe period of 
this study were examined for obvious drifts in the baseline 
values. When they were detected, a straight line was fitted 
through the data for the period during which drift occurred, 
and the data were then referenced to this line before 
removing the quiet day baselines. For example, from Day 35/7, 
hogiletorDay) toyetio72) atimonte Smiths the X"-component ¥, showed 
coe approximateiy Jdameer qiritt: ofve25 nil ithe i" -componentiva 
arneet yore 150r ntivand) ¢ herve =component::a 6 draft 1of) =—60: nT. 


Similarly, the data from Leduc exhibited drifts of 670 nT in 
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the X*-+component, 150 nT in the Y*-component, and -150 nT in 
the Z-component over the time period from Day 341, 1971 to 
Day 350, 1971. Ali the other stations had stable baselines 
during the observation period. 
Tivethise=cChaprelr,emaanGescraption. of atypical Jaritude 
profiles as a function of time will be presented, and a 
phenomenological modei of the currents as inferred from this 


data will be developed. 


Examination of the: |‘entire data suite, led» to’ 2a 
decomposition of each day into several sectors, eaca 


characterized by its own type of latitude profile. 


3.2.1 The 2290 - 0200 MLT Regine 

The first such sector or regime spans the period from 
moOut 2200 SET to about O200RMLT, facross the) region... ofy the 
Harang discontinuity. Representative profiles from each of 
the 4 hourly intervals are shown in Figures 3.1(a,b,c,d). It 
is evident from these profiles that there is a great 
Variability in the magnitudes of the average perturbations. 
In spite of this, in general, AxX* and AZ inthis regime 
are indicative of either westward flowing current, or both 
eastward and westward flowing current. It is probable that 
the magnetometer line senses end effects of electrojets 
which are not directly overhead, so that a Signature such as 


that shown in Figure 3.1(a), where AXX' undergoes a sign 
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profiles vfor the hours: 

(a) 2200-2300 MLT (0600-0700 UT) 
(b) 2300-2400 MLT (0700-0800 UT) 
(c) 0000-0100 MLT (0800-0900 UT) 
(ad) 0100-0200 MLT (0900-1000 UT) 
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change does not necessarily mean that both eastward and 
westward currents are flowing directly over the magnetometer 
line. Rather, for example, eastward flowing current may 
exist to the west of the line while westward flowing current 
flows directly overhead. However, regardless of the nature 
| Or the X*= and Z=component profiles in this time Sector,” the 
A Y' signature is distinctive in that it is negative across 
the entire auroral oval. Also, there is usually a step-like 
character to MAY' at the equatorward border of the oval, as 
is shown in each panel of Figure 3.1. ; 

mney SPEA (of all ithe dataptor this ‘regime. (Figure 3.2} 
demonstrates the above characteristics as well. AS will be 
discussed later, negative AY' signature is interpreted as 


the result of poleward flowing Hall current diverging 


poleward within the eastward and westward electrojets. 


3.2.2 The 0200 - 1000 MLT Regime 

The second distinctive regime spans the dawn meridian, 
lying between approximately 0200 MLT and 1000 MLT (about 
n000 . to /1800 UT).J This time period is domimated by the 
westward electrojet, although frequently there is an 
indication of a weak eastward flowing current equatorward of 
the westward electrojet particularly during recovery from 
periods of enhanced magnetospheric activity (Rostoker and 
Hron,, 1975). Figugesersesta)) Sthrougnv"s.a (h)reiare typical 
examples of latitude profiles from each hour of this regime. 


Sheminterval £rom 0200, to 03800 MLT . (Figure. 12.3.¢a)) is 
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Characterized by a predominately positive MA yY' profile and 
a weak or non-existent level shift across the electrojet 
region in this same component. The SPEA profile (Figure 
3.4(a)) demonstrates that, on average, there is a level 
Shift in the Y'-component, although arene 3.3(a) does not 
show this. The MAxX' profile is consistent with a westward 
electrojet combined with an equatorward flowing north-south 
current system. This latter is suggested by the positive 
AxX' in the poleward part of the profile. Note that the Z- 
component is asymmetric about the horizontal axis in that 
the positive extremum is greater than the negative extremnun. 
This feature is characteristic of this time sector, and is 
consistent with observations reported by 
Langel (1974a, 1974b) of the perturbation magnetic field 
observed at high latitudes by the polar orbiting satellites 
Ogo 2, Ogo 4, and Ogo 6. 

FIGULeG ss. 3(D) 1S "a aac for the interval 1100 to 
1200 UT (approximately 0300 to 0400 MLT). Qualitatively, it 
does not differ significantly from the previous latitude 
profile. The SPEA profile for this interval (Figure 3.4 (b)) 
is also similar to that for the previous hour although there 
is a more distinct level shift in the Y'-component. 

The next hourly interval (1200-1300 UT; 0400-0500 NLT; 
Figure 3.3(c)) shows a clear level shift in the Y'-component 
although it is), small. relative, to the maximum Ax 
perturbation. This level, shift, and its relative magnitude 


is preserved in the SPEA profile shown in Figure 3.4(c) 
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In fact, the remainder of the panels in Figure 3.3 are 
not unlike Figure 3.3(c) insofar as they all indicate an 
indisputable ievel shift in the Y'-component. As well, this 
shift increases in magnitude relative to the X*-component 
negative extremum as the local time of the profiles 
approaches noon. Also, the asymmetry in the Z-component that 
was described above persists as a general feature. The 
remainder of the SPEA latitude profiles in Figure 3.4 bear 
out this description as a general feature of this entire 


Magnetic local time sector. 


3.2.3 The 1000 - 1400 MLT Regime 

The third distinctive region which has been delineated 
spans local magnetic noon from about 1000 MLT to 1400 MLT 
fic00 UTeto 2200'UT)./ this region corresponds to, the “zone 
of confusion" identified by Harang (1946). An examination of 
‘latitude profiles from this sector bears out this 
description. Figure 3.5 shows a series of latitude profiles 
for the same day from 1000 MLT to 1400 MLT. The Z-component 
in Figure 3.5(a) is consistent with a weak eastward flowing 
current centered at about 67.59 N, and a stronger westward 
flowing current centered at about 73° N. However, the X'- 
component profile is not consistent with this. Similarly, 
the Z- and i commer ee ins Fiigures3.5(b) are not consistent 
with any simple eastward/westward electrojet interpretation. 
andeed) none of theipiprofiles. inp) Figure w3.5 can be 


interpreted in any simple way, and this difficulty exists on 
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ail days examined in this sector. The one feature that 
appears consistently in the latitude profiles from this 
region is the clear, relatively large positive-going level 
shift in the Y'-component. 

ihe SPEEA Of the profiles from 1000 to 1400 MLT (Figure 
Sia O(a, Dy Syrd)e) pa venyed pehicult to rexplain ianyterns Yof-any 
Simple east-west current flows. However, again the level 
Shatte -in Zi ome Lougit out clearly. ais isthe overran) 
Pegative bias i1n the xX*-~component. The level” shift “in the 
Y'-component actually appears more aS a ramp but this is 
Hore» an artefact or the SPEA than aj} feature of the real 
fate. Sliat is, ) in this time’ Sector, es ‘ a real level 
shift in AY", but the latitude range over which this 
Sccurs . Varies considerably from. ‘day, to day so.that «the 


Heaviside-like Yt-component becomes smeared out. 


3.2.4 The 1400 - 1900 MLT Regime 

The fourth regime which has been identified is the 
post-noon quadrant, lying between approximately 1400 to 1900 
MLT (2200-0300 UT). This regime presents latitude profiles 
On which ~@Xx" and *“A\Z can “be interpreted in terms of ~a 
simple eastward electrojet. For example, Figure 3.7(b) , the 
Potiuude, protaile, for 2300 to: 2400 UT (1500 to 1600 BLT) on 
Day 10, 1972, may be interpreted as an _ eastward current 
centered about 67.5° N, combined with a poleward flowing 
current to give the negative Jeter poleward of 70° N. 


Further similar examples are shown in the remainder of 
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Figure 3.6 Superposed epoch analysis of the hourly 
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Figure 3.7 Representative hourly average latitude 
profiles for the hours: 
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Figure 3.7. The Y'-component in this time regime is variable 
in that, part of the time, there is a positive-going level 
Ssiitt) (C.g., Figure 3.7 (b) and (c)), and other times there 
is no level shift (e.g Figure 3.7(d)). The SPEA of the data 
in this ‘sector indicates prior to approximately dusk a 
positive-going level shift in AY" (Figure 3.8). Near dusk 
(Figure 3.8(e)), the Y'-component tends to follow the x'!- 
component, indicating that the ionospheric current likely 
does not flow in a direction orthogonal to the magnetometer 
chain. Beyond this time in this regime, the SPEA gives 
peOfiles ‘in which the ¥'-component Shows a negative-going 
fever Shift although itvis' is not as clear a shift as 
observed in the pre-noon sector. However, some of the 
examples of the original data shown in Figure 3.7 show that 
indeed a negative-going level shift does occur in many 
instances although the majority of profiles up to 1900 MLT 


enow nO distinct jevel shift in AY". 


3.2.5 The 1900 - 2200 MLT “Regime 

The fifth and final region that has been identified 
from this study occurs in the evening sector ay about 1900 
Tom 200 MLT (0300 to 0600 UT) .) This regime ts .dominated by 
the positive NX" ‘Signature “of the \castward electrojet. 
However this tame) period ) ds “distinctive” in “that there 
exists a negative-going level shift in the Y'-component, and 
eumregion of negative: .A\ xX" poleward) of the positive Ax° 


regime. Although in the previous sector there also existed a 
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Figure 3.3 Superposed epoch analysis of the hourly 
averaged latitude profiles for the hours: 
(2) > 1900-75 00S NLP (2 200—25) Cau 7) 
(b) 1500-1600 MLT (2300-2400 UT) 
(Cc) 1600-1700 TET (OC0U-—0 10082) 
(d) 1700-1800 MLD (cO100-0206" U7) 
(e) 1800-1900 MLT (0200-0300 UT} 
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Region of negative “Ai"> in the poleward ‘part of the 
profile, the X*"-component ere ot peak in ths region, and 
the Z-component was consistent wen only an eastward 
current. The AY' level shift in the 1900-2200 MLT regime 
occurs over a narrow latitudinal width across and poleward 
of the eastward electrojet positive X'-bay signature. 
Figures 3.9(a,b,c) are typical of each hourly interval in 
this time sector. The corresponding SPEA profiles (figure 
3.10(4a,b,c )) again serve to emphasize the general nature of 


the observations in this regime. 


3.2.6 Summary of the Five Regimes 

To summarize so far, hourly averaged latitude profiles 
of the perturbation magnetic field may be divided into 5 
distinct regions, characterized primarily by the behaviour 
of the east-west or Y't-component of the field. Around 
‘midnight, including the region of the so-called Harang 
discontinuity, the Y'= component is negative over the width 
of the profile, and the X'-component shows the signature of 
a westward electrojet and, at times, of an eastward 
electrojet as well. From 0200 MLT to near noon, there is a 
deStanct level shift in AY', positive-going with 
increasing latitude. This level shift is confined primarily 
to the region of the westward electrojet and increases in 
Magnitude (with respect to the maximum Li) a Sa noc. 42s 
approached. The third region spans noon and is characterized 


by great variability in the X'- and Z-components. However 
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PEG Ure 6:7 Representative hourly average latitude 
profiles for the- hours: 
(a) 1900-2000 MLZ (0300-0400 UT) 
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Figure 3.10 


Superposed epoch 
averaged latitude 
(a) 1900-2000 MLT 
(bjs 2000-21000 kL 
(C) rez O=22: 007 Bile 


analysis of the hourly 
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Ay consistently shows a large positive-going step or 
fever silit.. Fron avout 100) to, 1900 MILT, the profiles are 
Characterized by the signature of an eastward flowing 
current and either a small positive-going Step tn Zu or 
DOomscepeate dit. from 1900 to 2200 MLYT, the profiles provide 
evidence of both eastward and westward currents ( the latter 
being the most poleward current) sive With these, a negative- 
going level shift in the Y'- component. 

the remainder jot -this «chapter wild. deal with the 
interpretation of these profiles in terms of a worldwide 
three dimensional current system, as well as presenting sone 
other information derived from these profiles and other 


SOUrTCEeES. 


3.3 Interpretation of the Data Suite 

The most distinctive feature of the data is the level 
Shift which occurs in the Y'-component. As described in 
Chapter 3, a net field-aligned current produces as a ground 
Magnetic signature, a step in MAY'. Field-aligned currents 
were tars proposed by Birkeland (1908, 1913) and 
subsequently were reintroduced by Fejer (1961), Kern(1962) 
and Bostrom (1964) in three-dimensional current models for 
Auroral) Treqions.) It eWwas not until recently (2muda -ct al, 
1967; Armstrong and Zmuda,1970) that such currents were 
identified from in situ measurements by magnetometers on 


board polar orbiting satellites. Several more recent studies 


(Yasuhara et al, 1975; Sugiura and Potemra,1976;Tijima and 
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Potemra,1976) show that there is a tendency for there to be 
het downward current flow in the morning sector of tke 
auroral oval, and net upward flow in the evening sector. If 
this is indeed the case, and the net field-aligned currents 
are uniformly diverging in the ionosphere, then no ground 
based signature would be observed (Vasyiiunas,1972). 
Fukushima (1969) has given a simple description cf this 
phenomenon by considering a line current flowing into the 
ionosphere (Figure 3.11). If the total inflow is I amps, 
then at a radial distance r freom the current, on the ground, 
the, magnetic frield’'" ys given by I/r,.<ande—-its'. directed 
clockwise (looking along the current). If this current 
diverges uniformly radially outward (in the horizontal 
plane), then the overhead ionospheric current density at the 
same observation point is I/277Tr, and the magnetic effect on 
the ground is given by i/7r, but in the counter-clockwise 
direction. The net effect is thus zero magnetic perturbation 
below the ionosphere. However, the ionosphere provides a 
highly conducting channel in the form of the auroral oval, 
so that one would expect to see distinctive magnetic 
Signatures at ground based observatories. These will have 
the basic characteristic of a level shift in the east-west 
component of the magnetic field across the region of net 
field-aligned current flow. fTfhis face together with 
Satellite observations of field-aligned currents has been 
used to aid in the interpretation of the latitude profiles. 
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Pa gure. 3.'11 


ZERO 


Schematic drawing to demonstrate how a 
vertical current which diverges uniformly 
into the horizontal plane produces no net 
Magnetic perturbation. See text for details. 
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Signature of a net downward field-aligned current, and a 
Begative-going, AVY‘ step to be the Signature of upward 
flowing current, then a picture of the diurnal variation of 
the disturbance component of the magnetic field associated 
with the field-aligned current flow may be obtained. Such a 
p1cture is shown in Figure 3.12. This is a histogram of the 
frequency of occurrence of inward and outward net field~ 
aligned current, as well as the number of cases when no such 
cases were observed in the ground data. The cross-hatched 
region delineates that region of the auroral oval known as 
the Harang discontinuity. 

There is a remarkable similarity between the behaviour 
of the net field-aligned current as inferred from the 
latitude profiles, and the diurnal variation of the average 
electric field observed at auroral latitudes (Mozer and 
Lucht, 1974; Iversen and Madsen, 1977). Figure. 3.13(a) is a 
plot of the east-west and north+rsouth components of the 
average auroral zone electric field as a function of local 
time obtained by Mozer and Lucht (1974). Most of the data 
used in compiling this plot was collected at Thompson, 
Manitoba, and Churchill, Manitoba, where local time is 
approximately the same as magnetic local time. Figure 
Bets(D) is Similar data an a polar plot, obtained: by Iversen 
and Madsen (1977). 

The first region discussed above (2200 toe 0200 Ski) 
features an electric field which is westward and poleward 


(towards dusk), and westward and equatorward (towards dawn). 
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Figure 3412 


Histogram showing the distribution of net 
field-aligned~ current) as "4 9fLuncize, of 
Universal Time, as inferred from the 
ground-based hnagnetic measurements. 
Approximate magnetic local time 1S shown 
across the top of the figure. 
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Figur e--3 7. 13.(a) 


The diurnal variation in the ionospheric 
eléctric fi6¢1d in the region of the aurora 
oval. The time scale is local time, which, 
for, these data, is approximately equal £0 
magnetic local time (see text). (Arter foeer 
and bOcnt, Vo Tae. 
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Figure 3.13(b) 


LO2 


Average ionospheric electric field 
ini «the. region. of thes auroral zone, 
in a geomagnetic latitude, magnetic 


time _ frame. (After Iversen and 
Madsen, 1977). 
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This field can drive a poleward Hall current and a westward 
Pedersen current. In this region, the latitude profiles show 
a Y‘'-component which is negative across the profile, and 
this signature can be produced by poleward flowing current. 
The second region, or dawn regime has a dominantly 
€quatorward component of electric field and the east-west 
component, although small, is primarily eastward. This 
equatorward electric field can drive a westward flowing Hall 
current, consistent with the magnetic Signatures observed in 
this» sector. The noon sector (1000 to 1400 MLT) is 
Characterized by a weak or zero east-west electric field, 
BPuoagasutransitionaiw inorth=south" electricoefielid.* "In: “this 
region, the ground based magnetic observations are confused 
With respect to electrojet flow, although the Y‘'-component 
Signature is that of a strong inward net field-aligned 
current flow. In the post-noon quadrant, the data of Mozer 
‘and Lucht indicate an essentially zero east-west component 
toethevelectracitield;. butiva, Strong poare vane component. 
This, too, is consistent with the latitude profiles which 
show in this sector an eastward electrojet that would be a 
Hall current driven by a poleward electric field. The data 
of Iversen and Madsen (1977) show a small eastward component 
in the electric field. However, the near absence of an east- 
west component in the electric field is consistent with 
essentially no north-south component of Hall current. Thus 
one would not expect the electrojet Hall current to diverge 


Significantly°*along field ‘lines’ in this region. Finally, 
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toom 1900 to 2Z000MLT the. electric, field \ has “a ‘westward 
component 1ncreasing toward midnight and a poleward 
component increasing toward dusk. These fields are also 
consistent with the observed eastward electrojet flow that 
is observed in this sector. Table 3.1 is a summary of these 
comparisons. 

As pointed out above, the level shift in the Y'- 


component ( AY! is believed to be related to net 


step! 
field-aligned current flow into and out of the ionosphere. 
The positive-going AY step in the dawn and noon sectors 
is then interpreted as the Signature of net inward field- 
aligned current flow. To maintain current continuity, it was 
hypothesized that this net inward current feeds the westward 
ionospheric electrojet. Similarly, the net inward flow in 
the immediate post-noon sector could feed the eastward 
electrojet. In the dusk sector, where AY! step as 
hnegative-going and apparently related in position to the 
negative X' latitudinal regime, the inferred upward field- 
aligned current has been hypothesized to be the result of 
the electrojets bleeding upward along the field lines into 
the magnetosphere. 

To test this hypothesis, it was assumed that the peak 


values of the X'-component « AX! were related to the 


peak ) 


magnitude of the electrojet current strength, and further 
that AY ctep was an indicator of the magnitude of the 
unbalanced field-aligned current. Subject to these 


assumptions, one would then predict that the norm of the 
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ratio. “of ZN to TD EL 


step peak from 0200 MLT to noon 


would increase and from noon to dusk, the same ratio would 
decrease. in the. dusk uadrant AY! xe 
would increase for negative \X', and decrease for positive 


A Xt. Figure 3.14 shows this ratio obtained 
erLOuUgh Othe SPEEA, plotted as. daLunction of magnetic Local 
tae. “Lt 01S) apparent: gehbat qualitatively, the above 
expectations Ware realized. Near moon, 1t was difficult to 
establish a value for A X', so no values are plotted for 
this regime. 

This analysis is important in that the results form a 
basis for a world-wide current system model which will he 
discussed in detail in Chapter 4. That is, it is believed 
Giat theydistinctive positive, level shift in™ AV Y' of the 
hourly averaged latitude profiles across local noon and in 
the pre-noon sector is a signature, at ieast in part, of 
unbalanced field-aligned current flow into the ionosphere, 
which feeds the ionospheric Hall current electrojets. fhe 
negative level shift in AY' in the pre-midnight quadrant 
reflects field-aligned current flow out of the ionosphere. 

An alternative way to examine these data is to consider 
AY Step and ee separately as functions of 
magnetuc local time. Figure 3.15 .1S the plot of the average, 
rot a . OGLven hour, o© AY. and AxX' k versus MLT. 


rep pea 


Note that for the dawn regime, Ax! nk. increases from 0200 
pea 


-1000 MLT (i.e.,it becomes less negative) ,by about a factor 


Gp three, while Avy? ‘ increases by a tactor of 2 to 3. 
step 
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Figure 3.14 


The ratio of AyY' tors ky! plotted as 
: g st : : 

a function of local magneti@@ “tine. LO. a 

first Tapprox imation; /Vy! Ts eenudicative 


of unbalanced field-aligiS@ current flow, 
and AX'peak is a measure of the ionospheric 
current flow. 
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Pitot! Wiof AY! and Ax? asta ee 
of magnetic 7581 time. peak. 3S. 

meanS a positive-going pete cee “Step 
south to north; negative means a 
negative-going level Se “gtep SOULhe to 
north. Positive AX peak is the nny ee of 
eastward flowing current, negative 

is the Signature of westward =r iBeek. 
current. 
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Simtlariy, (ih the “region of the ‘eastward electrojet, 


NX? - increase <4 : ‘ 

peak S;-until about dusk while AY Seep 
decreases. Beyond this time, Av step becomes increasingly 
negative in coincidence with a decreasing Bo ie Pen west 


these observations are consistent with the interpretation 
given in the previous paragraph. 

In the diurnaliveriation .of eet for the interval 
w200)tO 1000 MLY, there is a ‘distinct inflection ‘on the 
curve near local magnetic dawn. This inflection is thought 
to be due to a decrease in the total westward flowing 
current near dawn. Some of the current flows up field lines 
hear dawn due to a conductivity gradient between the sunlit 
and dark hemispheres, This observation will be discussed in 
some detail in Chapter 4. Finally it is noted that OnE" nak 
fae a negative extremum during the hour 0200-Q300MLT. This 
is in agreement with Alien and Kroehl (1975) who found that 
AL (indicative of the strength of the westward electrojet) 
tended to reach peak values approximately three hours after 
local magnetic midnight. 

The time sectors in which upward and downward net 
field-aligned current are found from ground based data are 
in good agreement with those described by Sugiura and 
Potemra (1976) (Faqure. 3.116) “except in the 1200-1400 (MET 
sector. These authors find net upward field-aligned current 
Te this sector, whereas in the present study, the reverse 


Current flow was observed, and this latter is consistent 


with a growing eastward electrojet. Most of the published 
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Figure..3.16 


Percentage occurrence of level shifts in the 
east-west component of the magnetic field 
observed by “the - Triad: satellite, i sacuea 
functivon) of MLT. The dotted regions 
represent net upward field-aligned current, 
the clear regions net downward field-aligned 
current. (After Sugiura and Potemra, 1976). 
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data from the Triad satellite have been from the summer 
months, whereas all the ground- based data in the present 
study come from winter months. An explanation of the 
discrepancy between the Triad satellite and ground 
magnetometer data will be offered in Chapter 4- 

Because of the role that the interplanetary magnetic 
field (IMF) plays in magnetospheric and auroral physics (see 
Chapter 2), it was useful to examine the correlation between 
the IMF and the presence of the AY' level shift observed 
in the latitude profiles. To this end, data were selected 
from two regions showing pronounced level shifts in the Y'- 
component. The hours 1600-1700 UT and 1700-1800 UT 
(approximately 0800-0900 MLT and 0900-1000 MLT) were chosen 
aS representative of positive-going level SivEts hnetihe ls 1.) = 
component, and the hours 0400-0500 UT. and 0500-0600 UT 
(2000-2100 NLT and 2100-2200) MLT) were. chosen ‘to represent 
the negative-going exis level Shirts. WFagure. «seas 
(a,b,c,and d) are plots of AY' as We BS, (1-5, 9the 
component of the IMF normal to the ecliptic plane). Although 
in some cases there were no data for the IMF for thoses 
Cases ..where Gata did: exist, there 12s. a clear Jinear 
relationship between AM step and B, for the pre-noon 
sector. That is when B, >0- (iee., the IMF points northwards, 
Guets parallel to the earth's ficid in the ecliptic, plane) 


there is a weak AY. ; however, if B, <0, then AYt 


cep 


increases linearly with increasing negative B,. If as 


cop 


Suggested above, the level shift in AAY' is related to a net 
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Figure 3.1/7 PI Ot S010 ties Ze Yas as a ‘function of “the 
component of the IMF normal to the ecliptic 
plane (5 1,0 for the hours 
(a) 0400-0500 UT 
(b) 0500-0600 UT 
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Figure /36 i 


PLOESROLG 12! as a function of the 
component of SheP rue normal to the ecliptic 
plane (B ), for the hours 

(Cc) TO00= C200 TUE 

(d) 1700-1800 UT 
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component of field-aligned current which in turn is related 
to the strength of the westward current, such a relationship 


between AY't and B 


step z would be expected. In fact, 


Hirshberg and Colburn (1969) have shown a good correlation 


between Kp, °a measure of auroral activity, and B when 


& A g 
Be is directed southward. 
ingsthe Sevening= sector;,<"the relationship between 


AY 


"step and B, 1Seenotivaseiclearpoalithougnhaabpattern 
Similar to that observed in the pre-noon sector is observed 
here. It would be of considerable help in determining 
unequivocably the relationship between AY step and B, if 


there were more cases of iarge AY step in this local time 
interval. However, for the period of this study, there was 
IMF dataceforekabout? only*h50Zavofrthe, days injeach hourly 
interval, and of this 50%, only about half of the cases 
provided large level shifts in AY'. The result of this 
‘unfortunate circumstance is a paucity of data involving 
large level shifts in the Y'-component. However, the results 
are consistent with the observations of others (Hirshberg 


and Colburn, 1969) in view of the interpretation of the Y'- 


component level shift presented in this thesis. 


3.4 Inversion of Latitude Profile Data near Dusk 

It) was ~pointed out earlier in this chapter (that there 
was evidence for penetration of the westward electrojet 
poleward of the eastward electrojet as far westward as 1800 


MLT. Rostoker and Kisabeth (1973) have shown that eastward 
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and westward electrojets exist simultaneously in this time 
sector during polar magnetic substorms. More recently, 
Kamide and Akasofu (1976) have published a similar result. 
As the existence of a westward flowing current in this time 
sector during quiet to moderately disturbed times is 
important to the understanding of magnetospheric processes, 
a detailed study of the current flow in this regime was 
undertaken using the linear inversion theory of Backus and 
Gilbert (1970). In particular, the formulation developed by 
Oldenburg (1976) was used to determine the latitudinal 
distribution of current in this time sector. Appendix II 
gives a brief outline of the theory involved in this 
technique. 

Besides carrying out the inversion, which leads to the 
only unique soluticn to the problem, a model of current 
intensity that “fits the data" was generated. That is, out 
for *the antinite ’ ‘set ~"or “particular? models, ‘one which 
satisfied the relation [3 6148 = minimum, and which was 
constrained to fit the observations that fit the data as 
closely as possible, was developed as an aid in interpreting 
the results of the inversion. There is not necessarily any 
physical reason for choosing such a model. However, such a 
choice is analytically and numerically easy to handle, and 
does satisfy the intuitive belief that the current density 
of auroral current systems, On average, will vary 
latitudinally in a smooth manner when averaged over the time 


Span of one hour. 
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To demonstrate the results of the inversion, one 
example from each hour from the sector 2300 to 0600 UT (1500 
to 2200 MLT) is presented. For each latitude profile, there 
Bre three accompanying figures. fhe first is a plot of 
Standard deviaticn vs polatitudee contoured for constant 
values of averaging function width. These plots are in 
effect a graph of the trade-off between error in the 
calculation cf the current density, and the resolution for 
the same calculation. This figure is required to interpret 
the second pilot which shows average current density as a 
function of colatitude contoured for a constant standard 
deviation. Finally, the "“flattest" model of height- 
integrated current density that fits the data is plotted to 
assist in the evaluation of the current density estimates. 
As discussed in Appendix Ii, the forward model used in this 
study employs both east-west current flows and north-south 
‘Current flows. Furthermore, in this regime, the electric 
field is approximately northward, so that the simplified 
relationship between the Hall current and the Pedersen 
current, as developed in the Appendix, has been used. In all 
cases, the ratio of height-integrated Hall conductivity to 
height-integrated Pedersen conductivity was set to two in 
accordance with the results of Brekke et al,(1974). The 
results of the inversion are shown in Figures 3.18 through 
B20. it is evident thatiprior to'1800:° NLT ,oF only castward 
current is detected, as shown by the curves of <J(9)>. It is 


interesting to compare the profiles in Figure 3.20 to those 
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Figure 3.18 Results of inverting the hourly averaged 
latitude profile LOL 1500-1600 MLT 


(2300-280 0 UT) / Day's, 1972s 

(a) The ati tude® protalevaswi ae function (of 
colatitude. 

(b) The flattest height-integrated current 
density model. 

(c) The standard deviation of the current 
estimates as a function of colatitude 

(d) Estimates of height-integrated current 
density as a tunction of colatitude, 
contoured for constant standard deviation. 
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Results of inverting the 
latitude profile for 
(0000-0100 00), Day 17,, 1972. 
(aj Thesisaticudens profile*asea Lunction of 
colatitude. 

(b) The flattest 
density model. 
(¢) The «standard ‘deviation of the current 
estimates as a function of colatitude 

(d} Estimates of height-integrated current 
density as a function. Of Scoletitude, 
corptoured for constant standard deviation. 
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Results of inverting the hourly averaged 
latitude profile for 1700-1800 MLT 
(010.0=0200) UT). Day 15, 1972. 

(a) The latitude profile as a function of 
colatitude. 

(b) The flattest height-integrated current 
density model. 

(c)) The: "standard deviation of the current 
estimates as a function of colatitude 

(d) Estimates of height-integrated current 
density as a function of colatitude, 
contoured for constant standard deviation. 
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(bPETH Sar lateest etheight~integrated” “current 
density model. 

(epi The “standard? deviation of the current 
estimates "as#a"tunction of colatitude 

(d) Estimates of height-integrated current 
density as a function of  .colatitude, 


contoured for constant standard deviation. 
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Pagure, 3.22 Results of inverting the hourly averaged 
latitude protile for 1800-1900 MLT 


(0200-0300 UT), Day 333, 1971. 

ta) The, latitude profile. as ,a .function, of 
colatitude. 

(b) The flattest height-integrated current 
density model. 

(cy) Dheastaucand deviations ,Ol ethe. current 
estimates as a £unction of colatitude 

(qd) -Estimates..of height-integrated current 
density as a function of colatitude, 
contoured for constant standard deviation. 
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Figure 3.23 Results of inverting the hourly averaged 

latitude profile fOr 1900-2000 MLT 


(0300-04004 UT), Dave 18,6 19725 

(a) The latitude profile as a function of 
colatitude. 

(b) The flattest height-integrated current 
density model. 

(ein Whe ostondards deviation,of the Current 
estimates as a function of colatitude 

(d) Estimates of . height-integrated current 
density as a function of “colatatude, 
contoured for constant standard deviation. 
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(a) (The Llaritude prorile as a  functvon “of 
colatitude. 

(b) The flattest height-integrated current 
density model. 

(c) The standard deviation of the current 
estimates as a function of colatitude 

(d) Estimates of height-integrated current 
density as a function of colati tude, 
contoured for constant standard deviation. 
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Results of inverting the hourly averaged 
latitude profile for 2100-2200 MLT 
(0500-0600 UT), Day 335, 1971. | 
(a) The latitude profiie as a function of 


colatitude. 
(b) The’ flattest height-integrated current 


density model. 

(cyjee the @Standard deviation of the current 
estimates as a function of colatitude 

(d) Estimates of height-integrated current 
density as a furnction-" of) *colatitude, 
contoured for constant standard deviation. 
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Results of inverting the hourly averaged 
latitude profile TOL 2200-2300 MLT 
(0600-0700 UT), Day 338, 1971. 

(a) ‘The: Latitude profliesas a «function ‘of 
colatitude. 

(b) The flattest height-integrated current 
density model. 

(c) The standard deviation “of” the ~current 
estimates as a function of colatitude 

(d) Estimates of height-integrated current 
density as a function of colatitude, 
contoured for constant standard deviation. 
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Piel Ouve nS 21 /LOmetno n0Ury 1700-1800 (NLT) in Figure: 3.24, 
there is a possibility of a westward current poleward of the 
iain eastward electrojet, -or:-this simply may be due toa 
westward current to the east of the magnetometer whose 
strength is’ large enough, so that the resulting magnetic 
perturbations can be seen. at the magnetometer line. For 
times milatem i than 1900.5 MULT | fithe “contours of <Ji(6)i> show 
unambiguously that there is a westward electrojet flowing 
poleward of the eastward electrojet. The validity of the 
forward model used is somewhat open to question as evidenced 
by the oscillatory nature of the flattest model. However, 
Such behaviour of the particular model does not detract in 
any great measure from the results of the actual inversion, 
as shown in Appendix II. Table 3.2 is a Summary of the 
results of the inversions of latitude plrofimes, <n) this 
sector. All cases included in this table are for times when 
no substorm activity was present in the magnetometer chain 
sector. The positive X'-component extremum varied from 10 nT 
to 100 nT, while the negative X'-component extremum varied 
EroOh1.)-nT-to-105-1'T. 

Since the examples analysed in this study were for 
quiet or only moderately disturbed periods, as evidenced by 
the range of the X'-component peak values, the westward 
electrojet that has been detected in the dusk sector is not 
the substorm westward electrojet. Rather, these results show 
that the convection westward electrojet may penetrate as far 


westward as 1800 MLT. Previous workers (Kisabeth and 


okay) a Oe TR Ni UGS maria Siete 
HA | as a i) Nie) ny ; ui iH si ig a ni — < 


whee i whl ih ab ie ol 


j 
Vi KM ey 
a He i 


sudan f slam asi aae ig Lie didi ‘ 
wane on + he ean a wa Mek ve 
Ree NO ee ey i Sal ipa io ane we we | 

. ~ ae “eye si 


4 ay Ne we) } int Ne iy : ay, i x 


i 


att i ie | abel wide i meet) ay hie. 4. tne 
| eae oe. pra e A get 


aks K 
ay 


: i nt hed a Fy hike s a rit es ie ih 
arn ( i ee bein’ O® f 
fk ae Gu | are re) in: He. “ahh a 
ra nN Rabie eh ett * \ ona nN i hp 7 PP een ay 4 Tepeeriges mn \) Me y 
ihe MO LE his: ar he srishes 
te f t i aoe * ae ‘ia aa aay’ ia) ie % a t- 5! ie, AD i ene ee bai 
say lf ir , ( } Ms ahs nt ak oe , i cents ong cee oti He Ren : viet 
, uy he | a atin i ; ry 
yh AN he | ye Sie ere earn Mf ‘i 
sit ys eae meh ies, 2 fe ae Oe ei wba one , | 
Bi S yh iN Dh oe 
ep a a a Ae dpi yf hh a a wa” i ae ‘ptoaot 4 | ; 
hie ak) ee co ey ae Aaa 5, woman ent: ag +f 
iy oor, isa ae ‘epithe Wnt ce bp betan re ‘poi: - 
en; cag eon tey MEU eh ae a + sinha new “ge bw rend tnt 
gr BP Aol Mate neha: ih wi sud, patecadg: ole 


vr Cit Papin TRY vn salts, ‘wh pbk oe wae: phate: - a 
3 , ait bt. Ho! Ot wo 


~ i 

Uj 
wie : +] 
hy We 7 an 


tak ok ee He TT a, - ubetbaee winged on 


TEAL de) a oy em rh . ial ha i. 
baddesinnws ohn sgael se, may. pias ee ue Py 
ae py) ae. ae iea 


it) a FI eve aa ink Lap tenity a salt 
_ it ni 4 Tomer dienes am oes 


bane Getanon, ann eat sgecanes: Aesniancndl 


an ‘ 


a a il aos ti aera dota 
ag ie ae ane ba 
ee iy A : eT 


| 


ies RGA LSD SB 8 


130 


ee ee ome a ee ne ee a) es ae ee eee eee ees ee oe 


ED ee rte ee ee ee ee ee ee ce SP ee re ee 


aT ee Oa A. CM es La Gee, ok : Tey eT Lee el 
{ Interval i Number of { CurrenteDirection, | 
{ (MLT) { Cases { (Eastward: E { 

{ { Westward: W) { 
-———— ne ber mae | 
i { { { 
i 1500 - 1600 {| 4 { E { 
{ { { { 
{ T6008 = 1700: > 4 4 { E { 
{ { { i 
in4700.=1 11800) j 3 | E ' 
{ { 1 I Et i 
i { i { 
{ 1800 - 1900 |] 3 | E+wW i 
I i { { 
919.00.<= 20009) | 4 { E+W { 
j { { i 
fy 2000 = 2100 4 3 { E+wW { 
i { { 4 
iz 1002= 2200" * | 3 { E+W { 
{ { { | 
te 200e- 23001 5 | E+ 
{ { { { 
isa De TF head AUER Rese eA Cee ial Ne ey 


an Ae if FY Pos fa 
ae 1 ee 
I am) aa 1 lh 
i ‘its ; nan 
A 2 ¥ : 
, 4 Mh 
VI Lise i a 
hy; LA es hy! ‘ 
2 bs Bet tie a ie ee = a 
pate itil ek mylene 
hi j 
es gd. 
an 
Ps 3 iy | oy eh 
. ke Lo i A Ae 
+ et 5 i oe ” ict “ 
y f oa j if 7 1 ha 
; : i jaar i 
: i ; Pe reeves i nee 
} : i! iM t i 
j by , 
ayy: ; 
Vela as . “ 4, ot ah oy poe e oe 
AL eae ce Tash Mime I) OT ices 2a aI eet bee Cie Hag "I oe 
f : 1 h 


wer onthe hoe hah : aa Jaen ays a eat sides vr ae : 

eat ee hd ier D a) ae ‘aes 
u. «brawP eae) H 59 i 
(W roe ngeer » | a i | 


are 


, i he } 
ena meets diy beni en ee rg ote elec lnethig salable 


@ 


i ad a oe ee ee oe ee ee oe 
— Se oe oe x 


g yy 1 
¥ wt Wy aD i ( 
nl ; j | : 
aed I, *) Na 
; Ae eae oe, 


= 
J 
rig 


: 


[ ye rerre vou hipaa dl tear oy eaehipee amt 


f 
\ 
¥ 
ei 
‘ : , 
ij : j 
| \ 
b 
il i 
a haz 
\* 
} j : rl 
4 
J bt 


| . | any i eT ey I eee: 
OT OLR ahi cee mn ane a a 
‘ a ' : 0 ( ? ; : 


3 a 


Rostoker,1973; Kamide and Akasofu,1976) have found substornm 
westward electrojets in this sector. The results of the 
inversion, together with these other results, make it 
imperative to include the westward current poleward of the 
eastward electrojet in the pre-midnight sector up to local 


dusk in any world-wide current nodel. 
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Chapter 4 A WORLDWIDE CURRENT MODEL 


4.1 Introduction 

in-\ this;) chapter, the,,,data,;;and its , interpretation 
outlined in the previous chapter, are synthesized to provide 
the input parameters for a world-wide three dimensional 
current model. The model is tested against the data inasmuch 
as the gross features of the data are reproduced. That is, 
no attempt is made en say, least squares fit a nodel 
latitude profile to an observed profile. Rather, the 
comparison that has been made is one of comparing the 
Magnitudes of the AAX' extrema and the level-shift in AyY'. 
This is justified in ao the model is designed to reproduce 
only the gross features of current flow in the ionosphere. 
However, it is believed that suitable variation of the model 
parameters will lead to current models that will reproduce 


specific latitude profiles. Several examples of such current 


models are described to demonstrate this. 


4.2.1 The Auroral Oval 

In order to develop a model of worldwide current flow, 
one must determine the direction of current flow and _ the 
geometry of the current systems. Chapter 3 has outlined the 


longitudinal distribution of eastward and westward current 
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flow. However, at is also necessary to define the 
ratitudinal extents of the .curnent.. systems. From the 
superposed epoch “analysis, and. using the ZZ extrema as 
indicators of the latitudinal boundaries of current 
(Kisabeth 1972), it is found that in general the east-west 
current systems are confined to a region about 5° wide in 
latitudinal extent. Furthermore, there is evidence that the 
eastward and westward ionospheric currents do not flow in a 
direction normal to the magnetometer line. That is, the Y't- 
Component. pin. the) sslatitude.” profiles: cannot) be attributed 
SOC LY. 8 tO pt ieid—alagned |.current, ~ibut 1s,in part .due) to) a 
fists O14 ' the-,smain electrojet Withigupespect) metro” elimes) for 
constant latitude (see section 2.2 and Figure 2.6). Indeed, 
the electrojets are observed to flow along the auroral oval 
(Feldstein, 1963, Akasofu et al, 1965). Kisabeth (1972) 
modelied substorm current systems flowing along the auroral 
oval by representing the current path by a parabola given by 

ab Avail 
(Bicosty + a’sin'g)* 


Where O is the colatitude of the boundary, (G7 is the 


longitude (measured counter-clockwise from midnight), @Q is 
the midnight colatitude (Y =0) and 6 is the dawn or dusk 
colatitude ((P=90° or W=270° respectively). However, he was 
interested in currents of relatively short Aone ea ne 
extent (up to 20°), whereas the work presented here deals 
with currents of global longitudinal extent. Kamide and 


Fukushima (1970) have examined current flow along a path 
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given by 

| 0= 0 +Z(0,-8,)(1+cos¢) 4.2 
where O ana P are jas) given above, and “O. ais) the 
colatitude of the path at noon (P=180°%) and O, is the 
colatitude at midnight (QY@=0°). Since small scale current 
Structures are of little concern in this thesis, equation 
4.2 1s ideally suited for describing the auroral oval since 
a single specification of 8 and 8, suffices to describe an 
entire oval boundary. 

In the midnight sector, it is believed that there is a 
component of poleward flow and that both eastward and 
westward electrojets coexist. To accommodate such currents, 
a separate specification of boundaries is required, but 
equation 4.2 may also be used to describe the path. Figure 
#21) is apolar plot of the auroral oval boundaries” used in 
this thesis and tabie 4.1 gives the parameters that describe 
‘these boundaries (Since the polar plot is in terms of 
atitude,, | the table, Jvsts Blatitudes 27 X%) instead -of 


colatitude (@)). 


4.2.2. Location of Currents 

The nature of the ionospheric current flow is based on 
Ene. interpretation: ,of the» Latitude profiles *givenmin the 
previous chapter. Westward current flows from noon around to 
midnight where it flows poleward and westward across the 
Exvdnaghts sector £0 92200 ML. This flow then continues 


westward to 1800 MLT. Ractward current flows) from Toon 
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Figure 4.1 


Polar plot of the auroral oval) useq aman 
modelling purposes. The boundaries follow 
the locus| 0 = 0 +44 (0,-8,)(L+cosx) where @ 
is colatitude and @Q iongitude, measured 
counterclockwise from midnight. OQ and & 
are the colatitudes at noon and midnight 
respectively. The values of @ and 8, are 
given in Table 4.1. 
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af0und (through dusk» to 2200 MELT and then northward and 
eastward to the equatorward boundary of the midnight sector 
westward electrojet, terminating at 2400 MLT. Both the 
eastward and westward electrojets are one segment of a three 
dimensional current loop similar to that used by Kisabeth 
(1972). However, the field-aligned currents connecting the 
ionospheric electrojets to the magnetosphere are not 
confined to sheets at the ends of the ionospheric currents 
but are in general distributed over the length of these 
currents. Thus in this model, the electrojets are assumed to 
Grow ‘and decay over a finite longitudinal extent, as 
described in Chapter 3. 

From the nature of the electric field configuration in 
ane auroral zone (see figures 3,12 and 3.13), it 1S: apparent 
that the eastward and westward electrojets are essentially 
Halil currents, so that for a complete current a aie, one 
Must also account for the currents in the direction of the 
electric field, the so-called Pedersen currents. Indeed, 
recent observations by the Triad satellite (Iijima and 
rOocenra, 19/6) andicate that the auroral oval is bounded by 
regions of field-aligned current (Figure 4.2). These 
currents close in the ionosphere through north-south 
currents. In the modei which has been developed for this 
thesis, all east-west ionospheric currents are bounded by 
infinitesimally thin sheets of field-aligned current which 
are in turn connected by north-south ionospheric current. In 


regions of westward current flow, there is downward field- 
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Figure 4.2 


The location of field-aligned current sheets 
as deduced from magnetometer data from the 
polar-orbiting Triad satellite «From Iijima 
and Potemra, 1978 ) The panel on the left is 
for low level activity (AL<10CQ} , an@ thazeon 
the right for higher a@ctivircy Cal wom 
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aligned current along the poleward boundary and this 
diverges equatorward in the ionosphere, across the 
electrojet, and up field lines at the equatorward boundary 
of the electrojet. In regions of eastward ionospheric 
current, the north-south current system is reversed. 

As well as these large scale currents, two other 
current systems are involved in this global model. Rostoker 
and Hron (1975) have shown that there is a weak eastward 
current equatorward of the main westward electrojet in the 
dawn sector. -This current system is believed to bea 
Signature of electron precipitation associated with a 
reconfiguration of the magnetosphere after a period of 
Stcong activity. On a Statistical basis then, the “strength 
of _this current system would average out to be relatively 
weak, and it has been included in the model only as a very 
weak current system. An additional current system which may 
‘be; Similar to Ss or DP-2 has also been added to the model 
and the reasons for including this system will be discussed 
at length later in this chapter. 

Current in the polar cap (the region from the north 
pole to the poleward boundary of the auroral oval) has not 
been included in this model. All the magnetic data which has 
been used in formulating the current model is winter data, 
andy durings this. period, thes polar cap is essentially dark 
throughout the day. Figure 4.3 shows the sunlit-dark 
terminator for 1900 UT at its extreme positions during the 
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Figure 4.3 


The extreme positions of the sunlit-dark 
terminator for the period from Day 33225 1978 
to Day 24,<1972, at) an altitude: Gf 417 .meeee 
The solid lines oriented from left to right 
indicate the terminator, and the dotted line 
is the position of the poleward border of 
the model auroral oval. The approximate 
location of the magnetometer line is shown 
for 1700. UT, and (they -Llocet20n  -Oneere 
Greenwich Meridian is shown at 0° longitude. 
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if any solar ultra-violet radiation reaching the polar cap 
upper atmosphere and therefore there is little if any photo- 
ionization. This in turn means that there is little if any 
polar cap conductivity, and therefore no significant polar 


cap current. 


4.2.3 The Electric.Field and Conductivity Model 

Mozer and Lucht (1974) and Iversen and Madsen (1977) 
have published data on the auroral zone electric field. When 
the research for this thesis was initiated only the former 
data were available and these data were used as the electric 
field model required in modelling the currents. The electric 
field consistent with the giobal current model aS 
Shown in figure 4.4. This polar plot shows only the unit 
vectors of the electric field because to define the electric 
ficid explicity, one requires an explicit. model. of the 
COMGUCLIVItyY. Thisiwil) bevdcalt within detarl “in section 
4,3. 

Figure 4.4 is most easily compared with the data of 
Iversen and Madsen (1977) (Figure 3.13). In the region of 
the auroral zone the model electric field is in reasonable 
agreement with the real data. Note that the model electric 
field immediately equatorward of the auroral zone in the 
dawn sector is essentially eastward. This is in keeping with 
the work of Rostoker and Hron (1975) in which it is shown 
that the dawn eastward electrojet is a Pedersen current. The 


east-west component of the low latitude electric field is 
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Figure 4,4 Polar plot of -the -unit vectors! (ota. 
electric field used in the auroral current 
fiow model. 
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directed towards noon in all instances, and the north-south 
component is equatorward in the morning and poleward in the 
afternoon. This is qualitatively in agreement with the ion 
drift velocity measurements of Heelis et al (1976). As well, 
engrelectticijaaelid, of «this iconfiguration;! for a matio of 
height-integrated conductivities (Hall to Pedersen) of 2, 
drives currents which are similar to the Sy systen. 

The conductivity was not absolutely defined in this 
model. However, the ratio of the height-integrated Hall 
conductivity to Pedersen conductivity was taken to be 2 
everywhere. Brekke et al (1974) have shown that for low 
Magnetic giactivity,. this!- ratios holds.in the:auroral zone. 
Yasuhara et al, (1975) Similarly have used this same height- 
integrated conductivity ratio in the auroral zone and _  sub- 
auroral latitudes. Since the current model developed for 
this thesis is a model for average conditions, the ratio of 


2 is well justified. 


As described in the previous chapter, a basic concept 
of what east-west and field-aligned currents were flowing 
had been established and it was also pointed out that there 
exist north-south currents which significantly influence the 
Magnetic perturbation pattern. The approach taken in 
developing this model was to generate magnetic perturbations 
from all currents believed to be present in the ionosphere 


and magnetosphere and combine these to produce a total 
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percturpatron [lLeld.) Thie- currents could not ali be Chosen 
andependently of one “another; that is, current continuity 
had to be maintained. This approach is similar to that used 
by Kisabeth (1972) to model substorms, but it is unlike that 
Moea DY Vdsunard ct al (1975) in which “a “global” model” of 
height-integrated Hall and Pedersen conductivities was 
assumed, along with a field-aligned current model, to solve 


the equation 


eal rhe 


where [ is the ionospheric height-integrated current 
density and Jy the field-aligned current density. 

In this thesis, J in the east-west direction has been 
assumed, and J, has been derived to afford current 
continuity. To determine the north-south component oy al ta 
the electric field model and ratio of height-integrated 


conductivities has been invoked. Consider a coordinate 


system in which x is directed northward, y eastward and zZ 


down . The horizontal height-integrated current density is 
given by 
6*E 
HE aS ag a B 4.3 


where Pas and paps are the height-integrated Hall and 
Pedersen conductivities respectively, E is the horizontal 
Srectr iG field Vector, and 8 the magnetic induction field 
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Dividing equation 4.5 by 4.4 


SEiage 1 - (275 \(%e) sail a via 


Thus, if one component of the horizontal current is assumed 
known, then one need only specify the ratio of the electric 
field components and the ratio of the height-integrated 
conductivities. The value of Ey /E, is the tangent of the 
angle that & makes with the x-axis, and because of this, 
only unit vectors have been shown in figure 4.4. 

ie skouldjbe spointved, out dthatesince the ylocus of j yeast= 
west current flow is not along parallels of latitude (see 
equation 4.1), before equation 4.6 may be used, a correction 
must be applied to the east-west height-integrated current 
density. That is, the auroral oval has an azimuthal 


deflection angle (Mj from a latitude circle given by 


-1/6,-0, .; 
«x = tan (es sin) ae 


and care must theretore be taken in specifying the true 
east-west current density. 

When, from the round based data, it was determined 
that the east-west current was varying in longitude, this 
Variation ‘has been ‘assumed, (to. be, linear. There jis no 
evidence.,available. concerning, the exact nature ,.of the 
longitudinal variation of the convection electrojets, and so 
the ese variation possible,’ but . one “Sstillvebie to 


reproduce the observations, has been chosen. 
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In the model developed for this thesis, ionospheric 
GULTENt 1s assumed "to flow “in a Sheet of “infinitesinal 
vertical extent; that is, the current is height-integrated. 
One must therefore decide at what altitude this ionospheric 
current sheet must be placed. Kamide and Brekke (1977) have 
examined radar data and found that cn average, the westward 
electrojet is situated at 120 km, and the eastward jet at 
100 km. For both electrojets, Kisabeth (1972) has used the 
value of 115 km. In this thesis, it has been assumed that 
the ionospheric sheet currents fliow at an altitude of 115 
km. Any error in the results due to an incorrect value of 
the current altitude will be less than 10% (Kamide and 
Brekke 1977). 

To produce the most reliable quantitative current 
model, it is desirable to include the effects of currents 
induced in the earth. AS discussed in the previous chapter, 
induction ° effects have been included by placing an 
infinitely conducting sphere at depth. For calculation of 
model substorm fields, Kisabeth (1972, 1975) and Kisabeth 
and Rostoker (1977) used a depth of 250 km, based on the 
fact that this depth gives the best fit between model 
Calculations and “realOfdata.. In| the modellang for > this 
thesis, a depth of 600 km has been used throughout. This too 
as tlaniNéenpiricateresiit, giving thesbest qualitative fit of 
the model to the data, when cases dominated by pure overhead 
Picessouet flow!are* considered. Indeed, since this work” is 


concerned with only slowly varying pertubations inasmuch as 
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all the data is hourly averaged, one might expect that the 
induction effect would be less than for the case of substorn 


nodeliing. 


In this section, the main current systems employed in 
the model will be described at some length. Basically there 
are 5 main currents; i) the westward flowiny current, which 
flows from noon through midnight to dusk; ii) the eastward 
flowing current, which flows from noon to midnight; iii) the 
north-south current associated with each of the above east- 
west currents; iv) the current which flows equatorward in 


the morning sector and poleward in the afternoon sector, and 


v) ‘an eastward current in the dawn sector. 


4.4.1 The Westward Current System. 

Getewail Dn bes smecaililed © £fromAcchaptervsthat. tthe rground 
based data for the time period noon through dawn to 
approximately midnight had been interpreted basically as a 
growing westward current. However, near dawn (approximately 
0800 MLT for winter) there was a marked decrease in the 
Magnitude of the peak of the superposed epoch Ax! profile. 
It is not perfectly clear what causes this decrease, but the 
following is offered as a possible explanation. 

Currents “‘flow © in “the auroral ‘oval because this region 


of the atmosphere has an enhanced electrical conductivity 


relative to the rest of the atmosphere. This conductivity is 
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brought about by the precipitation of charged particles from 
the magnetosphere which, through collisional processes, 
produce charge carriers in the ionosphere. As well as this, 
where the upper atmosphere is sunlit, solar ultra-violet 
(UV) radiation may produce charge carriers in the atmosphere 
through photo-ionization processes. Thus, up to dawn (and 
Symmetrically, up to dusk), there are two sources of 
conductivity, solar UV and particle precipitation. Away from 
dawn (and dusk), towards midnight, only particle 
precipitation is availabie TO produce electrical 
conductivity so that near the dawn and dusk terminators 
there is a gradient in conductivity (probably over a narrow 
region). AS well, there is a gradient in the electric field 
near dawn and dusk, (see Figure 3.13(a)). in order to 
accommodate these gradients in the conductivities and the 
electric field, current flows up field lines. Thus there is 
‘a reduction in the magnitude of the ionospheric current, and 
a corresponding reduction in the peak of the X*'t-component, 
as observed, 

In the global current model, it has been assumed that 
both the Hall and Pedersen conductivities vary in such a 
Manner that their ratio remains constant. In section 4.6.1 
Or thisrchapter, the effect of this reduction in’ .electroyet 
strength will be examined in more detail. Figure 4.5 is a 
PhotsoL current sntensity as) a function of time) foryboth ‘the 
eastward andsweetward. iclectrojets., It is, seen “that. the 


westward jet (the curve with values everywhere less than 
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Figure 4.5 


The model height-integrated current density 
for the east-west currents, as a function of 
NLT. The . current densities £0r botheew. 
eastward current (values > 0) and the 
westward current (values < 0) are Shown. . 
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zero) exhibits variations over its entire length. From noon 
to about 0800 MLT, the current density increases linearly 
PoOmO. CetOr U.221.-ANmt. on this ‘sector,) “the ‘width vis 
constant at 5° of latitude, so the current here ranges from 
0.0 A to 1.25x105A. Near 0800 MLT, there is a 20% reduction 
in the total current, with a concomitant reduction in the 
height-integrated current density. As a result, 2.5x10%A is 
fed out of the ionosphere in a field-aligned current sheet. 
This field-aligned current is most likely distributed over a 
finite longitudinal range, but this range cannot be 
determined from the available data. In lieu of this, the 
upward flowing current is replaced by an equivalent upward 
flowing current sheet. 

From 0800 MLT to 0100 MLT, the westward current 
continues to grow back to its peak value in the sunlit 
sector. That is, the total current increases from 1.00x105A 
sau OS00 MLTeto 1,25x105A at G100 NLT. 

Over the entire region discussed so far, since the 
westward current is growing, there must be a downward field- 
aligned current feeding the westward jet, (apart from the 
upward current sheet at 0800 MLT). Since the ionospheric 
current grows most rapidly in the sector from noon to 0800 
MLT, this is obviously where the greatest inward field- 
aligned current flows. Figure 4.6 is a plot of field-aligned 
current density as a function of magnetic local time. It is 
clear chat between midnight and 0800 MLT, very little field- 


aligned current flows into the ionosphere. This field- 
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Figure 4.6 


The model unbalanced field-aligned current 

densities that connect with the east-west 

current systems, as a function of HLiae 

denotes downward field-aligned current, and 
© upward field-aligned current. 
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aligned current which feeds the electrojet is a local net 
field-aligned current inasmuch as if one were to integrate 
along a meridian, there is an imbalance in the field-aligned 
current flow. 

Keferring to. figure 4.5. again, it is seen that the 
westward current intensity generally decreases from 0100 MLT 
fOn1000 HLT. fn fact, thertotal ‘current; | flowing ‘across a 
meridian in this sector decreases linearly from 1.25x105A at 
0100 MLT, to 0.0 A at 1800 MLT. The non-linear variations in 
the current intensity come about from variations in the 
width of the westward current in this regime. In this 
sector, field-aligned current must *ilow .out tof WwW the 
ionosphere, and the outward field-aligned current density is 
Shown by the open circles in figure 4.6. The variations in 
Pas  GCULVe! -are)) again’ \;due')to the varying width of the 
current. In fact, the upward flowing current as a function 
of longitude is constant. The high current density at 1800 
MLT is due solely to the extremely narrow (0.5°) latitudinal 
extent of the electrojet at this point. 

Superimposed on this entire westward electrojet system 
is a north-south current system. In accordance with the 
procedure outlined earlier, the current flowing mii netne 
north-south system was obtained directly from the electric 
meld and conductivity models. Figure 4.7 isva plot of the 
global north-south current system plotted in terms of the 
Piet alcaned current contribution. The curves of immediate 


interest are the upper two curves in the pre-midnight hours 
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Figure 4./7 


The field-aligned current strength for the 
north-south current systems as a function of 


ML?. © denotes upward field-aligned 
current, and @ indicates downward 
field-aligned current. fhe dashed line 


represents current flow at the equatorward 
border of the auroral oval, and the solid 
line indicates current flow at the poleward 
border of the auroral oval. 
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and their continuation into the post-midnight hours. Note 
that the downward current density is greater than the 
upward. This is due to the fact the southward flowing 
current flows along meridians and therefore diverges, so 
_that for a constant current, the current density varies as 
the sine of the colatitude. 

The 24 panels of Figure 4.8 are the model latitude 
profiles for the westward flowing currents, including the 
associated “net" field-aligned currents, but excluding the 
north-south currents. The features that one might expect 
from knowledge of the current distribution are evident in 
these profiles. There is a decrease in the X*f-component 
extrema from 0900 MLT to 0800 MLT (Figure 4.8(j) and Figure 
4.8(1i) respectively). Because of the presence of unbalanced, 
downward field-aligned current from 0100 to noon MLT, one 
would expect a positive going level shift (north-to-south) 
in the east-west component. This in fact is seen only up to 
0300 MLT; the model profiles for 0100 and 0200 MLT show 
@ither a weak negative going shift in AY‘, or no shift at 
all. This is due to the strong upward field-aligned current 
in the midnight sector, completely masking the weak downward 
current in the post-midnight sector. Indeed, the positive- 
going level-shift is seen beyond noon, spilling over as an 
end-effect into the 1300 MLT and perhaps 1400 MLT sector. 
However, the magnitude of all these level-shifts appears 
small relative to those suggested by the real data. fhis 


could be due to one or both of two possibilities: (1) there 
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Figure 4.8 


Model latitude profiles due to the complete 
current system (including the 
ionospheric and 
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Figure 4.8 Model latitude profiles due to the complete 
westward current system. (including the 
ionospheric and field-aligned currents). 
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Figure 4.8 Model- latitude profiles due to the complete 
westward fou) Pap ax =} thn « system {including the 
ionospheric and field-aligned currents). 
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Figure 4.8 


Model latitude profiles due to the complete 
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is more unbalanced field-aligned current than has been used 
in the model, or (2) there are other Ouewene systems which 
will produce a AY' level-shift. The first of these requires 
an ionospheric current system other than.a westward current, 
since increasing this latter will not change the size of the 
Y'-component level-shift relative to the X'-component 
extremum. Since, aS pointed out earlier, it is not believed 
that significant currents flow in the winter polar cap, this 
first possibility suggests the presence of currents flowing 
equatorward. Such a current system has been included in the 
global model, and will be returned to later in this chapter 
Mee tion 4.4.3). The second possibility is also a viable 
one, insofar as the north-south current system has not yet 
been considered. 

First, however, the remaining features of the westward 
system aione should be considered. Note that the Z-component 
‘exhibits an asymmetry not unlike that in the data, in the 
post-midnight hours. Also, the east-west component is not 
negative across the entire profile, as it is in the data, in 
the midnight sector. Finally, the east-west component does 
show a negative-going step, in the pre-midnight hours, 
Similar to that observed in the data. 

When the appropriate north-south currents are added to 
cis westward current system, the model profiles are 
im proved considerably when compared: with sthe realy data. 
Figure 4.9 (24 panels) are plots of model profiles in which 


the profiles of figure 4.8 have been added to the southward 
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Figure 4.9 Model latitude proziles due to the complete 
westward current system, to which has been 
added the corresponding north-south current 
systen. 
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Figure 4.9 Model latitude profiles due to the complete 
westward current system, to which has been 
added the corresponding north-south current 
Systen. 
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Model latitude profiles due to the complete 
westward current system, to which has been 
added the corresponding north-south current 
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Figure 4.9 Model latitude profiles due to the complete 
westward current systen, to which has been 
added the corresponding north~south current 
system. 
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flowing meridional current system consistent with the 
electric field and conductivity model. In the post-midnight 
sector, to noon, the size of the Y'-component level-shift is 
enhanced. As will be demonstrated later in this chapter, 
this levei-shift is now in close agreement with the data. 
Thus, as suggested in the penultimate paragraph, the 
possibility of additional currents producing a AY' level- 
shift is a reasonable one. Although level-shifts in the 
east-west component may pe due in part to unbalanced field- 
aligned current, they cannot be explained totally on this 
basis; a north-south current system of great longitudinal 
extent will also produce a similar Y'-component. However, it 
is not possible to model the Y't-component levei-shift solely 
by . north-south current flow, and still produce the observed 
Variation in the X*-component. AS already discussed, it is 
evident that the westward electrojet grows over a region of 
some longitudinal extent, and this demands the presence of 
unbalanced field-aligned current into the ionosphere. 
Significant end effects of the north-south system 
associated with the westward electrojet are evident in the 
model profiles for times immediately after noon (Figure 4.9) 
when they are compared with the corresponding panels of 
Figure 4.8. This effect is seen mainly in the 4&-component, 
and is due primarily to the ionospheric portion of the 
north-south three-dimensional system. A Similar effect may 
be seen in the 1800 MLT model profiles, (Figures 4.8(s) and 


4.9(s)) and is due to the north-south current systen 
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associated with the pre-midnight westward electrojet. 


4.4.2 The Eastward Current System 

In Chapter 3, it was shown that the ground based 
magnetic data in the post-noon sector were indicative of an 
eastward flowing electrojet. In fact, based on the magnitude 
of the X'-component positive extremum this eastward current 
was found to grow from noon to near dusk, where there is a 
small decrease in the apparent current magnitude. This is 
Similar to the behaviour of ene westward current flow, 
although this decrease in current strength is not as 
pronounced. However, the same argument as developed for the 
westward current has been applied to the eastward current in 
the post-noon sector. AS well, there is a marked decrease in 
the magnitude of the positive AX' extrema from about 2000 
Brrecoe2200 ALTERS Figure 984.5) "is. ‘a * plot “of “the height= 
‘integrated current density for the eastward electrojet (the 
curve with positive values in this figure). The current 
grows linearly from 0.0 Am-? at noon to a maximum of 0.116 
An=05(6257xT0*A“totaljl at 1600 °MLTS At this time, “there tsa 
reduction in current, with 20% of the current flowing out of 
the ionosphere in a field-aligned current sheet. The current 
remains constant until 2100 MLT, at which time it begins to 
flow poleward. As can be seen in figure 4.1, in this regime 
(2100 to 2400 MLT) the eastward current is constrained to a 
wedge-like region of the auroral oval. As the eastward 


flowing current diverges poleward, it encounters the 
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equatorward border of the westward current, where the 
eastward current diverges up the field lines. Thus, in the 
model, the eastward electrojet flows out of the ionosphere 
in a field-aligned current sheet, the locus of which 
Separates the eastward and westward electrojets in the 
region of the Harang discontinuity. Figure 4.6 shows the 
current density of the field-aligned current which feeds the 
eastward electrojet (1200 to 1600 MLT). It is considerably 
less then the main feed to the westward electrojet (0900 to 
1200 MLT), in keeping with the relatively smaller total 
current that flows in the eastward electrojet. As in the 
case of the westward flowing current, there is a north-south 
current system associated with the eastward electrojet. In 
the post-noon sector, this current is directed poleward in 
keeping with the electric field observations. Figure 4.7 
includes a plot of the field-aligned current density 
associated with this north-south current. Although this 
current density is considerably less than that used for the 
westward current, it is consistent with the electric field 
model and conductivity ratio described earlier. Figure 4.10 
(24 panels) are plots of the model latitude profiles for the 
eastward current together with its associated fieid-aligned 
Currents. The growth of the eastward electrojet 1S apparent 
in the post-noon profiles up to 1600 MLT (Figures 4.10 ({m) 
through 4.10(¢qg)). The slight decrease in the peak value of 
AX' is visible at 1600 MLT, and from 1600 MLT to about 2100 


MLT (Figures 4.10(q) to 4.10(v)) the magnitude of the X'- 
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Figure 4.10 Model latitude profiles due to the complete 
eastward current system (including the 
ionospheric and field-aligned currents). 
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Figure 4.10 


Model latitude profiles due to the complete 
current system {including the 
ionospheric and 
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Figure 4.10 
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Figure 4.10 Model latitude profiles due to the complete 
eastward current System (including the 
ionospheric and field-aligned currents). 
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component peak is essentially constant. Beyond this time to 
midnight, the strength of the eastward electrojet decreases. 
piceecrtt) Of; the current, flow jwith respect todparallels iof 
datitude that has been built into the model is evident in 
Bercy, model” profiles, for 1700  MLT to 1900 J/MLT.. (Figures 
4.10(r) to 4.10(t)). In this time period, the tilt in the 
Oval ais greatest, and in these profiles, it is clear that 
the east-west component follows the north-south component. 
The east-west component in the immediate pre-midnight sector 
(Figures 4.10(a), 4.10(w), 4.10(x)) Shows a negative going 
level-shift (from south to north latitudes) consistent with 
the presence of the upward field-aligned current sheet in 
this regime. 

Outside the regime of the eastward current, some 
effects are visible. For example, there is a small, broad, 
positive AX' signature seen in the model profile for 0000 
meLt (Fig. 4.21 (a)). Similarly, an end-effect is evident in 
the two hours preceding nocn (Figures 4.70(k) and 4.10(i)). 
In general however, the eastward flowing current produces 
little effect in regimes outside the actual location of the 
current. 

The result of adding the north-south current systen 
which flows in conjunction with the eastward current is 
shown in Figure 4.11. The effect of the system on the 
profiles shown in Figure 4.10 is not as large as was seen in 
the case of the westward electrojet. The main effect is 


observed in the profiles for 2100 MLT through midnight to 
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Figure 4.11 Model latitude profiles due to the complete 
eastward current system, to which has been 
added the corresponding north-south current 
systen. 
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Figure 4.11 Model latitude profiles due to the complete 
eastward current system, to which has been 
added the corresponding north-south current 
systen. 
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Figure 4.11 Model latitude profiles due to the complete 
eastward current system, to which has been 
added the corresponding north-south current 
systen. 
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Figure 4.11 Model latitude profiles due to the complete 
eastward current system, to which has been 
added the corresponding north-south current 
system. 
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2300 MLT (Figure 4.11(v) to 4.11(x). In Figure Gotucvy, (2100 
NLT) the Y‘'-component is altered by the poleward flowing 
current system so as to almost completely mask the effect of 
the tilted eastward current system. The same effect is seen 
in Figure 4.11(w) (2200 MLT), where as well, the level-shift 
in the Y'-component is reduced by the addition of the north- 
south current system. Also, in the 2200 MLT profile, the 2Z- 
component is changed so that the well-defined negative 
extremum shown in the profiles for the eastward current 
alone (Figure 4.10 (w)) is all but lost. The profiles for 
2300 MLT (Figure 4.10(x) and 4.11(x)) show that the north- 
south current system produces a marked change in all 3 
components; the X'-component develops distinctly positive 
warues, at high Jatitude;, the Y'=-component ‘takes on the 
character of a strongly tilted east-west current system; and 
the Z-component gains a marked positive extremum. Apart from 
amis tine period, littie change is observed in the profiles 
which result from adding the poleward flowing current systen 


to the eastward electrojet. 


4.4.3 The North-South Current System 

For the sake of completeness, the entire north-south 
current system is shown in Figure 4.12 (24 panels). This 
system is frequently neglected in auroral current models, at 
least in the case of polar magnetic substorms. Indeed, in 
these cases there may be some justification for considering 


only the east-west current systems inasmuch as the ratio of 
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Figure 4.12 Model latitude profiles due to the entire 
north-south current system. — 
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Figure 4.12 Model latitude profiles due to the entire 
north-south current systen. 
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Figure 4.12 Model latitude profiles due to the entire 
north-south current system. 
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Figure 4.12 Model latitude profiles due to the entire 
north-south current systen. 
S)utSOG BU Om te 00 wha 
op pee LS HORG) Geley PAO MMOR IER 
Uu)-200 0. TomZtIO0eN Ge 
Vp U2100 C0 22200 GALT 
wp 2200 CO 3 00 SE 
x) P2300 tC oOreeg Gena T 


2CHR O MIN. O SEC. 


FIELD INTENSITY (NT) 


'SS 60 65 70 75 80 85 
CENTEREO OIPOLE LATITUDE 


1SHR O MIN. O SEC. 


(NT) 


*eses. 


eee, 


FIELD INTENSITY 


‘SS 60 65 70 75 80 85 
CENTERED DIPOLE LATITUDE 


18HR O MIN. O SEC. 


ELELD INTENSITY (NT) 


SS 60 65 70 75 80 85 
CENTERED OIPOLE LATITUDE 


ETECO, INTENSITY (NT) 


(NT) 


FIELD INTENSITY 


(NT) 


FIELD INTENSITY 


207 


23HR O MIN. O SEC. 


ISS 60 65 70 75 $0.85 
CENTERED SOIPOLE LATITUDE 


22HR O MIN. O SEC. 


SS 60 6S 70 75 80 85 
CENTERED DIPOLE LATITUDE 


21HR O MIN. O SEC. 


es 
‘SS 60 65 70 75 80 8s 
CENTEREO DIPOLE LATITUDE 


rage Hed oe net eld 4 be 


Res 


len th : a he 
Hig) THY ih Saye ‘CASE ADT 
Utd TH. vere to oe Bi 


ah | 
§ ay Sago ta fp nee’ 


itn rm “ny, “ es Rats 
ay } ‘ 


Higgnrgen gel opm eter 


Fi ' ; a : ene +A hye ‘ a it ae ' 
' r eh 
ree . ae 018 Li ne 
, . By a aint * 4: re 
i 
| ; i = . 
r ar o*. sy 
TYP Vntaedy (te 


| ' 
F L; A fi 
' ‘, 
! " \ 
1 res “ 2s 70 t x F 
H : : ‘ a ae aabeale Dae 
ee Se | % 
aaives 2. i “ni, swe ‘ Mite 
gees . ¥ we 
a e 2 “4 Esa = 
F vat * ia alt ak eo, + enen siete baeita ie 
‘ haere ian , y ; 
a Pe eae ae ve | NA: ion 4 
| } s 
| 3 
er » Ka 
ae 
| e 
> i ‘ if 


a e ; 
4 pom pols » — re, ? rh - 
- ¥ an ond Vales 


208 


aly to oe increases to about 4 during periods of high 
Magnetic activity (Brekke et al, 1974), so that the role of 
the north-south currents is relatively less important during 
Substorms. Since the north-south currents are primarily 
Pedersen currents, their role in substorn modelling is 
possibly relatively small. Kisabeth (1972, 1975) and 
Bannister (1976) have successfully modelled substorms while 
ignoring north-south current systems. However, in this 
model, representing quiet to moderately disturbed periods 
the north-south current system becomes relatively important 
as May be seen in the model latitude profiles of Figure 
4.12. In the first two profiles (Figure 4.12(a) and (b); 
0000 and 0100 MLT), the effect of the northward flowing 
current in the midnight sector is seen as a strong positive 
Zao perturbation. “' It) will ‘be “recalled that in the “SPEA 
profiles from the post-midnight quadrant, the Z-component 
‘exhibited an asymmetry, with the positive extremum being 
greater than the norm of the negative extremum. In fact, 
this positive AZ profile continues in the model north-south 
System until about dawn. From 0400 to 1100 MLT (Figures 
Geriz tej to-- 4s 12 tly ys “the**yY*-component™ *is* “ot greatest 
Gnterest. Tt is’ evident that the total north-south "system 
enhances the level-shift in AY‘, as was seen in the 
comparison of Figure 4.8 and 4.9. In the 1100 and 1200 MLT 
model profiles (Figure 4.12(1) and (m)), the Z-component is 
distinctly negative. As will be seen in section 4.5, this 


effect is highly influential in the total model. As already 
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noted, the north-south current system has only a minimal 
effect in the region of the eastward electrojet, from noon 
to dusk. However, the total system produces a fairly marked 
effect, in the .post-dusk '.sector, on” the “ZN Yt profiles 
guces. fe I2ct) to. 4.12 v)<...1900 . to.4 2106) MLE ele this 
regime, the. total north-south current system markedly 
enhances the negative going level-shift in the Y'-component. 

Finally, the last two panels of Figure 4.12 (w and x; 
2200 and 2300 MLT) show a negative Y'-component across the 
entire latitude profile, in keeping with the poleward 


flowing current in this region. 


4.4.4, Other Contributing Current Systems 

In addition to the current systems discussed up to this 
poant, it was necessary.to add two other current. systems, 
although these are of relatively minor importance. The first 
of these was the eastward. flowing current located 
immediately equatorward of the westward electrojet in the 
daw. Sector. This current was.made to flow anti-parallel. to 
the westward electrojet in a 5° wide latitudinal strip whose 
poleward border coincided with the equatorward border of the 
westward electrojet. The current was limited to the 
hong tucessrange (or .~50° +0.1659; that is; ,castwand curpent 
flowed trom 0400 MLT tov 1100\MLT. For modelling purposes. it 
was assumed that this current is fed from a downward field- 
aligned sheet current that diverges eastward. At dawn, the 


conductivity gradient described earlier (Section 4.4.1) is 
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encountered, and additional Current fiows into the 
Pon ose tere at this point. At 1100 MLT the current diverges 
back up a field-aligned current sheet to the magnetosphere. 
The maximum current flowing in this electrojet is 10¢ aA, for 
a maximum height-integrated current density of 0.018 Am-!,. 
Latitude profiles for this current system possess no unusual 
features and are not presented here. The maximum positive 
AX* perturbation is about 5nT, occurring in the profile for 
0900 MLTI. 

The final current system is also a relatively weak 
System when compared with either the westward or the 
eastward electrojet systems, but one which possesses useful 
features, from the point of view of the global model. 

When this study was first undertaken, it was observed 
that the AyY' level-shift was not symmetric about the field 
intensity origin. (See, for example, Figure 4.3) However, an 
‘unbalanced, distributed field-aligned current produces a 
AY' level-shift which is essentially symmetric about the 
Grigin (figure 3.77). As well, it became ‘evident that~ the 
inferences drawn from the ground-based magnetic data 
regarding unbalanced field-aligned currents were not in 
total agreement with the results obtained from the Triad 
Magnetometer data (Sugiura and Potemra, 19/6). That is, 
interpretation of the satellite data showed that net upward 
field-aligned current flowed in the post-noon sector, the 
exact opposite of the results of the ground-based data. 


Hughes and Rostoker (1977) offered an heuristic argument in 
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an attempt to explain this discrepancy. This argument is as 
follows. Rostoker and Hron (1975) have pointed out that in 
the dawn sector there exists an equatorward-flowing Hall 
current, and it is assumed that this Hall current circulates 
through the low latitude ionosphere to the post-noon 
quadrant and then flows up field lines into the 
magnetosphere. Consider two infinite anti-parallel plane 
sheet currents separated by a finite distance (Figure 
4.13(a)). The ionospheric equivalent current LOL the 
downward-flowing sheet current (I) is two horizontal sheets 
Ores strength 5/2 flowing: toward the vertical sheet. 
Similarly, the equivalent current for the upward-flowing 
current sheet (I) is two horizontal current sheets (1/2) 
flowing away from the vertical sheet. (This is simply an 
extension of the technique used by Fukushima (1969) for line 
currents). Connecting these two vertical sheets by a 
‘horizontal sheet current results in an equivalent current 
for the system of zero, so that no magnetic perturbation is 
observed on the ground. 

If, as suggested above, there is a horizontal current 
flow into the region of Birkeland current flow associated 
With the auroral oval in the post-noon quadrant (Figure 
4.13(b)), then there will be an unbalanced upward current 
flow (as observed by friad), but no level-shift will occur 
in the east-west component as observed on the ground. 
However, Figure 4.13(b) shows that there will be a bias in 


the east-west component. A similar argument to this may be 
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Figure 4.13(a) 


Schematic drawing of equivalent ionospheric 


currents LOL infinite vertical current 
Sheets and a pair of anti-parallel vertical 
current sheets connecting through a 


horizontal sheet. Equivalent currents from 


the downward-flowing current are shown by 


the dashed lines, and those from the 
upward-flowing current are shown by the 
dotted lines. Real currents are shown by the 
solid lines. 
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Figure 4.13 (b) 


Schematic drawing of equivalent 
Currents for 
upward current (1, 


downward current, LOLS because of 
ionospheric current ‘I, flowing into 
system. The equivalent “current ) 172 
uniform across the current systen 


directed to the right. 
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applied to the pre-noon quadrant to produce a bias in the 
east-west component. That is, downward current which 
diverges equatorward will preduce a positive bias in the Y'- 
component. Thus, it is possible to explain the discrepancy 
between the ground data and the satellite data, as well as 
explain the observed bias in the Y'-component, with the same 
apqument.s "Of “cowEse; yin reality; «this arguments is not 
totally valid. The closure of the vertical currents is not 
complete, since closure currents in the outer magnetosphere 
are neglected, and also, the real currents are of finite 
extent. However, the argument was successful enough to lead 
to further investigation of such a current systen. 

These currents will be referred to as the low-latitude 
current system. The total current flowing in this system is 
0.1875 x 105A, only 15% of the maximum westward flowing 
current. The low-latitude current system consists of field- 
‘aligned current flowing into the auroral oval within the 
0600 to 1100 MLT sector. This current diverges equatorward 
and eventually into an eastward flowing current located 
between 0° to 40° latitude. In the post-noon quadrant, this 
eastward current in turn flows poleward and then diverges up 
field lines at the poleward border of the eastward 
electrojet. Note that in the pre-noon sector, the field- 
aligned current is distributed over the latitude range of 
the oval, whereas in the post-noon sector, the field-aligned 
Current is confined to a sheet. From 0800 to 1100 MLT, the 


ionospheric current has an average height-integrated current 
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density of 6 x 10-*Am-1, from 0800 to 0600 MLT, the height- 
integrated current density decreases to zero. At sub-auroral 
latitudes, the ionospheric current intensity is anti- 
syMmetric about noon. 

Schematically, this current system bears some 
resemblance to the Sq current system. However the model 
system is not completely closed in the ionosphere, as is the 
Sq system. This model system bears an even closer 
resemblance to the DP-2 system of Nishida (1968) although 
DP-2 is not symmetric, in its geometry, about noon. Akasofu 
et al (1973) have suggested that Nishida's DP-2 system may 
be due to an intensification of the Sy current system 
together with an expansion of the auroral oval (see Chapter 
1). It is possible then that the low-latitude current systen 
corresponds to DP=-2 or SY for the winter ionosphere, since 
Oy polar cap current flows in the model. This topic wili be 
‘discussed somewhat further in section 4.7.1 of this chapter. 

The 24 panels of Figure 4.14 are latitude profiles 
generated by the model low-latitude current system. It is 
evident from Figures 4.14 (a) to 4.14(f) and Figures 4.14(t) 
POmeletu(xjeu thats whis System has little etiect an tnevdark 
hemisphere. Note that in the 0600 MLT to 1800 MLT sector 
(Figure 4.14(g) to 4.14(s)) the Y'-component is biased as 
described in previous paragraphs, i.e. AY' is positive in 
the pre-noon sector, and negative in the post-noon sector. 
However, unlike in the theoretical argument put forth 


earlier, there is an additional feature to the Y'-component, 
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Figure 4.174 Model latitude profiles due to the 
low-latitude current systen. 
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Figure 4.14 Model latitude profiles due me the 
low-latitude current system. 
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in that the profile does not exhibit only a bias, but shows 
some structure as well. However, AY' in these profiles is 
not easily interpreted in terms of unbalanced field-aligned 
current, so that in fact, as suggested by Hughes and 
Rostoker (1977), it is possible to have oe unbalanced field- 
aligned current without there being a readily recognizable 
Signature on the ground. In spite of this, the strength of 
the upward fieid-aligned current in the post-noon sector due 
to the low-latitude current system is insufficient to 
produce an overall net upward field-aligned current, so that 
the discrepancy between this model and the Triad 
observations remains. This will be dealt with further in 


section 4.7.1. 


4.5 The Global Current Model 


Figure 4.15 is a polar plot showing sSchematicaiiy the 
-current systems used in this model. For the sake of clarity, 
the north-south current systems are not shown here. Neither 
relative widths of the arrews representing ionospheric 
currents nor the diameters of circles representing 
unbalanced field-aligned current are necessarily drawn to 
scale. However, most of the general features of the current 


systems used in this model are represented. 
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Figure 4.15 


Schematic drawing of the complete model 


current system. Downward field-aligned 
currents are shown by 7 JUpWard “by ° 
The width of the. arrows is’ only Van 
approximate INatCcali1on of relative 
ionospheric current strength. Balanced 
field-aligned currents, and the 
corresponding north-south ionospheric 


currents have been omitted for the sake of 
clarity. 
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4.5.1 Qualitative Comparison of Model Latitude Profiles with 
SPEA Profiles 

Figure 4.16 contains 48 separate panels in which the 
model latitude profile for a given hour is displayed next to 
the SPEA latitude profile for the same hour. 

it is re-iterated at this point that the model has ‘not 
been constructed to produce an exact fit to the data. The 
variations in the data from day to day are large. However, 
certain features do persist from day to day, and it is these 
features that the eerie has been designed to reproduce. 
Thus, a only a semi-quantitative comparison will be nade 
between the model profiles and the SPEA profiles. In section 


and Ay! LOL 


; : . : 
4.4.2, a statistical comparison of AX Hele 


peak 
the model and SPEA profiles will be given. 
Comparison of the 0000 MLT profiles shows a high degree 
of similarity between the model profiles and the data. The 
‘Y*-component is negative across the profile in both cases 
and has a Similar step-like character. The X'-components are 
also similar. The model Z-component is asymmetric with the 
positive extremum being larger than the absolute value of 
the negative extremum, as is the case in the data. However, 
at high latitudes, the observed Z-component maintains 
relatively large values, whereas in the model, the Z- 
component decreases quite rapidly with increasing latitude 
beyond that at which the peak occurs. This is a persistent 


feature in the comparison of the model with the data in the 


post-midnight sector. 
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Figure 4.16 


Comparison of the model latitude profiles, 
for the complete model, with the SPEA 
profiles. Note that the SPEA profiles are 
labelled in Universal Time, and the model 
profiles in Magnetic Local Time. ° Local 


Magnetic midnight is at approximately 0800 
UT. 
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Comparison of the nodel latitude profiles, 
for the complete model, with the SPEA 
profiles. Note \that’ the SPEA protitessage 
labelled in Universal Time, and the model 


profiles in Magnetic: ‘Local Times; Loca 


magnetic midnight is at approximately 0800 
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Figure 4.16 


Comparison of the model latitude profiles, 
for the compiete model, with the SPEA 
profiles. «Note that the: SPEA profiles var. 
labelled in Universal Time, and the model 
profiles in’ Magnetic’ Local ‘Tims. "ilocas 


magnetic midnight is at approximately 0800 
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Figure 4.16 


Comparison of the model latitude profiles, 
for the complete model, with the SPEA 
profiles. Note that the SPEA profiles are 
labelled in Universal Time, and the model 
profiles Zn Magnetic “Local fime.) locas 


magnetic midnight is at approximately 0800 
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Comparison of the model latitude profiles, 
for the complete model, with the SPEA 
profiles. Note that the .SPEA profiles are 
labelled in Universal Time, and the model 
profiles in Magnetic Local Time. Local 
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Figure 4.16 


Comparison of the model latitude profiles, 
for the complete model, with the SPEA 
profiles. Note that the SPEA profiles are 
labelled in Universal Time, and the model 
profiles in Magnetic Local Time. Local 


magnetic midnight is at approximately 08006 
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Figure 4.16 


Comparison of the model latitude profiles, 
for the complete model, with the SPEA 
profiles. Note that the SPEA profiles are 
labelled in Universal Time, and the model 
profiles in Magnetic. Local. Time. "jLogar 


Magnetic midnight is at approximately 0800 
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Figure 4.16 


Comparison of the model latitude profiles, 
for the complete model, with the SPEA 
profiles. Note that the SPEA profiles are 
labelled in Universal Time, and the model 
proitizes in Magnetic. (Local (Times. piocar 


magnetic midnight is at approximately 0800 
Ut. 
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At 0100 MLT,all 3 components in the model profiles are 
in fair agreement with the character of the components in 
the SPEA profiles. 

The 0200 MLT profiles also show good agreement in the 
X'- and Z-components. The Y'-component in the data shows 
evidence of a small positive going level-shift, perhaps 
superimposed on a slight poleward tilted oval. (AAY' appears 
to follow MX"). The model AY profile shows evidence of 
the opposite tilt, as indeed such is built into the model, 
and essentially no level-shift. However, in both cases, the 
Maghitude of the Y'-component is relatively small, and these 
discrepancies are not major cnes. 

For 0300 MLT, the agreement between model and 
observation is quite good. The asymmetry in 4\Z is similar 
in both, and the values of the peak AX' are essentially the 
same in both cases. The Y'-component also shows the same 
‘qualitative features of a slight level-shift across the 
width of the electrojet. 

The fifth set of panels in Figure 4.16 (0400 MLT) also 
show a good qualitive agreement between the Observations and 
the model, although the magnitude of eer in the model 
is less than that observed. This difference is difficult to 
ascribe to the net field-aligned current component of the 
model, as increasing this necessitates increasing the 
westward flowing current. The westward current, however, 
produces peak values of AX' which are in good agreement 


with the observations. It is likely that this difference 
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arises due to the choice of model north-south current. If 
the electric field during the observation period is not well 
modelled by the field shown in Figure 4.4, or if the height- 
integrated conductivity ratio shows variations away from the 
MoneceyOLry 2hiin this sector): <than bie model north-south 
current system would differ from that used. However, lack of 
Specitic~ information: \about »sthese factors“ for ‘the data 
collection period have forced the use of the present north- 
south current model. Further refinements may improve the 


agreement between the observed and model AY'. ,abut?ecCas 


cep 
will be shown in section 4.4.2, the overail agreement is 
acceptable. 

The comments of the above paragraph also apply to the 
comparison of model profiles with SPEA profiles for 0500 MLT 
ponough 0800 "HLT. \Insaddition; it willetbe snoted a that the 
model produces a more step-like Y'-component than is 
‘Observed in the SPEA profiles. Indeed, as described in 
Shapteres3; (thes *nodel se AY profilesPin’ this sector jbehave 
more like those observed in the individual hourly averaged 
latitude profiles than that in the SPEA profiles. 

As the noon sector is approached, it becomes clear that 
the global current model fails to reproduce some features of 
the ‘datac cIt*éwill “be**recalled “that "the “main ‘feature 
described in the noon sector was the large positive-going 
level-shift in the Y'-component. This feature is reasonably 
well produced by the model. However, the magnitude of the 


X'-component is not well modelled by the current systems 


ie ce ve ; ! i) iy ae t i 
Se | eres: marti 


et a 


Pr’ dy \ i eds e: 80 a oni 


aot Pe i. int OR.  coaae’s Asia si § ia 
odenen i aeeaota ce > ot verbaus: al Nee grea: oid i 
el a a Diag | paler » ebro 9 


Aa svt ep 
y "4 4 ’ % fl 4 mn ‘ 
ce an el ee fe yeas a aa Dont inte aA. a rea x 
¢ Lik é es 
ro , rp ee a ist & oye sm was Fe lh ig ih lad ae a 
. i PE vine ee } at Oe 
Pe NE ae eee LOU. lees a Nie Na) 4 A ie ACTF Fete Tk peewee. 
a 


rh AN 
add oo eheie onke ane apackaq anode he no; acta i 


ean Tek at 70% sadex Or ig) éRA2 dase) aakenedee oboe: Wie o 
ve (bebe: re haw: #, perenne. 


y ne 


ae ena NG oe paanedods atin ae Updo 

Hee 
gi: “pod prone stab maiLekoy nate: sone i ‘ano ' 
ieeded “soe aeR SD dy ‘em anatneny +e ilo waist ve “ 


7 


a yal anced iF buninial wee npueandcick!: “a 


: Linen aete aad Wea sau Mis: me | a ruth. 
a at: wrt 
+23, ae eO aangond | fh abaaasorque: ray sos ee “aed er | 


ae 


' (ip aah Ippon: 2930 


vite a past ith ‘alg: deol 


ee 


otaew th autem yaad tase enka pit exe” 
pusourennedeot onde git. ee. aérind) aSea tent, ab 
(ditenor s67; 2h edprne: onan somenagacomty eds ak. a“ 


i 
a 


7 


whe to’ ebusuapee aad) caval utehow: fae ire ri ts ™" 
aadaoya teers ots 7 Sotinbon ane rag i st oe snug nqngoe aE, 


1g ate et yy fo 
. eer yew aa . ele 
, bt : am 
Y Am ole 
c j h igh » 
: oN Ae 


248 


already described. From the superposed epoch analysis, it 
appears that there is little if any electrojet current 
flowing in the noon sector, whereas the end effects of the 
model westward current are evident in the model profiles as 
relatively large MAX" perturbations. ne well, the SPEA 
profiles in the noon sector show evidence of the Sq current 
system at the lower latitudes. That is, AX! is biased 
negatively at low latitudes, consistent with what would be 
expected for the Sq system. Although, as noted in section 
4.4.4, the model low latitude current system might bear some 
Tesemblence\ to the Sq current, Sq circulates in»the high- 
latitude sub-auroral ionosphere as a broad westward current. 
This feature has not been built into the model, and it would 
appear that the atmospheric dynamo current system is 
required to satisfy all the observed perturbations. 

The sinedusionn of isaneSq currentssystenixrilinot soive 
-the entire problem in this sector. From 0900 to 1200 MALT, 
the observed JAX! values across the profile are in general 
less than those in the model profiles, and at noon, the Ax! 
SpeAGr profile sadoes) snotatindicate®) the presence of an 
electrojet. In an attempt to improve the model fit in this 
sector, the model westward electrojet was modified Slightly, 
in that no current flow was permitted from noon to 1020 MLT. 
However, the westward electrojet was to grow to a maximum of 
1.25 x 105 A by 0800 MLT, as in the original model. Thus, 


from 1020 MLT to 0800 MLT, the net downward field-aligned 


current was more intense than in the original model. The 
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effect of this modification was limited to essentially the 
nG@O0; MLT to 1300 MLT sector, as shown in Figure 
Beeiia; b,c; 4). It is apparent that the model AX* profiles 
are improved, although the magnitude of ee is 
increased. It is suggested that the best possible model lies 
between the original and this modification. 

Another possible explanation of this difficulty is the 
fact that the SPEA profiles in the noon sector have been 
constructed differently than elsewhere. At noon, the 
profiles were referenced to AY' (see Chapter 2). It is 
possible that, had the noon sector data been sorted into 
cases which had distinctive eastward or westward electrojet 
Signatures, then the correspondence between model and data 
would have been improved. 

One additional alternative is that there may exist 
current systems in the noon sector other than those which 
‘have been modelled. Recently, Iijima and Potemra (1976) have 
published resuits showing field-aligned currents additional 
to and poleward of those shown in Figure 4.2 in the noon 
Sector, It is not clear at ‘this tine ‘How these ~ currents 
connect, at at all, in ‘the “ionosphere, and TLurther 
observations of the electric field and fieid-aligned 
currents in this region are necessary before conclusions may 
be drawn. At present, recourse to the variability in the 
data, from profile to profile, from this sector must be made 
to minimize the importance of this discrepancy between the 


model and SPEA X'-component profiles in the noon sector. As 
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VOOO MELT to -1300 NLT ‘Sector, athat., pesuLe 
from moving the origin of the westward 
current system from noon to 1020 MLT. 
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Will be shown in section 4.4.2, if this variability is taken 
into account, the overall fit of the model to the data is 
quite acceptable. 

From 1300 MLT onwards, the model reproduces the 
features of the observed profiles quite oe itt The positive- 
going step in MAY' is evident in the model profiles up to 
1500 MLT, and the growth of the eastward current is visible 
through the increasing magnitude of the peak of the X'- 
component. The nature of the Y'-component for the model in 
ene “profiles for 1500, “1600 and 1700 MLT is” not. exactly the 
same as that observed. This is not believed to be a_ serious 
fault in that the model has been computed only at a single 
longitude for each epoch, whereas the SPEA profiles are 
Superposed averages of hourly averaged profiles. Thus a 
boundary which changes from day to day (for example, the 
longitude beyond which the net downward field-aligned 
Gurrent ceases to flow) will introduce certain errors «which 
cannot be accounted for in the model. The presence of the 
Westward-ctrrent in the, post-dusk sector of the model 
profiles is in good agreement with the observations. 

Finally, the signature cf poleward flowing current in 
the 2300 MLT SPEA profile (the negative Y*-component across 
the profile) is essentially reproduced in the model profile 
Tor 2 500 Mir. 

Another qualitative comparison may be made by 
considering a polar plot of the model results. Figure 


4.18(a) is the polar plot of the model magnetic pertubations 
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j a lot in which the model horizontal 
meapeced 'SS) eae perturbation vectors are rotated 
90° clockwise to represent an equivalent 

current flow. The length of the vector is 

proportional to the horizontal magnetic 


perturbation. 
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rotated 90° clockwise to represent equivalent current flow. 
This figure is to be compared with Figure 4.18(b), taken 
from a recent publication of Friis-Christensen and Wilhjelm 
Meio). In this figure, high latitude average equivalent 
current vectors for ten winters of 1966-1968 are plotted, 


separated according to the value of By, andy By for the, IMe. 


7 
Although there is much more information in Figure 4.18(b) 
than in Figure 4.18(a), it is apparent that the model 
equivalent current vectors are very similar to those shown 
in the bottom three panels of Figure 4.18(b). Since the 
model has been developed from data during periods when B, 
was both positive and negative, it would be expected that 
the model equivalent current vector field would agree with 
observations for which B, was less than zero. 

A final comparison may be made with the data shown in 
Figure 4.18(c), taken from Chen and Rostoker (1974). 
Parthough this’ plot is of data from a relatively disturbed 
period of time, the overall agreement between the model and 


the equivalent Current. vectors of Figure “4. 18{a)) 7 as 


remarkable. 


4.5.2 Quantitative Evaluation of the Model 

As discussed previously, the model has _ not been 
compared to the data by, for example, doing a least-squares 
analysis on the latitude profile curves. Recalling, -that 
AX neak is related to the magnitude of the ionospheric Hall 


Currents, and AY beep is related to the strength of net 
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Figure 4.18(b) 


Polar plot of the ;-average high-latitude 
horizontal magnetic perturbation vectors, 
for winter data, ina format like that of 
Figure 4.18(a). The top row shows data for 
cases when BZ =0, the bottom row shows data 
for By, <-1. The data have also been arranged 
according to the azimuthal component of the 
IMF, as noted in the Figure. (Prom 
Friis-Christensen and Wilhjelm, 1975). 


250 


Lt <> bn Ot Fe ee f AUS Fe 


ren Hea aa as ie Pak ad 


— pli itil puta peepee tnrin ene ore 
ei Ce TV eam a a ee) Wis rs add 
fA ai a ere DM 


om J F ; 
4 (re Fete ai . ere te ha: a 

’ a bad at OLE Oe A RN 
: i} 


f 
j 
Ny 
i 
{ 
ia ae 
1 ee a 
a 


j |  . f Bo i 5 . 

(ey ME aH) i Liat ee et 
he si snl igaatepe t jananate veyed, 
yee ot: aes Re» oon ‘ah "nes ae 
tome . Pigure marie PS Ed yy oy 
i Me Oe Dy Ce ‘potty 
‘ ait tox th eT, ee, ok. alae 
oP ee a ge I, od Y hia an 4 < 
Fick 2 Wigs EME: fy OM NG poked. | « 

ce a een ay Lreemet 3 


pened: ab ORmypettindionie aim PAK 


256 


OVER-HERD IONOSPRERIC CURRENT VECTORS UNIVERSAL TIME 
1 


DAY 290 7 HA 12.2 HIN. 
YERR 1970 


? —-- 20 
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MAGNETIC LOCAL TIME 


Figure 4.18(c). Overhead ionospheric current vectors for Day 
290, 1970. The 22 hour interval shown starts 
at O7 262 UT. (From Chen and 


Rostoker, 1974). 
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Z1eld-aligned current, the magnitudes of AX peak and 
AY step produced by the model current system are compared 
to those observed. Figure 4.19 is a plot of Ot oak and 
AY step obtained from the model as a function of magnetic 
local time compared with the observed values. As described 
in Chapter 3, values near noon are not shown as it was 
difficult to determine a value for AX' in this regime. In 
the near midnight sector, Beeics for AY' were not well 
defined, and so data from this regime is also absent from 
Figure 4.19. The inset shows the values of ea for the 
pre-midnight sector westward electrojet. 

The agreement between the model values of Oe and 
AY! step and the observed values appears reasonably good in 
this figure. However, to demonstrate the agreement in a more 
qualitative way, Figure 4.20 (a,b,c,d, and e€) shows plots of 
the model values as functions of the observed values. For a 
‘perfect fit, the points plotted in these figures would fall 
OR a straight line with unit slope. The straight line drawn 


through the points is the best fit straight line, except in 


the case of Figure 4.20{a). In this plot, the model Se ear 


is compared to the observed AXx' for the § "0200-1000 2st 


peak ’ 
Sector. «it isein this sector that, Large discrepancies arise 
in the north-south component near noon, but it is also the 
sector in which the greatest variability in the X'-component 
was observed. To take advantage of this variability, a 


weighted linear least squares analysis for the best fit 


Straight line was used, in which the weight of each data 
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Figure 4.19 


Comparison of AX'peax and AY'step from the 
model and the superposed epoch analysis. 
Mi——m represents AY ' step as determined 
from the SPEA data, O--N shows AY'step 
for the model. @—@ represents AXfeak aS 
found in the SPEA data, O---O is AXheak 

for the model. The insert is the comparison 
of AX'pyeak for the pre-midnight westward 
electrojet. 
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Figure 4.20(a) 


Figure 4.20 (b) 


Correlation betweenMX'beak for the model 
and AX!" neak for the data for the period 
0200-1000 MLT. The error bars indicate the 
standard deviation. of the SPEEA)” /\ xe 
values. The solid line is the weighted 
regression line, as described in the text. 


Correlation between AY'step for the model 
andAY'step for the data for the period 
0200-1000 MLT. All data have been equally 
weighted. 
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Figure 4.20(c) 


Figure 4.20 (4d) 


Figure 4.20 (e) 


Correlation betweenAX'peak for the model 
and AX'pheak for the data for the period 
1400-2200 MLT All data have been equally 
weighted. 


Correlation between AY Step for the model 
and A\y' ep for the data for the period 
1400-2200 MLT. All data have been equally 


weighted. 


Correlation between AX", ak for the model 
and Nit ey eho amu data for the period 
1800-2200 MLT. “These data are for the 
pre-midnight westward electrojet. All data 
have been equally weighted. 
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point was inversely proportional to the standard deviation 
of the data value. These standard deviations are shown as 
error bars in Figure 4.20(a). For the other comparisons, all 
data points were uniformly weighted. 

Table 4.2 summarizes the results of pie analys?rsi in 
euemacases, the Pslope vor the best fit) line ds unity Within 
the error of the estimate. As well, ali cases have a high 
correlation coefficient. In all but one case, Student's t- 
test indicates that the correlation is significant at the 1% 
level or better. In the case of any for 1400-2200 MLT, 
only 4 data points are available, accounting for the low 
Significance of the correlation. 

On the basis of the statistics, as well as the good 
qualitative agreement of the model profiles with the SPEA 
profiles, it is concluded that the global Current model is a 


reasonable model of the average three-dimensional current 


flow in the high latitude regions. 


4.5.3 Preliminary Summary and Discussion 

in’ Summary “to this point, the global \current Systerm 
consists of 5 ionospheric east~west current systems, each 
with appropriate field-aligned current connecting it to the 
outer magnetosphere, and appropriate north-south current 
systems. The eastward and westward convection electrojets 
grow in the noon sector, being fed by relatively intense 
inward field-aligned current distributed along the 


electrojet regions. In the sunlit auroral Oval, the 
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conductivity arises due to solar UV radiation and particle 
precipitation. At or near dawn and dusk, the UV source of 
conductivity disappears, the electrojet strengths are 
reduced, and current dues upwards to the magnetosphere. The 
westward electrojet continues to grow siege untid<0400GNLE, 
at which time the current diverges both slightly poleward 
and up field lines. This westward current continues as a 
westward current in the region poleward of the eastward 
electrojet until 1800 MLT, at which time it has completely 
diverged up field lines. The eastward electrojet continues 
at constant strength beyond dusk to 2100 MLT, at which time 
it too diverges poleward. When this current reaches the 
eguatorward boundary of the westward current, it flows up 
field lines. 

In the pre-midnight sector, the = tsee due. etotga 
combination of the upward current from the westward 
‘electrojet, the upward current from the eastward electrojet 
and the north-south current. Because of this mixture of 
contributing effects, it is thought that the correlation 
between AY and B, of the IMF in this sector (Chapter 
3) is not as strong as might have been expected. However, in 
the dawn sector, the Seer is due essentially only to the 
Single current system feeding the westward current. As has 
been shown by Hirshberg and Colburn (1969), the strength of 
the westward current is correlated with B, of the IMF, so 


that the good correlation Boe ce Oeerer and B, is not 


unexpected. 
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As well as the two convection electrojets, a very weak 
eastward electrojet between 0400 and 1100 MLT has been put 
into the model, in accordance with the results of Rostoker 
and Hron (1975). Evidence of this current is not visible in 
the SPEA profiles. However, it is not always detectable in 
Single averaged profiles so in all likelihood, its presence 
is averaged out in the SPEA data. The effect of this current 
system is minimal, although a slight positive perturbation 
is visible in the X'-component for the 0900 MLT model 
profile (Figure 4.16) 

Finally, there is the so-called low-latitude current 
system which contributes slightly to the level-shift in the 
Y'-component but whose main purpose empirically is to bias 
that component in accordance with the observations. Although 
certain features of this current system are similar to Sq, 
ateis-more likely that the low-latitude current systen 
‘represents winter ¥ ObEDe=2. 

A qualitative comparison of the profiles produced by 
the’ Nodela»uwithsyntheo’S PEA osprofiles; (-profisseiebyeepuotize, 
indicates reasonable agreement except in the noon sector, 
where the AX profiles show some disagreement. It has been 
Shown however, that by incorporating a small change in the 
geometry of the westward electrojet in this sector, 
qualitative agreement in the X'-component is improved. 
However, even in this sector, the Y'-component shows good 
ee se ith. A statistical comparison between the SPEA and 
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excellent agreement. 

As has been emphasized before, this model, as presented 
to this point, has been designed to reproduce only gross 
average features of the perturbation magnetic field as 
measured at ground level. The pus ston tere eae arises as 
to whether variations of this model may describe specific 


events. The next section gives three such examples. 


To investigate the general applicability of the current 
nodel described in the previous sections, several 
instantaneous latitude profiles have been modelled using 
variations of the average model. It will be demonstrated in 
the following sections that the average model is well suited 
to act as a basis for modelling specific magnetospheric 


events. 


4.6.1 Modelling of Two Unusual Dusk Sector Latitude Profiles 

During a study of energetic particle precipitation into 
the high latitude ionosphere (Rostoker et al, 1978), several 
latitude profiles in the dusk sector attracted attention 
because they deviated narkedly from the average latitude 
profiles in this sector. The deviation was most pronounced 
in the vertical component as may be seen in Figures 4.21(b) 
and 4.22(b). Figure 4.21(b) is the instantaneous latitude 
Hootile for Day 19, 1972, at 0332 UT. This) profiie wonid 


fall into the hourly averaged profile range of 1900-2000 
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Figure 4.21 Latitude profiles for Day 19, to i275 02s 7a 
(a) the result of the model calculation 
(b) the actual latitude profile 
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Figure 4.22 Latitude profiles for pays, boi 27 a0 aia 
(a) the result of the model calculation 
(b) the actual latitude profile 
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NLT. The feature of note in this profile is the behaviour of 
the positive 4~component at low latitudes. It is unusual .in 
that the extremum is at least as large as the AX' extremun, 
and its location is displaced equatorward of Ciewasda? 
location of the AZ peak. In more usual latitude profiles, 
Az peaks at the latitude where AX' is at about half of its 
peak value, and as well, the peak value of the Z-component 
pepavout halt'thetpeak AAX* value. Also in this profile, two 
anomalous inflections occur in the X'-component, at 
datitudes of about 62°N and 67°N. 

Figure 4.22(b) is the instantaneous latitude profile 
for Day 38,1972 at 0214 UT (corresponding to the averaged 
haere tude “profile - range’ of 1800-1900 MLT)S THis profivevis 
unusual because of the low latitude AAZ behaviour as well. 
Although the values of Z-component extrema fall within 
normal ranges, AZ possesses a double peak at the 
equatorward side of the profile. 

The above two profiles have certain characteristics in 
common. Both occur “in ‘the dusk ‘sector, ‘both are Tlow 
amplitude profiles (AAX' < 50nt) indicating weak electrojet 
currents, and both exhibit abnormal behaviour near the 
equatorward edge of the electrojet. The nature*orechys 
behaviour suggested that a current other than a simple 
eastward electrojet was involved in producing these 
profiles. 

To investigate the possibility of an anomalous 


structure in the eastward flowing current, Rostoker -et al 
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considered the electron and proton fluxes precipitating over 
the auroral zone in this sector. Soft particle Spectrometer 
and energetic particle detector data were available from the 
ISIS-2 polar orbiting satellite during 0328-0336 UI on Day 
19, 1972, when the satellite was within 40 = 3° of longitude 
of the meridian line of magnetometers, and during 0209-0217 
UT, Day 38 1972, when the satellite was passing directly 
over the station line. | 

On Day 19, energetic electron fluxes (up to. several 
keV) were observed in the latitudinal EangesoLto9ls toms st.52) 
in close correspondence with the eastward electrojet. These 
Beuxesed dropped® (sharply veat*)69°N,= so sthatt ino electron 
precipitation was present to produce conductivity at lower 
latitudes. There was also insufficient flux of high energy 
peosonss (EP>nAS5S0nkeN)y rtotiaccount >) forrarconductivity high 
enough to produce ionospheric currents equatorward of 69°N. 
sOnee Davyes387hga (.Simitar pattern was observed, with electron 
fluxes (15ekeV SE >&10eeV) confined to) the’ latitude “range 
T4H9N to 68.5°9N. These latitudes correspond well to the 
latitudes at which the AX' perturbation is at half maxinun, 
which in turn correspond to the approximate latitudinal 
limits of the eastward electrojet. Equatorward of this there 
was again insufficient precipitation to account for 
ionospheric currents. 

It will be recalled from section 4.4.1 that the average 


current model contained provision for sunlit hemisphere 


electrojets which exist by virtue of conductivity produced 
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by solar UV radiation. This feature of the model has been 
expanded upon to produce model latitude profiles very 
Shuiubar: “to those of Figure 4.21 (b) and 4.22 (bd) - 
Specifically, the unusual behaviour exhibited by Az in 
these two profiles is believed tobe ite, result of current 
which flows in the sunlit sector due to ionization generated 
by solar UV radiation. This current will be called the UV. 
electrojet. This electrojet is superposed on an auroral 


ee eee we ae ee wen ee ewe ee eo ee 


(eastward) electrojet whose: ~existence "depends "one the 
electric field and the conductivity generated by particle 
precipitation. The UV electrojet extends equatorward of the 
auroral electrojet because there is a significant poleward 
BEfectric field (@.e., <ohe “which will drive eastward Hall 
current) egquatorward of the southern border of the auroral 
oval as defined by energetic particle precipitation, 

This model of current flow is appealing because of two 
ae The A\X' associated with the. UV electrojet may 
have values of the same order of magnitude as AX' for the 
auroral electrojet during quiet times. Thus, the enhancement 
of AX*' seen in the eguatorward latitudes of Fig 4.21(b) and 
4,22(b) may be understood in terms of proximity of the 
observations to the UV electrojet. More importantly, the 
unusual AZ behaviour outlined above can be explained as an 
edge effect of the poleward Pedersen current flow which 
connects the Birkeland current sheets associated with the UV 
electrojet system. 


To this end then, the model as already developed was 
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modified to include a wider UV electrojet than was included 
in the average model. In addition, since the sun sets on the 
F-region about one aour later than it sets on the E-region, 
and the former is the location of significant Pedersen 
currents, the UV electrojet associated north-south systen 
was extended one time zone further into the evening sector 
than the eastward UV electrojet. Figure 4.23 is a schematic 
drawing of the model as used for these examples. 

Tables 4.3(a) and 4.3(b) display a complete summary of 
all the parameter values used. Note that equation 4.2 has 


been recast in the forn: 


Sle es /2, AG) (1+ cosy) 


where Dr and Nes are the latitudes of the borders at 1200 
-and 0000 MLT respectively. 

Figures 4.21 (aj and 4.22(aj) are “the model) Jatvtude 
profiles obtained by calculating the magnetic perturbations 
Oh a meridian one hour towards midnight from the dusk 
fertinator for’ Day 19, while for Day 38, thesdataiwere 
modelled on a meridian 49° away from the dusk terminator 
towards noon, the positions being those of the station line 
at the times that the anomalous profiles were recorded. The 
agreement between the observed and model profile is 
excellent in most respects. The positive AZ recorded at 


Resolute Bay on Day 19 is the main exception. This value was 
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Figure 4.23 


Schematic of the current systems used to 
model the profiles of Figure 4.20{b) and 
Pigure 4.21 (bd). ~Table ""4.3'¢ajy" Pand Table 
4.3(b) describe the model parameters in 
detail. 
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maintained steadily Over several hours and therefore 
Suggests the presence of a relatively Steady state polar cap 
current system which is not accounted for in the model. All 
other differences are relatively minor and could have been 
minimized by further manipulation of the border locations 
and longitudinal and latitudinal variations in height- 
antegrated current density. In view of the fact that no 
model based on a finite number of data is unique, it was 
felt that such "fine-tuning" was not warranted. It is 
contended that the model does account for the observed 
features, and it is concluded that a significant portion of 
the ionospheric current flow in the afternoon sector is 
diverted up the field lines at a conductivity discontuity 


between the dark and sunlit ionospheres. 


Bao. 2 The Substorm of Day 23, 1972. 

Although the current model that has been developed in 
the previous sections is intended to be a generalization of 
the currents flowing during quiet to moderately disturbed 
times, occasionally certain features in data from disturbed 
periods appear to be similar to features in the hourly 
averaged data. An example of such an event is provided by 
data from 0107 UT on Day 23, 1972. Although this day fails 
into the day range from which data has been taken to 
construct the hourly averaged latitude profiles, at wasmor 


used because it was an extremely active day. Figure 4.24 is 


a display of magnetograms from each station, with the 
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Magnetograms for the stations in the 
University of Alberta magnetometer chain, 
for Day 23,-9/2, CO0C stone 200T UT. 


Figure 4.24 
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stations arranged north to south from the toprof thei page, 
for the UT hour 0100 to 0200. Note that near 0107 UT that 
the north-south component becomes increasingly positive at 
jow latitudes, but increasingly hnegative-going at higher 
latitudes (north of SMIT). During the time that this is 
occurring, the east-west component is strongly positive- 
going in the mid-latitudes of the magnetometer line. Figure 
4.25 (a,b,c,d) shows four latitude profiles, beginning at 01 
hours, 7 min, 34 sec UT and spaced 123 sec apart. These show 
Q@m9E0Wth of the X'-componentY indicating, an. increasing 
eastward current flow, and in Figure 4.25(c) and (d), an 
increasingly negative AX' poleward of about 66°N. As well, 
the Z-component is consistent with an increasing eastward 
current and, as the negative AX! develops, the Z-component 
is consistent with an increasing westward flowing current in 
the poleward part of the profiles. The »¥"*-component is 
‘remarkable for its double peak in the first two profiles 
(Figure 4.25(a) and (b)) which merge into a single well 
defined peak. during the later two profiles (Figure 4.25(c) 
and (d)). 

Figure 4. 26e@,b,c,d) show linear mappings™@of all-sky 
camera (ASCA) photcgraphs for 0900 to 0112 UT, taken from 
Fort Smith. These indicate an auroral structure (a westward 
travelling surge) which advances westward towards Fort Smith 
and then remains overhead (Figure 4.26(c) and (d)). While 
overhead, the structure moves slightly poleward. Figure 4.27 


is a north-south magnetogram from College, Alaska. College 
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Figure 4.26 Linear mappings of the auroral forms visible 
on all-sky camera pictures taken at Fort 
Smith on Day i223, Vig72 {slat Gehewetimes 
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Figure 4.27 Magnetogram frem College, Alaska, for Day 
23, 1972. The vertical dashed lines 


delineate the time period of interest. 
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is situated about three time zones west of the magnetometer 
line. The magnetogram indicates an increasingly positive 
north-south component, indicative of a growing eastward 
current. This actual growth of the eastward electrojet is 
emphasized, because, as shown by Kisabeth (1972) it is 
possible to generate a north-south component of the magnetic 
field which appears to indicate both eastward and westward 
current flow by introducing a shear in a purely westward 
flowing current. Such is not the case during the period of 
interest here. 

Up to this point, then, from the AX' and AZ profiles, 
combined with the ASCA data and the observations from 
College, it appears that the profiles of Figure 4.25 are 
produced by an eastward current that is growing in strength, 
combined with a westward flowing current that is advancing 
from east of the magnetometer line to a position overhead. 

The Y*-component is not unlike those seen in the post- 
1800 MLT averaged latitude profiles. That is, the Ax! 
Profiles erare®*consistentsowith “the “existence “of * Locally 
unbalanced upward-flowing field-aligned currents in that 
Ay! shows a negative-going level-shift. However one 
additional feature is present in that the earlier AyY' 
profiles (Figure 4.25(a) and (b)) are double-peaked, and the 
last two (Figure 4.25(c) and (d)) show a pronounced peak 
superimposed upon the level-shift. The proposed explanation 


for this structure is given in the following description of 


the nodel. 
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In order to model these profiles, differential latitude 
profiles were constructed. That is, it was assumed that the 
phenomenon of interest was primarily the structure shown in 
Figure 4.26 (a) through (d), and that this structure was 
Superimposed upon an existing system of currents. Thus, the 
Prot res tor) (Days. 23, 71 hour, 7 min,’ 34 ‘sec wasitakenvas a 
reference level, and was subtracted from the other profiles. 
This has the advantage of removing uncertainties in the 
model of the relatively quiet profile of Figure 4.25(a). The 
resulting three profiles, then, were modelled. 

A schematic diagram of the model is shown in Figure 
4.28 (a,b,c,d). The small circles represent the location of 
the station line, i.e., the meridian along which model 
calculations were made. Figure 4.28(a) represents the 
baseline current system, and consists of a basic eastward 
current (46,125 A) with boundaries at 62°N and 66°N, and 
longitudinal boundaries identical with the average model. 
East of the observation meridian is a westward current, 
which terminates 15° east of the obec ween points. This 
current, as in the average model, is linearly decreasing in 
strength as a function of longitude, (as shown by the shaded 
area), with a maximum integrated current of 2.05x105 A 
flowing at 315° longitude (i.e. 60° east of the observation 
meridian) and between 66°N and 72°N latitude. Superimposed 


on both of these systems is a north-south current system 


consistent with the electric field and conductivity models 


used in the average current model. As well, an intensified 
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Figure 4.28 


Schematic diagram showing the currents used 
to. model the event of Day 23, 1972. (a) 
represents the baseline system, and (b), (c) 
and (d) the development of the currents, in 
time, aie order to model the latitude 
profiles of Figure 4.25. The open-circles 
represent the locations of the magnetometer 
Sites. The shaded area at the leading edge 
of the of the westward electrojet indicates 
the region of intensified equatorward 

current flow. . 
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region of equatorward current flow has been added to the 
western-most 5° of the westward current. The total current 
flowing equatorward in this region is 2.03x105 ADs somthat;, 
in effect, the current at the leading edge of the model 
westward electrojet is primarily eguatorward. This intense 
north-south system is necessary in this model to reproduce 
the double-peaked nature of the AY profile. That is, the 
Y*-component profiles exhibit a double-peak when the 
observations are nade at long distances beyond the 
longitudinal extremes of a north-south current system 
(Kisabeth, 1972). 

Figure 4.28(b) shows the time development of the model 
representing the changes seen in the profile of Figures 
4,.25(b) (1 hour, 9 min, 37 sec). The eastward electrojet has 
increased in strength by a factor of 3.5 to 1.614x105 A. The 
westward electrojet has not changed in maximum strength, but 
‘simply advances 11° westward. Figure 4.28(c) is the current 
system corresponding to the time of 1 hr, 11 min 40 sec 
(Figure’ 4.25¢c)). The eastward ‘current "has increased:an 
strength again, by a factor of 1.6 from the previous figure, 
to »4.214x105 A. The westward current has continued —to 
advance, untilyvaits® Leading’: edge’’ is’ mnowe7S® Pwest. oft the 
observation meridian, so that the observation points now lie 
Slightly to the east of the center of the north-south 
intensification. As well, the westward current increases in 
strength to a maximum value of 8.2x10° A, while the total 


current in the intensified north-south system fremains 
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constant. 

Finally, Figure 4.28 (d) represents the currents giving 
pace.) LO the latitude (profile at 1 hour, 13 min, 43 sec. The 
eastward current again grows in strength to a total current 
of 3.23x10SAmps. No further westward expansion of the 
westward electrojet occurs, but this system expands poleward 
2°, so that the poleward boundary is at 74°N latitude. 
Simultaneously, the maximum current strength increases by a 
factor of 2 to 1.64x10°%Amps. The intensified north-south 
system expands poleward also, but undergoes no increase in 
total current strength. Table 4.4 summarizes the model 
parameters. | 

In this model, the UV electrojet plays no role inasmuch 
as it undergoes no changes in time and therefore does not 
appear in the differential Mvey suber The pre-noon westward 
Current, and the current: flow across the midnight-sector 
have been constructed to be commensurate with the main 
systems described in the preceding paragraphs. 

The model =.latitude ) profites ™ generated @r1ron. > taecce 
current models (and referenced to the baseline system of 
Figure 4.28(a)) are shown in Figure 4.29(a,b,c). (in this 
figure, the symbols H,D refer to Xx*,Y* respectively). It is 
evident that the fit is extremely good. In particular the 
unusual behaviour of the Y'-component is reproduced very 
well. This is achieved due to two features of the model. 
First, the double peak in Ay' is the distant effect of the 


north-south intensification associated with the westward 
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Figure 4.29 


Comparison of the differential profiles with 
the model calculations. The data points are 
indicated by H,D, and 2 (corresponding to 
X', Y', and Z), end. the, model results by che 
continuous curves. 

{a) Comparison of the differential profile 
for. 012092337 UT with the moce!l caiculatveru, 
(bd) Comparison of the differential profile 
for 01:11:40 UT with the model cakculations 
(c). Comparison of the differential. profile 
for 01273243 UL with the model catcugationce 
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current. This intensification is also responsible for the 
positive peak in the sane ccmponent of the later profiles 
(Figure 4.29(b) and (c)). Second, the negative-going level- 
shift in the Y'-component is due to the upward-flowing net 
field-aligned current associated with the west ward 
electrojet. Thus it appears that this disturbance is a 
pertubation of the quiet-time model current systen. 
Basically, the quiet-time model is altered by 
increasing the current in the auroral electrojet and both 
increasing the current in the westward electrojet and 
letting it expand westward from the nominal position 
assigned to the westward electrojet in the quiet-time model. 
The westward current probably expands due to an enhanced 
particle precipitation into the region poleward of the 
eastward electrojet thus providing a Hall conductivity where 
little or none existed before. The intensification of the 
intense equatorward flowing current at the leading edge of 
the westward current is probably due to an enhanced Pedersen 
conductivity in this region, iee., an increase in relatively 
soft particles, although an enhancement of the equatorward 
electric field in this regime would also bring about a 
north-south current enhancement. Such a change in the 
electric field would also cause an enhancement of the 
westward flowing Hall current in this same sector, although 
the magnetic field has been successfully modelled without 
invoking such a localized enhancement in the westward 


Current. However, in view of the lack of corroborating 
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particle and electric field data, it is not possible to 
apply a unique physical interpretation to this disturbance 
model. One can only interpret: the results “in terms. ef 
-equivalent current flow. This model has been presented here 
only to show that, indeed, a variation in the model for 
average conditions is remarkably successful in modelling a 


variation in real substorm current systems. 


4.7 Comparison of the Model for hverage Conditions with 
Satellite Observations of Field-Aligned Currents 
Reference has been made in preceding sections of this 

chapter to the fact that the glcbal current model of section 

4.5 has a field-aligned current distribution that does not 

compare well in all aspects with the interpretation of some 

polar-orbiting satellite observations, specifically, data 
from, sthe Triad satellite (Sugiura and Potemra,, 19/0; Tijima 
aiieesroremra,: 1976). Statistical .charactertst1ce ot.) tie 

Piecld=aigned “current. flow as ~wuntierrnedg@ trom | them iriad 

satellite data have been shown<}in@Fagure 4.2.6 JAs “well, 

Figure 4.30 shows the most recently determined average 

characteristics of the upward and downward field-aligned 

current densities as a function of latitude. Combination of 
the data in these two figures indicates that although the 
satellite observations imply a strong unbalanced downward 
current flow in the pre-noon hours, there is a similarly 
strong net upward flow in the immediate post-noon hours. [It 


is in this latter region that the disagreement between the 
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satellite observations and ground based observations arise. 
Elsewhere, the two data suites are in good general 
agreement although the ground based data indicates that a 
pair of east-west aligned field-aligned current sheets may 
extend as far into the dusk sector as 1800 MLT, whereas the 
Satellite data shows such a structure only until about 2200 
pets this discrepancy 2s not fully understood. “Tijima ‘and 
Potemra (1977) present a figure similar to Figure 4.2 for 
more disturbed periods and show that even then, the field- 
aligned current pair that would be associated with westward 
current flow only extend up to 2200 MLT as well, although 
Kamide and Akasofu (1976) present Triad data in which the 
Substorm westward electrojet and the associated field- 
aligned currents extend to 1900" MLT.." As “well, ‘the 
latitudinal extent of this westward electrojet is small and 
the associated field-aligned currents occupy a narrow 
scegqion, So that, “On a statistical” basis, these currents “are 


not evident in average pictures of field-aligned current 


flow developed from Triad data. 


4.7.1 An Hypothesis to Explain the Post-Noon Discrepancy 

As Hae cue eea Garlier, a polar=-cap ‘CurccrenC)  Syotem wihac 
not been included in the model developed for this thesis 
because it is not believed that there is, on average, 
sufficient conductivity in the polar cap tO support a 
significant current flow during the winter months. However, 


almost all the Triad data that has been published has been 
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from the summer months when one might expect ionospheric 
polar cap current systems to exist. The average features of 
polar cap and auroral zone magnetic perturbations in the 
“Summer are markedly different from those in the winter. For 
example, when one compares Figure 4.31 showing the average 
Magnetic perturbation vectors {rotated into equivalent 
current directions) for cae Summer MONthS “Of 1966) through 
1968  (Friis-Christensen and -Wilhjelm, 1975), with. Figure 
#.17(b). it is evident that the high latitude equivalent 
current vectors are rotated sunward in the summer, relative 
to their direction in the winter months. 

As well as much of the published Triad data being from 
Summer months, these published data are also found to be 
confined to periods during which By of the IMF is less than 
PeLOc this |1S =" not. to say that ali the Trvad) daca thatihas 
been analyzed has been from times when By <Q. Indeed, the 
data have not been selected on this basis 
(Potemra, 1977(b)). This observation is pointed out because 
when BycOe in the summer, the Svalgaard effect is indicative 
of a counter-clockwise polar cap current on the dayside,. 
That ics, there is a, pre-nocon to post=noon current @ilow 
across the polar cap. 

A possible explanation for the difference between the 
Triad summer data and the ground-based winter data can _ he 
developed based on the existence of a cross polar cap 
current system, particularly when By <0.) SEIQUTC (eG. 32ers 


reproduced from the paper of Iversen and Madsen (1977). This 
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Figure 4.31 Polar plot of average high-latitude 
horizontal magnetic perturbation vectors, 
Tor summer data. (From Friis-Christensen and 
W2thjein, 975) 
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Figure 4.32 Average ionospheric electric field vectors 
for times when B, < 0. (After Iversen and 
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is a plot of the average electric field for By <0. The noon- 
sector, high latitude average electric field is consistent 
with a cross-polar cap current. Recent work by McDiarmid et 
al GIILS), 2usdng. data eafrom «cthe: pebSTS<2 polar-orbiting 
satellite, shows that, for By <O,0n there» dsj, a Pinegions cof 
downward current immediately after local magnetic noon and 
well poleward of the field-aligned current sheets observed 
by Triad (e.g. see Figure 4.2). (This current has also been 
observed by Iijima and Potemra (1976), but the relationship 
to By was not elucidated). This is also Ape Peat with a 
cross-polar cap current, of somewhat limited extent, but the 
Tight sense, if this field-aligned current connects to the 
upward current sheet at the poleward edge of the auroral 
Guat tin ithe ipost-noon jusector.,») Thus, «there «sis! sufficient 
evidence to justify examining the effect of a morning-to- 
afternoon cross polar cap current system. 

Thee cross-polarne cap-, current, systems chat,| wid be 
developed is very similar to the Sy current model of 
Kawasaki and Akasofu (1973). These authors used a flat-earth 
approximation, but allowed field-aligned current to flow 
along dipole field lines. In the present model, a spherical 
earth has been used, and also field-aligned current has been 
allowed to flow along dipole field lines. Cross polar cap 
flow is modelled by current flowing along great circles. 
(see Appendix I). This current is connected to field-aligned 
current sheets which are placed symmetrically about noon, 


and along lines of constant latitude. A uniformly 
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distributed current flows into the ionosphere on the morning 
Side of the polar cap, and out of the ionosphere on the 
afternoon side (see sFigure sAl.2). “As well, a strictly 
ionospheric current which circulates around the field- 
aligned current sheets has been included. The current 
densities have been adjusted to give a Hall to Pedersen 
conductivity ratio of 2, consistent with the value used 
throughout this thesis. It is emphasized at this point that 
this polar cap model has not been developed in an attempt to 
nodel real data, but simply to provide a possible 
explanation of the Triad/ground-based data discrepancy. 
indeed, it is believed that this discrepancy is probably a 
Summertime feature, and ground-based summer data were not 
available in sufficient quantity at the time of this study 
to permit a detailed analysis... 

For the purpose of this discussion, the polar-cap has 
been defined as the region poleward of 75°N latitude. The 
field-aligned current sheets were positioned on the 75°N 
latitude circle, over longitude ranges from 0700-1100 MLT 
for downward current, and 1300-1700 MLT for upward current. 
A total of 5 x 10% A flowed across the polar cap connecting 
the field-aligned current sheets (Figure NIDAGA yes 

The computed perturbation magnetic field for this 
current system is shown in Figure 1.33% -LaLS SYSteminas oe 
been added to the quiet-time model since the specification 


of the borders in the two systems is not identical. However, 


it is evident that a current system like that shown in 
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Figure 4.33 


Polar plot of the hagnetic field 
perturbation vectors due to the polar cap 
current system described in the text, 
rotated into the equivalent current 
direction. 
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Figure 4.33, if added to a system like that of the quiet 
time model (Figure 4.18 (a)), would rotate the equivalent 
current vectors of Figure 4.18 (a) into a pattern very 
Similar to that observed, for summer, by Friis-Christensen 
and Wilhjelm (1975) (Figure 4.31). As well, and this is of 
greatest importance, the polar-cap current system produces 
latitude profiles in the post-noon sector which could be 
interpreted as signature of net downward field-aligned 
current. The Y'-component shows a positive-going level-shift 
(Figure 4.34 (a,b,c,4)) araten ds occurs where it is known, 
from the model, that the field-aligned current is actually 
directed upward. This result occurs primarily because the 
Magnetic perturbations due to the ionospheric current flow 
are larger than the perturbations due to field-aligned 
current. 

It is suggested, then, that during the summer, when 
By<0, current flows across the polar cap from the morning 
Bector to the afternoon (Sector, and that this current 2s 
most intense when By <0. This current is connected fo ¥Ehe 
magnetosphere via field lines which lie along the poleward 
border of the auroral oval. During times of moderate 
activity, the field-aligned current associated with the 
cross polar cap current could be sufficiently strong to 
cause the satellite to measure a magnetic signature 


consistent with a net upward current in the post-noon 


quadrant. The existence of this current however does not 


preclude the existence of a downward field-aligned current 
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CHAPTER S -:" POSSIBLE SOURCE MECHANISMS 


Up to this point in this thesis, nothing .- specific has 
been Said regarding the mapping of the field-aligned 
currents to the outer magnetosphere. Indeed, this is still a 
topic of active ieee aiern Thastesectizon —waltie bree y 
Outline some of the thoughts from the literature concerning 
sources of field-aligned currents. 

Sato (1974) has given a theoretical overview of 
possible field-aligned current sources. Assuming frozen 
mreid Conditions, ises, E +V¥ XB =10, “whererkly, -andee Beecare 


the electric field, velocity field, and magnetic induction 


field vectors respectively, he obtains, 


Voor ye(y so) 7 O-(75i) 5.1 


Mes tirstetermsofithis: isreoften “zero,” butt’ hase non-zero 
Values if a’ shear. «@€xistS in the.magnetic field component 
normal to the plasma flow. Such can arise .on a» boundary 
separating regions of relative motion of open and closed 
field lines. Thus, V-(V~x8)#0O if a current flows in the 
Same direction as the plasma flow. 

The second term of equation 5.1 is thought, in general, 


to be the most important since there is a shear in the 


velocity component normal to the magnetic field at the 


boundary between open and closed field lines. 
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Thus, VEE =O on this boundary, implying a charge 


accumulation on this boundary. Sato states a theorem 


concerning field-aligned current, based on equation 5.1: 


"If the plasma convects around a point in the same 
sense aS the proton gyration, positive charges 
accumulate there, so that a field-aligned current 
BLOWS Out from that. points On the “other ew elekn yale 

ape) plasma, convection is in the same sense asthe 

electron gyration HOtLOn, negative 

charges (electrons) accumulate, a fielid-aligned 
current (thereby flowing into ‘the “center ot) the 

WOrvex. in other words, if the vorticity vector 26 

parallel to the magnetic vector, a field-aligned 

CUbLent PLOWS ln, DUE AL It) 1S: antieperaiiel,— sa 

Eeeid-aligned curpent flows out." (Sato, 11974). 

For the actuai case of the magnetosphere, in the region 
of open field lines above the polar cap where the field 
lines are directed downward, the above theorem predicts a 
downward field-aligned current on the dawn side of the polar 
Gap, and an upward ficld-aligned current on the, dusk side. 
This dynamo or magnetohydrodynamic (MHD) generator has been 
discussed by Akasoful (19747 1975,1977) as a possible 
mechanism for driving substorms. Figure 5.1 is a schematic 
of this, and in particular, panel (b) shows the operation of 
the dynamo. Solar wind plasma is driven in the ty direction, 
across the magnetic field (B), oriented in the +z direction. 
This results in the generation of the -¥ x B electric field 
in the -x direction. Current flows down the field lines on 


the left of the figure, and the entire system constitutes a 


generator. Akasofu (1975) has also pointed out how such a 


system can explain, at least in part, why auroral activity 
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Figure wo. } 


A schematic diagram indicating the processes 
associated with the solar 
wind - magnetospheric dynamo. Panel (a) 
Shows the location of the dynamo, {b) 
indicates the basic processes, and (c) shows 
the connecting Circuits, (Aiter 
Akasofu, 1977). 
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occurs in a particular region (the auroral oval) around the 
earth. For example, when the IMF is directed southward, then 
a “neutral line" surrounds the magnetosphere in the 
equatorial plane (Figure 5.2). The morning. half OL ties bane 
actS aS a positive terminal of the dynamo,the afternoon half 
Somecire, Negative termingl., Thus .cirrent “flowe (imto tre 
fonosphere ini the morning half, and out of the ionosphere in 
the evening half, as pointed out above. 

In the magnetosphere, the plasma sheet is filled with 
hot plasmas which may give rise to intense drift currents. 


The drift current can be described by 


aes ml ek +7 


where 
[A 
oie -V x Zr 8 
he AN V8 
vee oP 


pol 


aA 
vod being the charge density, b (= %) the unit vector in the 
direction of the magnetic field, and in and ¢. the 
components of the pressure tensor perpendicular and parallel 


to & respectively. Taking the divergence of equation 5.2 
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Figure 5.2 A schematic diagram showing how the 
"terminals" of the solar wind dynamo connect 
through the ionosphere. (After 
Akasora, 9197 5)). 


322 


my 
2 


4) RL 
3 Cae 
noe 
‘ t 
$y Aezy 
) oh wee ol me i a wp 7 
: i ; , 
YY 
Fi ' i 
Aa ' - 
\ 
LS 
al 
Pee ae: OP eee ' 
‘ is 
de M 
ri 4d Luby Cale bles 
oo ¢ 
E i : 
te im 
i fal 
‘i 2 le 
' i ¢ ms 


ae a PY ah la ‘| is 


csepcnig tea on ron owen | 
1 fi ry v i - 7 
i : 


i | 7 


) i His ." e 


(i) % pia hii ial ape Steteteteas 


‘low 


Jee 


leads to 

Se V+ Lay i Vas + O-Jool 5.3 
AL fet 7 iiohar re wimg Y dE 

ee Cg cat V8) -(b ~bV)b (V8 «Vale 


The polarization term in general opposes charge accumulation 
due to the other terms. Thus, a current divergence in the 
magnetosphere may arise if the plasma pressure has a 
gradient in the direction of either the curvature current 
anadyor the V7 Bucurrents rien interesting to note in passing 
that the vortex-like generator has the nature of a voltage 
generator, whereas the source due to curvature and gradients 
in Bis a current generator. 

Finally, aif, thejedonosphere phasjggradientcniny placna 
dencity, it no,dongervactsS. as de» passives closure! epathpator 
field-aligned,|,currents, »but.' behaves, as .ay,,field=aligned 
current source due to the development of polarization 
electric fields. Height-integrated ionospheric current 


density is given by 


Ne ee a ae 5.4 


~~ 


where oa and 24 are the height-integrated Pedersen and 


Hall conductivities respectively. If it is assumed that E is 


irrotationads, «then 
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In the case of an homogeneous ionosphere, field-aligned 
currents arise only because of gradients ERE, andeEehas 2s 
origin in the magnetosphere. In this case then, the source 
of Jy is the magnetosphere, and the field-aligned currents 
are connected to the ionospheric Pedersen currents. In the 


case of an inhomogeneous ionosphere however, gradients in 


z 


p and Pays Wil Avead Ato Jy also. 


As described earlier, the solar wind dynamo produces 
field-aligned current sheets ‘around the polar cap. These 
Sheets correspond to the poleward side of the north-south 
current system in the model of Chapter 4. Sato argues that 
the current sheets on the equatorward side of the auroral 
oval arise from ionospheric conductivity gradients, and that 
therefore, the ionosphere regulates magnetospheric 
convection. The gradients in ionospheric conductivity govern 
the field-aligned current closure with the result that 
polarization. ‘electrics > fields*® @arek set umupleihesesoLectric 
fields map to the outer magnetosphere and drive convection. 

Cole? (1961), 1974,1976) “has discussed ae =nodely pinury which 
field-aligned currents are connected to a dynamo which is 
located in’ a boundary iayer of plasma inside the 
Magnetopause. He has shown that solar wind protons may 
penetrate the magnetopause on the morning side, and _ solar 
wind electrons may do so on the evening side, thus 
generating charged boundary layers which may act as a source 
of ppchaeaeGuse currentse He also points out that the 


dynamo action of ionospheric winds may generate 
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electrostatic fields or gradients of LOUTZACL ON -andwrenat 
these fields may map to the hagnetopause (Cole, 1976). 
Eastman et al(1976) have also discussed a boundary layer MHD 
generator as the source of the field-aligned currents which 
bound the polar cap. 

Potenra (1977) has renarked that; due Comretie 
Statistical stability of the poleward field-aligned current 
Sheets, as inferred from the Triad sateilite magnetometer 
data, these field-aligned currents constitute the primary or 
driven current system, and are associated with boundaries of 
the plasma sheet far distant from the earth where buik 
plasma convection may act as a generator. Further, he 
Suggested that the equatorward field-aligned current sheets 
are secondary currents, which map to the inner edge of the 
plasma sheet. These currents exist in response to localized 
variations in, for example, ionospheric conductivity. 

An alternative driving mechanism has been discussed by 
Rostoker and Bostrom (1976). They confined their attention 
to field-aligned current flow in the dark hemisphere only. 
Based on ene work “of? Prank (1971)47 “Dassene(1970))7-e7ancd 
Bostokem= vet) al, 7(1975), they assumed thatethe plasia sreet 
maps into the auroral oval, and that the poleward edge of 
the oval maps to the boundary of the tail lobe and plasma 
sheet (see Figure 1.1). Thus, field-aligned currents map 
into the plasma sheet, and flow on closed field lines. 
Rostoker and Bostrom assume that static forces in the tail 


are unbalanced, and that these forces lead to an outflow of 
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plasma from the center of the plasma sheet to its flanks. 
This outflow leads to a potential difference between high 
and low latitude field lines, which in turn drives field- 
aligned currents in pairs. Figure 5.3 shows the resulting 
field-aligned current flow. 

It is clear from the above discussion that one would 
expect downward field-aligned current on the morning side of 
the polar cap and upward field-aligned current on the 
afternoon side. it is not, however, patently clear how the 
field-aligned current pattern determined in this thesis fits 
into the above source mechanisms. Certainly, in the model, 
the north-south current system with its associated balanced 
field-aligned current flow satisfies, in general, the 
pattern of downward current flow on the morning side of the 
polar cap, and the reverse on the afternoon side. However, 
the existence of a relatively intense inward current flow in 
the noon sector, and a similarly intense outward flow in the 
pre-midnight sector is not explained in terms of the above 
models. Figure 5.4 is a schematic in which only certain 
features of the model current system are shown. The eastward 
and west ward electrojets are indicated as Hall 
currents (I, ), connected to downward field-aligned current 
near noon, and upward field-aligned current near midnight. 
Across the midnight sector, the westward current is 
represented as a separate current systen, with the 
ionospheric part being a Pedersen current (I, ), and 


connected to separate field-aligned currents. Combined, the 
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Exrqgureito.. 3 


A schematic of the fieid-aligned current 
flow generated by the Rostoker-Bostron 
dynamo. The view is from down-tail towards 
the earth. The field-aligned currents are 
indicated by and and are confined 
to the plasma_ sheet. (After Rostoker and 
Bostrom, 1976). 
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A simplified schematic of the global current 
model. The main eastward and westward electrojets 
are shown as Hall currents(I,,). The westward 
electrojet in the Ai ee sector is shown asa 
Pedersen current (I r consistent with the west- 
ward electric pent in this. sector. The Hall 
currents are fed by downward field-aligned 
current ((@).. concentrated primarily in the noon 
sector. These Hall currents return to the 
magnetosphere as upward field-aligned current (@) 
in the midnight sector. The westward current in 
the midnight sector (I) is fed by a downward 
current to the east equal in strength to the 
upward field aligned current shown in the same 
sector, so that there is no net field-aligned 
current in the immediate post-midnight sector.In 
the pre-midnight sector, the upward flowing 
current connected to the eastward Hall current 
and the westward Pedersen current add to give 

the observed net upward field-aligned current. 
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three current systems represent the major currents of the 
model, less the north-south current system. 

This breakdown into a Hall and Pedersen component is 
somewhat appealing, at least in part, because it permits one 
to account for the known potental jade: acne the polar cap 
Bio to; d0 KV (AxTond ef a1,;1965)}*matithe model, there: is 
HOW POlar cap —curnbent- system for ithe wamter months so that 
the current flowing into the ionosphere: in the pre-noon 
hours, and around the oval into the revue don eeeseuor 
cannot be purely Hall current. However, in the midnight 
sector, the electric field has a westward component, as does 
the current. This, ai the midnight sector, the westward 
electrojet has the nature of a Pedersen current, and is 
therefore capable of dissipating power. Indeed,if tae value 
of 10 mV m-1 is taken aS a representative value of the 
westward component of the electric field, and it is assumed 
that this Pedersen current flows across four time zones at 
an average latitude of 70°N, then the potential drop is 
amproximately 38 KV. THis requites @thaty thegifreid=avigned 
current be connected to a magnetospheric generator, and the 
mechanism such as that described by Akasofu (1977) and as 
outlined earlier in this chapter, may be invoked. 

To explain the pattern of field-aligned current flow 
connected with the Hall current electrojets, the following 
model is proposed. Figure 5.5 is a schematic drawing of a 
possible plasma convection pattern in the equatorial’ plane, 


or equivalently, an electric equipotential pattern. This 
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Figure 5.5 


A schematic of the magnetospheric convection 
flow pattern required to provide the 
fieid-aligned currents that are connected to 
the 1onospheric Hall currents. The +" 
indicates the region tron where 
field-aligned current will flow into the 
ionosphere; the "-" indicates the region to 
which field-aligned current will flow from 
the ionosphere. 
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pattern is not unlike one given by Harel and Wolf (1976), 
and is reasonably consistent with the average high-latitude 
electric field pattern {see Figure Se LSD). Pilaswha, als 
convected earthward from the tail and diverges to flow 
around the earth. In the evening sector, most of the plasma 
is turned toward the nearest flank long before it reaches 
the dayside regions of the magnetopause. However, some of 
the plasma flowing around the dusk side penetrates into the 
noon sector where it is deflected tO flow towards “dawn. As 
well, all the plasma which flows around the dawn side of the 
@€arth returns to the tail along the dawn flank of the 
magnetosphere. In the regions where V is not curl-free, 
equation 5.1 predicts a divergence of the electric field 
such that a positive charge would build up in the region 
labelled + in Figure 5.5, and a negative charge would build 
up in the region labelled -. These charge concentrations are 
dissipated via field-aligned current into the ionosphere 
fLOueethe + -region, and “Out "Of “the 1 onespheresinatne, 
region. 

An equivalent way of describing this has been given by 
Vasyliunas (1972(a), 1972 (b)). He shows that azimuthal 
pressure gradients result in field-aligned current flow. In 
the regions where the plasma is being diverted back down 
tail, there is an enhanced particle pressure which result in 
field-aligned currents which flow in the directions 
described above. Vasyliunas (1972(b)) also Shows that these 


pressure gradient field-aligned currents may arise from the 
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interaction of the ring current (see below) with an electric 
field, and that these field aligned currents can be related 
to an ionospheric Hall current. 

In section 1.3.2 of this thesis, the presence of a 
westward flowing CML a= Ceri eGt has een lng current was 
mentioned in connection with Dst. Frank (1967) arse 
detected this current aireceay from an analysis of particle 
data from the Ogo 3 satellite. As well-as -this.-Symmetric- 
furEent,. there is an asynmetric ring currént (Cahill, 1966; 
Pian, 19/0), and | several authors (Fejer, 1961; 
Swift, 1967, 1968; Akasofu and Meng, .1969) have suggested 
that both the auroral electrojets are connected to the 
partial ring current. It has also been proposed that tae 
westward electrojet connects to a partial ring current which 
lies at a greater distance from the earth than the one to 
which the eastward electrojet connects (Crooker and 
Merpherron, 1972) 5 

Thus, it is suggested an his thesis) sthar he Preloid— 
aligned curnents connected to the Halil) currents (Figure 33.5) 
arise from within the magnetosphere in the manner outlined 
above. further, it is suggested that these ticid=alignued 
currents are connected to the partial ring current to forma 
closed loop. More details of the actual convective motions 
in the dayside of the magnetosphere than are currently 
available are required to test the consistency of the above 
hypothesis. In the absence of this, the above description 


provides a mechanism whereby the convection electrojets are 
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connected to the magnetosfhere through fielid-aligned 


currents. 


See ee ee ee ee ee 


A detailed study of ground-based mcneeomenen data has 
been carried out using the superposed epoch analysis 
technique of Chree and linear inversion techniques. Based on 
this study, a comprehensive three-dimensional model of 
ionospheric-magnetospheric current flow has been developed. 
The data and the results of the model calculations have been 
parameterized on the basis of the peak value in the Xt- 
component of the magnetic field, and the level shift in 
QY', and these parameters have been compared statistically. 
Within the margin of error in the data, the fit of the nodel 
to the data is excellent. 

Although many features of the model are not new, 
certain new results have come to light. The presence of both 
eastward and westward electrojet components dependent on 
conductivity caused by solar UV radiation has been detected. 
The existence of westward flowing current poleward of the 
e€astward convection electrojet in the pre-midnight sector 
has been demonstrated by using linear inverse techniques, 
and, by forward modelling, it is shown that this westward 
current is associated with upward-flowing, locally 
unbalanced field-aligned current. 

Finally, it has been shown that both the eastward and 


westward electrojets are connected, in the noon sector, to 
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the outer magnetosphere through locally unbalanced field- 
abigned currents) © It iis Suggested that these field-aligned 
currents arise from asyumetric convection of plasma on the 
dayside of the magnetosphere. It is further suggested that 
the ionospheric Hall currents are connected through these 
field lines to the ring current, as described by, for 
example, Vasyliunas {1972(b)). 

It is believed that the current model developed in this 
thesis is important for several reasons. First, it has been 
tested against a real data suite, and shown to be consistent 
with that data Oh average. Further, the model has 
demonstrated that, for low-level magnetic activity, the UV 
electrojet can be important to the complete description of 
ionospheric-magnetospheric current systems, Also, when 
perturbed appropriately, the model can be used to describe 
- substorm current systems. Thus, a generalized current model 
has been developed which may be used as a base system for 
modeliing real current fiow for many levels of auroral 
activity. Finally, field-aligned current in the noon sector 
to feed the eastward and westward electrojets has never 
before been included in an ionospheric-magnetospheric 
current model. The success of the model in reproducing the 
nature of the AyY' level-shift indicates that these currents 
are required. 

The model described in this thesis is not presented as 
the ultimate in ionospheric-magnetospheric current models, 


as all current models are subject to revision until the 
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2 


entire underlying physicssas understood. It “is; hoped that 
the work described in this thesis will serve the role of 


providing one more step in this search. 
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eee OLX & MODELLING OF THE CURRENT SYSTEMS 


AT 1 he magnetic field due to a three-dimensional 
Kisabeth(1972) has developed equations which facilitate 
calculation of a three-dimensional current system. This 
appendix summarizes the develcpment of these equations. 
Figure AI.1 defines the vectors used in spherical 
coordinates for the calculations..“TherBiot-Savart Law is 


formulated as: 


g = fel /f T=) d’r ee 


ioe = if 
where 
fo is the permittivity of free-space 
J is the veurrent deasaty i ngaticc 
we =(%, ~)= observer coordinate vector 
and t= (rf 8G) = sounce coordinate vector. 
The coordinate directions at the source (7.0,—) are 


Z : Ae 
different than those at the observation point (iC eer 


but they are related by: 


A 
A r 
£ x 
A 8 


AL.2 


pine AY 


where A\ is an orthogonal matrix with components: 
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Figure AI.1 Diagram defining thers vectors ~ inn spherical 
coordinates used in the magnetic field 
calculations. 
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A, = sin8, sin Ocos(-g) + cosO, cos O 
A = sin®, cos 8 cos (QZ -y)— sin cos, 
Ais = stxO, sin(g? -Y) 

= cos 8, sin 8cos(p -() - sin 8, cos 

Ay = cos®, cos Ocos (yp ~~) + $inB, sin 


> 
e 


Axx = cos®, sin(P -f) 
A, = Sin Osinlp —p) 
Ay, = ~cosOsin(g -¢) 
As; = cos(p-y) 


At a point (a, Of Q, ), the magnetic field froma line 


current can be expressed as 


gine res dCy ar 


where 5; aT k ig) COL as} takes on valuestof Vij 2A ore: 
respectively; and 2 is the curve defining the current 


path. The matrix ae has components: 


O “A3, Ly als, a Ax a,ds, 
FG = “ieah ¢, ds, ~ (rA,48, + 0,A4))d5, (7A 332, Azz 2 ads, 
i. Maes ds, (Ae, -7 Ai 2, )ds, (Ak, + rh Tore 
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radius of the earth 


Qt 
6 = radius of the superconducting sphere 
asa 
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{Rf ae 7r* uP T tC Ea cosxX (cosX =A,,) 
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6 ey fia BLP CGN 
ie TR’ R+r-Z‘cosxX 


R's r’—rcos 


T, (cosr <0) = rTcosh ~ 2 «nf Sores 


a r (1- cos) 
T, (0s > 0) = r[,cos R tIn SET ES -7, nO 


and ds = (ds ds,ds,) depends on the current path. For 


example, for current flowing along a path given by 


O= 6+ %(0,-8)(L+c5¢) AL 4 


i, oy) Ge) i 
; ee en. m, { 
heh e h : pai i 
ire 
ed eee = Le 2 Mabe aate nit VENChars 
7 oe i, os ee oe pale et me oy 
X | MC pray." 
teh an A A 1 u 
Py : \ 
= A b : a y 5 = re 
pat gta (UM Thee Ae ie a eee es.) “) fark ehicty. iow i 
wet 
i 
ot 
. 1 
ei 
a 
rat 
} 
=<! 2 PF, ca 4 . pean 4 


tan : faa prey 
fos a 


1 tense te . i 
Mecoes 8 tiie 2 Bm hei 10h i 


j i 


i o 
Ce ae ee ae Howie ee aa 
vd lilies ‘dug & lis Rae ‘Stesres ‘rot tec 


bs (edt) (S- ORs SB 8) 


5 he 


which describes the locus of the auroral oval as used in 
Chapter 4, ds is found as follows. Rewriting equation AI.4 


as 


6=c+deay 


where 
C =. pee) 
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ACL), 


and defining [| = @-c -doos¢ 
Boek, a unit vector perpendicular to the | current ‘path a1s 


given by 
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The unit vector parallel to ds ’ uy , is given by 


A 


L, 


I 


A A 
Uae 

L A ; A 
_dsingp@ + Sin8 


= je 
(sin*O +d‘sin*¢) 
The differential path length,a@s, is given by 
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ds =i,ds 
r(-dsinpé +sindp)dyp 


and 
ds =o 
aes Sen gs 
ay = rd sin pty) aa i (8,-8)sinepde 
ds, = rsinbde 
For current flowing along a field line, 
ds = Zrcotede 
ds, = rdo 
ds, = O 
Ar 2 A Cross-Polar Cap Current system 


In Chapter 4, the magnetic field due to current flowing 
across the polar cap is discussed, Tolcaiculate this fields 
the method outlined in the first section of this appendix is 
used, but first an expression for the current path must be 
formulated. It is assumed that the ionospheric part of the 
current is\‘confined to great”° circles across the) polar cap, 
and further, that these are symmetric about the noon 
meridian and perpendicular to it (Figure AlI.2). 


The equation of a geodesic is found by minimizing the 


arc length, i.e., minimize 
4 
Li Jas AI.5 
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where @$ is the differential arc length, and the integrai is 


carried out from points 1 to 2 on the curve. 
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FOr La sphere, ds=r(1+ $*sin*o) “do , where p= ty, 


Then, 
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Taking the variation §Z ana setting it equal to zero leads 


to the result (Goldstein, 1950) 
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ae Oe 4 
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or, after some manipulation, 
~ — c esc“ 

2ZIN7 
sin 6 (1-¢*—c*cot’9)* 


where Colts a-constant. 


Integration of equation AI.7 (see 


Murphy, 1956) gives 
p= X—sin (Reot 


which is the equation of a great circle. 


For the case shown in Figure AI.2 , 
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great circle are Qf) and (& -Y ), or (& 4%) and 
(O,-f). Substitution of these into AI.8 permits the 


definition of ™% and Rae with the result that 


col 8 = £0f8 cosy 


RA) 


Equations AI.9, together with 0 =0, , ‘define the 
boundaries between which the current flows across the polar 


cap, as discussed in Chapter 4. 
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APPENDIX TE LINEAR INVERSE THEORY 

One of the most important methods of analysis developed 
in recent years to aid in the definition of parameter 
Variations in regions not assessable by direct devituneneave 
is the technigue of linear inversion. While a complete 
gescription of this technique is~ not possible in’a few 
pages, a brief description of the principles of linear 
inversion theory will be cutlined here. For a more detailed 
presentation of the technique, the reader is referred to the 
pioneering papers of aAonus ander Gri berks t(19677-19 (0) "ee ech 
mecent srapplhications) . of | Oldenburg 7976)" + and=4a general, 
review by Parker (1977). 

The three-dimesional east-west current systen 
feabeth, 1972) (see Appendix I of this thesis also) is used 
as the forward model to obtain height-integrated ionospheric 
current densities. Once the parameters defining the geometry 
of, + ithe. current: system (longitude®of the caster and western 
field-aligned current sheets, and the northern and southern 
latitudinal boundaries of the ionospheric current) have been 
specified, the magnetic field observations at the surface of 


the earth {r=R), colatitude oP e and longitude f, are 


obtained from 


8, 
ap ay * de ee — 
G (2,22) -| £,,(0)G, (R,0,,2; 8) i L235 AII.1 
4 
ae He t Cape ea 
where B; is the f measurement of any componen or e 
magnetic field due to the east-west current systen, J-,(@) 
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2 


is the height-integrated current density, in Am-1, and 8 y 

and 6, are the colatitudinal limits of the ionospheric 
Ew - 

current. G; 1S the Frechet kernel or Green's function for 


the problem, and .is' a ‘function “of the (current, systen 


geometry and the position of the observer. In particular, 
et tox cle 
pols Y AII.2 
Te 


where hs and dG are defined in Appendix I. 
The magnetic contribution for a north-south current 


system may be calculated in a Similar manner, viz, 


a’ 8 2 NS 
F822) -[ TOG RO,g;0)d0  j-423 nee 
| 6, 
as has been shown by Kisabeth (1978). 
In any realistic modelling of current systems, both 
| north-south and east-west current systems exist 


—_ 


simultaneously in any region, so that the datum 6; is given 


by 
EW NS 
8: = Gye OP Gee 
i) J J AII.4 
The height-integrated ionospheric current is given by 
Jus in > > Ens ~ Ee 
ALI. > 
pe ie Pais ‘ei = 2 Few 
where S. and pay are the height-integrated Pedersen and 


(DEY eer me Lo a ae 
AAS see, ie. shpaioa eee oath 
te Pe stmoeteans gh Hosakon aap rae ol 


er ont aise 
Dos, <> 8 ne 


ete ny ae bons ab Fic 
Tener get. ge Pike atnen mat iele Poae ig at Acre: ‘on tvinddandos: ‘obsongas 
oe Y _ TTT VAyy tm be Py ls hed wasioites oi cal 


he 


iw 


(siete peice al Nl w tee % ‘towbdady yet sive fete 

Me Ne ge Ue Ne ded Bee abs partcondh ee a et at 6 
re dine % neveees Stem aa junvisytae | ‘bite, aL oa Oe 
ty a hi peer Gabd: ade 48 4 “sain yap me onic : , ta" 


7 & 


it ie ‘ y. | Ae rk 
by, LIM: oe | “es “e -e 


tv aovie at  toxxor ae seit 0d0, Kiicupeeatsuueas eae’ 


| ew | Ae ie? oo. a pve € ‘we ; ty S 
Ship's ee 7 Fates te a z “43 z = ae . a fe ia rh ee ; 
re e TiA A Fi a apace 


ae x + os - ae one 


is "ite 
aah Pr ae A a 4 ey uae a) ) ak: Ae ar $ 


Ya: aeexebed:. Nadie eas oie: ene He! salts ‘whede \ oe 


[om . a oe i hia 
10 ey ae ene 


359 


Halli» conductivities respectively. Ex and Gr are the 
east-west and north-south components of the ionospheric 
horizontal electric field, and where northward and eastward 
are taken a positive directions (Anderson and 
Voudra ky, 1975). 

Hipwelt Musa cassumed that the electric field istpurely 


northward, then 


and : 


ATT. 6 


Combining equations AII.1 and AII.2, and substituting for 


a. from equation AII.6 gives: 


G -[% “(G+ a8 hed sb Ly: 
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=| Je) nek as 30)de 7443 


Ew 


2 Zo. 
where J(e) has been written for J, (@) and Gi = G ia =P G; 


# 
is the complete kernel for the problem. This approach 
requires either that Za en pe known, or that 
. latter 
ZplXy may be treated as a free parameter. In the 
the ratio may be adjusted until discrepancies between 


case, 


the model and the observations are minimized. in this study, 


it haat awe ee ‘pa vi on 
TOT SA MS a ane ves 
Dicteist wusnahe Sas pends: 


Wiggiat as Launeiay’ eset re 


: z Ant Bei ie. oo Biel al Pete 
% be!) Seah ans : 
i) 

; 
ars ey 4 ‘ 4 " f | iy 
VR Y Bdliee ey ee att ate, ‘onal ba 

F | 
i ke i ¥ Kies ar Oa a 
eT AL a ty pI Tea eT beet Mn eeyinh TUSWMES SE) tla ee eae 
i teen e 
ao 
- 7 2 + - Vora oy 
fy ” 
3 % ha 4 CoN elo 2 Sa) My ere +4 
Ae Ry ’ “Ae 
| 
ee - y? » ah dee o 
ai i cA 
’ » f me che | eg 
x 
‘ 


= 
4 
- 


eR 6 Ont RP cea Bao tA fren i paseh, ane edu a 


(as | 


eee vay Te et kaos a8, ate 


on 


oa fe . H ne ia A a. a. i 
/ yi * 7 ne ill 
bie), Ohya). Tatas ah ta 
DI ae ie: ae 
1 a r+) 7 
a Ly 7 


vf i eR ih ne: mh) eae Ss ey 
ae re +, 2 Pp. ad ys ania ah od ia ‘os 
ee wi ® mio any : ie uy - Py nan . 


Pps Weck ted Qa Me we en: “quae: a wala ‘aod, oe “erwi 


ie es ha " va ae 
i } 


cro tte + AO, * i 
ny 


ie" ; 
- 


a 


360 


several values of Pie teks Were Used’ initially, canging wilnon 
2.0 to 0.5. The value of 0.5 was found to yield the best- 
fitting results, and this is a reasonable value based on the 
results of Brekke et al (1974). 

The inverse problem that was to be solved can be stated 
assfollows: “given n measurements,B; (J=1,2eeeee0), Qf ivwthe 
perturbation magnetic field, what can be said about Je) @ 
the height-integrated current density, in, sone. cOlatitudinal 
region, 6 <@<@, 2?" Linear inverse theory Shows that only 
linear combinations of the data are available to answer this 
question. That is, 
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cater = & (0.5 = | 
J 6 
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N 
where. 7A(O1,O, )i+= ZX 4;60,)G-(8) and rer erecn yi. is the 
<oh ft v 


6, 
(6) A(6,0,) 46 


AII.8 


estimate of J at 8, ° 

A(G&, @ \ 2s the “averaging funetion' and is 
effectively a window through which the height-integrated 
Current density, is vieweds»Notesthateifp,its pic) epossi ble) eta 


finds ea /aset (not a7zts such that A(Q,@Q,) is a Dirac-deita 


function, then, 


6, | 
Meine yore f T(6) S(0-6,.)\d@ = T() 
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estimate of the height-integrated density in the 


neighbourhood of Oe 


The important feature of this technique is a theorem 
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proved by Backus and Gilbert (1970) Min’ Qulich Tisai] 4 gehown 
that all possible current densities, J(@), which give rise 
to the observations will provide the same estimate, 
S$J{Oo)>, when averaged with a(9,Q@,). As Oldenburg (1976) 
has described, observational e€rrors, or those incurred by a 
deépanturesvof Gthe “real” current systen from the assumed 
forward model, may also be included in the analysis to 
provide estimates of the standard deviation,o@ ( O,)+ of 


<3 ( , )>. If it is assumed that the errors in the data are 


Gaussian distributed with zero mean, Then, 
e ‘sot a0 Zz 
o"/ <4 (8,)9/ = = 9 (0); AII.9 
t= 


where oF is the variance of the By datum. 

Ideally, it ‘is’ desirable to minimize this quantity, 
while at the same time, it is aiso preferable to minimize 
the width of Ree ¢ Cy to give the maximum resolution of the 
height-integrated current density. Backus and Gilbert (1970) 
fave. shown that it is. not. possible to minimize, both these 


quantities simultaneously. It thus becomes necessary to 


formulate the expression 
os , 
G = D(Q)cosp +o (Q,) sin ane 


where D(O,) is the width of a(O,@,) and W is called the 
trade-off parameter. The functional p is then minimized for 
a given value of Y . Thus, for Y =0, the a;'s determined 
from minimizing @ produce a very narrow averaging LiImcir. Ol, 


so the resolution of the estimate <J(Q@,)> is maximized, but 


Binds gh! adit ‘ee ee aa 
bd oo (oti pes mgeLe aes 
wane aeeie - ait, “at bineve 


dene ts el "Me ctidinagal oni 


ae 


i es 


co ge 2 Wo Nt aor DakeaD? Thee tsaive sth * eurtoant 


i begvdanael) DD ie i te “sae detent “ia eae pane 


ere ene Ss Ric oot ‘nidiaargle ati pith) 

he BOs Or iFecenate wiih ea pain on Epo ¥ «BPH Noe ds | 
ae aed - 

OCR Ny Shi ie baie shoe uieh on eae: soteaee bate 


im 


ava Mtonvg ame 2 tan kt ry Ok wanneg ait, ee aR woe awodn ’ MS) oo 


Xs — 
ihe “sb tener any ottlin ot 


Pere ma Ma si Gn = ginal ae MG | 


) i 
Di) 


ate tel rminctind al Shh seoreesinl aona ks 


1 ee Pine AN 
4 val ; 


pat tink Low bee a Bite» ee ae outed “ar | ol: 
i ie Pah A : je 7 vs 
or a hans oa ber ae ak pe nage ieee / 
Uy oe ty Vi A 
isin. tat vr up hh & " Ya til wa a re ; 


ee, 
; Dad 


i iaueioi bebe Lain Whe wg iy pea 
| a Me od , 
die | shots at: a ae ee 


‘, a if - 
‘i i i" . 7 ’ 
A ned oy | Ve gt " -, ( 7 , 
7 ’ i] , v 7 i) ; rf i ; HI Dt . - 


i : 5 wi . - | 
al Wid ; ’ iy es i 7 
q ry, 
7 fe 7 "a i ole 
ve I Wing Pen Arp ne - ‘ 


362 


at the same tine, O” is also maximized. Thus, equation 
ATI.10 expresses a form of uncertainty principle that 
Operates whenever imperfect data is inverted. 

Although there are infinitely many current densities 
that will recreate the observations of the Magnetic field 
(Backus and Gilbert, 1967), the estimate <7(@) > will not 
in general do so. There thus exists the possibility that a 
model current density which does reproduce the data might be 
or some assistance in the interpretation of the results of 
ah inversion. Such models are found by minimizing a specific 
functional, subject to the constraint that the observations - 
be reproduced. 

If the model is constrained further such that the 
current density is zero at the specified latitudinal limits 
Dee the) current 7 O, “and “@ ), then there ista) seteor 


constraints that can be used to define the current density. 


Specifically, from equation AII./7, 


G 
ra J(@) G (e)do j= 543 Le eee 


and 


ECO) eshte Are 


If the so-called "flattest-model" of height-integrated 


current density is chosen as the model that £its, the’ -<deta, 


then the functional 


a) ae ey: AL6 Alters 


rei 8 dager ye afer 


Ve ett) E 


at Or eAae Cert ‘ 


Plaga ‘onan fi 
Kink? ht ae ae pea aiye 


i 
A Pf ee 


Hawk detalles wii PoRidonaaiaeade inte dh 


ae LaAty 


iidn < CS aleale sane eR agudlen. meals 
Pt ht ti? 8 eh sa% ad Aout recbinds ge ange pape. oe om: 


Reva ce 


Pee ae ii sewing ve ere a hbeee bus ie a 


toner tre ea 20 91d at Oe Hort ‘oly ae ‘goadin 


* 


J . Cy coe 


KY hae oR ‘vindie kab Le ul Kioind ae etepow, ‘haa a 


be A rae eg ey ir ahd seas Weu + aml 
4 Pi sR, Sl ese MR ap bhi re bike i no ats an alee Rt 
ined ae 


- 


hashetiobinaa Lae, eel Sana totlalson be ent debinet ana” bal 
Wi RE Rk Be kM SR ee Peron ene wisi og ones Bh obonebale 
BO ena ne od Opener te Ck hws SB) ) ‘heats 
¥ Pca ee nna ee oo vukiinb oP poo al ate dade ad 


spi / 


win bte or math. 


* ‘ ’ 
wre A 
' yo id o ‘ “ 
fi Le 


Rk. PE | ae and ean if A ai Care, 


betes tod» by bath, 19 magpens es fy Lisarde, nt 2: i | 
y ball ‘nl wil rt REAZ nde nas an ng oh: ih, ‘Usteawt bliccaeat i. 


i Leek reba ous toa 


363 


must be minimized with respect to J(@) . 
To minimize equation AII.13 subject to the constraints 
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